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A  pulse  radiolysis  facility  for  the  study  of  fast 
chemical  kinetics  in  aqueous  solution,  based  on  a  Febetron 
706  electron  beam  accelerator,  has  been  established  at  the 
University  of  Florida  Radiation  Chemistry  Laboratory. 

In  previous  work  on  the  radiolysis  of  tetraphenylborate 
(TPB)  solutions  carried  out  in  this  laboratory,  it  was  found 
that  several  organic  products,  including  benzene,  phenol, 
and  biphenyl,  are  produced  with  substantial  yield.   However, 
the  reaction  mechanism  has  not  been  determined.   Since  the 
tetraphenylborate  anion,  TPB~,  is  a  reducing  species,  it 
should  be  more  readily  subject  to  attack  by  OH*  than  by  eaq~ 
or  by  H« .   The  lack  of  reactivity  between  TPB  and  eaq~  has 
been  confirmed  by  directly  monitoring  the  transient  signal 
due  to  ea  ~.   Concerning  the  reaction  with  OH« ,  two  schemes 
were  investigated:  (1)  a  rapid  electron  transfer  from 


B(C6H5)4"  to  0H-;  (2)  OH-  addition  to  B(C6H5)4".   Comparison 
of  transient  absorption  spectra  resulting  from  the  two 
different  schemes  above  suggests  that  the  0H«  addition  is 
the  dominant  reaction  under  conditions  of  N20  saturation, 
with  an  experimentally  determined  second-order  rate  constant 
of  6.2  x  109  M-1s-1.   A  mechanism  based  on  an  initial  first- 
order  self-decomposition  of  the  OH*  adduct,  (C6H5)3BC6H5OHT , 
is  proposed. 

Reactions  initiated  by  OH*  radicals  or  eaq"  in  agueous 
IrCl63~  solutions  were  studied.   The  rate  constant  for  the 
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respective  reactions  were  found  to  be  4.9  x  10"  M  s   and 
6.1  x  109  M~1s"1.   The  oxidation  product,  IrCl62~  disappears 
rapidly  in  N20-saturated  basic  solution  or  in  either  neutral 
N2-saturated  or  aerated  solution,  but  is  nearly  inert  in 
neutral  solution  with  N20  present.   The  complex  IrCl62" 
reacts  rapidly  with  hydrogen  peroxide  in  basic  media,  as 
confirmed  on  the  benchtop  and  by  stopped-flow  kinetics.   It 
is  therefore  inferred  that  reaction  with  H02~  may  account 
for  the  loss  of  IrCl62"  under  basic  conditions.   Loss  of 
Ir(IV)  in  neutral  N2-saturated  solution  without  added  N20 
may  involve  electron  transfer  from  Ir(II),  and  loss  of 
Ir(IV)  in  aerated  solution  is  attributed  to  reduction  by 
superoxide  ion,  02~. 

Kinetic  modeling  on  the  respective  mechanistic  schemes 
of  the  above  systems  gives  good  agreement  with  our  experi- 
mental results. 
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CHAPTER  I 
INTRODUCTION 

1 . 1  Foreword 

The  pulse  radiolysis  technique  offers  a  convenient  way 
of  generating  a  variety  of  unstable  or  transient  species 
under  well-defined  conditions.   This  technique  involves  the 
initiation  of  a  chemical  reaction  process  via  a  short  pulse 
of  ionizing  radiation.   Of  major  interest  is  the  kinetics 
of,  as  well  as  both  the  qualitative  and  quantitative  meas- 
urement of  these  transient  species  involved  in  the  fast 
elementary  steps  of  the  reaction  process. 

The  modification  of  a  gas-phase  pulse  radiolysis  system 
to  a  liquid-phase  pulse  radiolysis  system  was  undertaken  to 
expand  the  chemical  kinetics  facility  of  the  University  of 
Florida  Radiation  Chemistry  Laboratory.   Efforts  to  obtain 
kinetic  information  related  to  various  aqueous  solution 
systems  were  made  in  order  to  supplement  existing  aqueous- 
solution-based  studies  carried  out  in  this  laboratory  in- 
volving techniques  such  as  steady-state  Co-60  7-radiolysis 
and  solar-assisted  photoelectrochemistry. 

To  study  fast  chemical  reaction  systems,  two  necessary 
requirements  are  the  initiation  of  the  process  on  a  time- 
scale  that  is  shorter  than  or  comparable  to  the  process 
being  observed,  and  the  detection  of  a  physicochemical 
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change  associated  with  the  process.   Our  liquid-phase  pulse 
radiolysis  system  was  developed  around  a  Febetron  706  elec- 
tron beam  (e-beam)  accelerator  capable  of  delivering  a 
short,  3  nanosecond  full  width  at  half-max  (FWHM)  pulse  of 
600  KeV  electrons  with  a  peak  current  of  8  x  103  amps.  The 
monitoring  technique  employed  was  time-resolved  absorption 
spectroscopy . 

A  microcomputer-based  data  acquisition  system  was 
integrated  into  the  pulse  radiolysis  system  to  acquire, 
analyze,  and  store  data  obtained  from  the  study  of  fast 
chemical  processes. 

The  kinetic  studies  involving  aqueous  solutions  of 
sodium  tetraphenylborate  (NaTPB)  were  undertaken  to  aid  in 
developing  a  plausible  mechanism  concerning  the  radiolytic 
degradation  of  the  TPB~  ion.   Examination  of  the  iridium 
chloride  system  was  carried  out  to  expand  the  solution 
conditions  under  which  this  system  had  been  previously 
studied  by  pulse  radiolysis.   Due  to  the  similarities  of 
flash  photolysis  and  pulse  radiolysis,  specifically  in  the 
case  of  radiation-induced  redox  reactions,  previous  photoly- 
sis work  on  the  iridium  chlorides  served  as  a  reference  to 
the  pulse  radiolysis  data  collected  in  this  study. 

Various  chemical  dosimeter  methods  applicable  to  pulsed 
electron  sources  were  studied  to  establish  the  proper  dose 
delivered  to  the  liquid  sample  in  the  Febetron  706  system, 
as  modified  for  liquid-phase  work. 


1.2  Review  of  Previous  Work 




1.2.1  Pulse  Radiolysis 


At  the  University  of  Florida  Radiation  Chemistry  Labo- 
ratory, Febetron  706  studies  have  been  completed  on  the  gas 
phase  reaction  kinetics  and  mechanisms  of  various  chemical 
systems  such  as  the  alkyl  and  perf luoroalkyl  iodides, 
oxygen/ozone,2  and  OH  reaction  dynamics  with  small  mole- 
cules.3-5 The  identity  and  reactivity  of  important  tran- 
sient species  generated  in  the  pulse  radiolysis  of  gaseous 
systems,  at  atmospheric  pressures  or  less,  was  of  previous 
interest  in  these  studies. 

Historically,  pulse  radiolysis  was  developed  to  eluci- 
date postulated  transient  behavior  in  radiation-induced 
chemical  reaction  systems,  and  apparatus  was  first  brought 
into  operation  in  about  I960.6'7  Early  pulse  radiolysis 
systems  used  visible  absorption  spectroscopy  to  monitor 
intermediates  produced  from  intense  single  electron  acceler- 
ator pulses  lasting  typically  around  one  microsecond  or 
less,  and  much  of  the  developmental  work  on  pulse  radiolysis 
involved  liquid  phase  systems.8"10  Development  of  the  pulse 
radiolysis  technique  has  since  brought  about  variations  of 
the  method  to  include  stroboscopic  techniques  capable  of 
picosecond  resolution11  and  multiple  time-resolved  observa- 
tion techniques  such  as  electron  spin  resonance  (ESR), 


conductivity, 15~17  emission  spectroscopy,18  and  resonance 
Raman  spectroscopy,19-21  to  monitor  reaction  intermediates. 
The  pulse  radiolysis  method  is  the  radiation  chemistry 
analog  of  flash  photolysis,  in  that  the  photoflash  and 
photodissociation  of  molecules  is  replaced  by  a  short,  high- 
energy  pulse  of  electrons,  causing  excitation  and  ionization 
of  molecules  by  electron  impact.22  The  technique  of  pulse 
radiolysis,  used  to  study  fast  reactions,  has  been  applied 
in  areas  of  aqueous  inorganic,  organic,  and  biochemistry. 
Other  fast  reaction  techniques  such  as  flash  photolysis, 
stopped-flow  techniques,  and  temperature  jump  or  pressure 
jump  relaxation  methods  have  been  designed  to  monitor  con- 
centrations and  to  measure  rate  coefficients.23   However, 

the  pulse  radiolysis  technique  was  the  first  method  capable 
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of  studying  processes  occurring  at  or  below  10   seconds. 

The  characteristics  of  the  electron  pulse  used  in 
electron  beam  pulse  radiolysis  are  its  time  profile,  peak 
current  and  energy,  cross-sectional  homogeneity,  and  diver- 
gence.  For  adequate  time  resolution  of  the  observed  tran- 
sient signal,  the  time  profile  of  the  pulse  should  be  very 
short  compared  to  the  relaxation  of  the  intermediates.   The 
pulse  current  determines  the  number  of  interactions  between 
the  chemical  system  and  the  e-beam  which,  in  turn,  deter- 
mines the  concentration  or  amount  of  species  formed.   The 
maximum  electron  energy  determines  the  penetration  depth  of 
the  e-beam  into  the  reaction  cell  and  should  be  high  enough 
to  adequately  irradiate  the  volume  of  the  chemical  system 


being  monitored.   The  cross-sectional  homogeneity  and  diver- 
gence of  the  beam  determine  how  uniformly  irradiated  the 
optically-sampled  reaction  region  is. 

Transformations  over  the  years  of  the  standard  pulse 
radiolysis  installation  have  led  to  increasingly  sophisti- 
cated experimental  setups.25-28   Spectrographic  methods 
involving  oscilloscopes  and  Polaroid  cameras  have  been 
largely  replaced  by  digital  signal  processors,  or  transient 
digitizers  that  have  an  interface  for  computer  control.1 
Signal  processing  has  been  aided  by  the  use  of  backing-off 
circuits  which  can  measure  small  transients  on  top  of  the 
analyzing  light.29   Lasers  and  the  capability  of  commercial 
xenon-arc  lamps  to  supply  pulse  current  on  top  of  DC  current 
have  provided  intense  analyzing  light  sources.30"32   Lastly, 
on-line  computer  systems  have  been  developed  which  not  only 
greatly  reduce  the  time  reguired  for  data  analysis  but  also 
improve  both  guantity  and  guality  of  the  data 
processed. 29 ' 33~35 

1.2.2  Radiolvsis  of  Water  and  Aqueous  Solutions 

Early  investigations  on  the  radiolysis  of  water  have 
played  a  significant  role  in  the  development  of  radiation 
chemistry.36-39  In  addition,  sophisticated  pulse  radiolysis 
technigues  achieving  the  time  resolution  of  radiation  in- 
duced chemical  events  have  been  applied  more  often  in  the 
radiolysis  of  water  than  in  any  other  liquid.40   As  a 


result,  most  of  the  aspects  of  the  radiation  chemistry  of 
water  and  aqueous  solutions  are  now  reasonably  well  under- 
stood. 

A  few  practical  purposes  for  water  irradiation  studies 
arise  from  1)  the  importance  of  understanding  the  effect  of 
radiation  on  biological  systems  in  which  water  is  the  main 
substrate,  and  2)  various  aspects  of  nuclear  reactor  tech- 
nology concerning  water's  use  as  both  a  coolant  and  modera- 
tor.41'42 Water  and  aqueous  solutions  have  been  studied,  in 
particular,  for  three  reasons:  (1)  their  role  as  solvents  is 
central  to  chemistry  in  general  and  to  radiochemistry  in 
particular;  (2)  they  are  readily  available  and  easy  to  work 
with;  (3)  water  is  a  polar  liquid  that  responds  in  charac- 
teristic ways  to  radiation. 

Important  features  of  the  radiation  chemistry  of  water 
are:  (1)  a  description  of  the  primary  radical  and  molecular 
products  and  how  they  evolve  in  both  time  and  space;  (2)  the 
effect  of  pH  on  the  radiation  chemical  yields  of  primary 
products;  (3)  the  track  structure  and  related  phenomena 
associated  with  the  passage  of  a  charged  particle  through  a 
liquid;  and  (4)  a  detailed  account  of  all  the  radical  reac- 
tions involving  the  primary  radicals  and  radiolysis 
products.43 

Upon  initial  ionization  of  a  water  molecule  by  a  fast 
primary  charged  particle  such  as  an  electron,  the  liberated 
electron  from  the  ionization  event  will  typically  have 
sufficient  energy  to  further  ionize  several  adjacent  water 


molecules,  which  leads  to  the  creation  of  clusters  of  ions 
along  the  track  of  the  ionizing  particle.   These  clusters 
are  called  spurs;  detailed  descriptions  of  their  structure 
and  distribution  along  the  particle  tracks  can  be  found  in 
the  literature.44  Once  the  dry  electrons  have  been  slowed 
to  thermal  energy,  they  undergo  a  hydration  (or  solvation) 
process  which  involves  the  water  molecules  becoming  oriented 
about  the  charged  species.   This  process  occurs  on  the  order 
of  10-11  seconds,  which  is  the  relaxation  time  for  dipoles 
in  water.   Alternatively,  some  water  molecules  are  excited 
to  upper  electronic  states  from  which  they  can  autoionize, 
dissociate,  or  fall  back  to  the  ground  electronic  state. 
Figure  1  from  Buxton43  shows  the  time  scale  of  events  initi- 
ated by  the  absorption  of  energy  by  water  from  an  incident 
ionizing  radiation,  such  as  an  energetic  photon  or  charged 
particle.   Also  shown  are  the  various  routes  to  formation  of 
the  primary  radical  and  molecular  products. 

At  approximately  10~10  seconds  after  the  ionizing 
event,  all  the  physicochemical  processes  are  complete  and 
the  chemical  stage  begins.   Initially,  this  process  involves 
diffusion  of  the  radiolysis  products  out  of  the  spurs  where 
they  either  react  together  to  form  molecular  and  secondary 
radical  products,  or  they  escape  to  the  bulk  solution, 
becoming  homogeneously  distributed.   Appendix  B  lists  the 
elementary  reaction  steps  in  the  radiolysis  of  liguid  water, 
along  with  their  appropriate  rate  constants. 
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Fig.  1.   Scheme  for  the  radiolysis  of  water, 


At  approximately  10"7  seconds,  the  spur  expansion  is 
complete,  resulting  in  the  so-called  "primary  yield"  of  the 
radiolysis  products  eaq" ,  H30+,  OH,  H,  H2 ,  H202 ,  and  H02 . 
Radiation  chemical  yields  are  expressed  as  G-values,  as 
shown  in  Table  l,38  where  G  is  the  number  of  species  (repre- 
sented by  X)  created  or  destroyed  per  100  eV  of  absorbed 
dose,  and  are  written  as  G(X). 

Table  1 
Radical  and  Molecular  Product  Yields  in  Irradiated  Water 


Primary  yield     eaq"  H+    OH    H    H2   H202    H02 

G(  molecules  ^    263      2<63   2.72   .55   .45    .68      0.026 
100  eV 


The  yields  are  essentially  independent  of  pH  except  for 
extreme  acidic  conditions  (pH  <  3).   The  G-values  are  also 
dependent  on  the  quality  of  the  ionizing  radiation,  or  the 
mean  linear  energy  transfer  (LET).   Under  extremely  intense 
irradiations,  at  dose  rates  so  high  that  the  spurs  overlap, 
yields  may  also  depend  on  the  dose  rate.40   In  the  case  of 
high-LET  radiation415  such  as  a   particles,  heavy  ions  or 

g 

fission  fragments,  or  very  high  dose  rates  (  >  10 
rad/sec) ,46-49  lower  radical  and  larger  molecular  yields 
result. 

In  the  radiolysis  of  dilute  solutions,  most  of  the 
energy  is  absorbed  by  the  solvent,  and  effects  due  to  energy 
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deposited  directly  in  the  solute  are  unimportant  for  solute 
concentrations  less  than  about  0.01  M.   Thus  all  the  ob- 
served chemical  changes  in  irradiated  dilute  aqueous  solu- 
tions are  brought  about  indirectly  by  the  molecular  and 
radical  products  of  the  water  radiolysis;  this  phenomenon  is 
referred  to  as  the  "indirect  effect."40   In  the  presence  of 
a  solute,  transient  species  and  final  products  result  from 
attack  of  H,  ea  ~  OH,  and  H202  on  the  added  substrate. 
Radicals  derived  from  the  substrate  may  also  react  to  form 
final  products.   As  evinced  by  the  G-values  of  the  primary 
radicals,  the  radiolysis  of  water  produces  approximately 
equal  numbers  of  powerful  oxidizing  and  reducing  radicals. 
The  standard  reduction  potentials  for  the  hydroxyl  radical, 
OH,  and  the  hydrated  electron,  eaq" ,  are  2.72  V  and  -2.9  V, 
respectively . 40 

For  pulse  radiolysis  chemical  applications,  well-de- 
fined conditions,  i.e.  either  totally  oxidizing  or  totally 
reducing,  are  desirable,  and  these  conditions  can  be 
achieved  by  interconversion  of  the  primary  radicals. 
Briefly,  for  a  convenient  way  of  obtaining  almost  totally 
oxidizing  conditions,  N20-saturated  (N20  =  2.5  x  10"2  M) 
solutions  are  used.   N20  converts  e  ~  to  OH,  or 

e~  +  N20 >  N2  +  0" >  OH  +  0H~        129 

Acidic  conditions  of  pH  <  3  should  be  avoided  due  to  the 
interconversion  of  eag~  to  the  hydrogen  atom,  which  is  a 
major  reducing  species  in  acidic  media, 

*aq 


e   ~  +  H+ >  H     103 
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For  the  creation  of  almost  totally  reducing  conditions, 
addition  of  an  organic  solute  to  an  inert  gas-saturated 
solution  leaves  ea  ~  as  the  predominant  reducing  primary 
species.   Reactions  involved  are 


OH  +  ROH >  R(CH.)OH  +  H20       1-1 

H  +  ROH >  R(CH.)OH   +  H2        1-2 

In  general,  the  organic  radicals  formed,  R(CH.)OH  ,  are  less 
powerful  reductants  than  ea  ~  and  are  rather  unreactive.   In 
particular,  N2-saturated  solutions  with  added  tertiary 
butanol  were  used  in  our  work. 

1.2.3  Aqueous  Solution  Studies  of  Aromatic  Systems 

In  aqueous  systems  containing  organic  solutes,  reac- 
tions occurring  between  solute  and  the  primary  products  of 
water  radiolysis  are  predominately  abstraction,  addition, 
and  one-electron  reduction.   Organic  radicals  produced  by 
these  processes  can  further  react  before  being  converted  to 
stable  products,  thus  leading  to  complex  radiolysis  mecha- 
nisms,  which  have  been  reviewed.  -L_° 

Mechanisms  can  be  somewhat  generalized  within  families 
of  organic  compounds  containing  the  same  functional  group 
due  to  the  tendency  of  the  primary  products,  eaq~ ,  H  and  OH, 
to  react  with  these  functional  groups  as  opposed  to  the 
molecule  as  a  whole.40   Particularly,  benzene  can  be  uti- 
lized as  an  analogous  system  for  many  aromatic  compounds, 
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with  addition  to  the  unsaturated  bonds  being  the  character- 
istic reaction  path.   Reaction  of  the  primary  radicals  from 
water  radiolysis  with  aromatic  systems  exhibit  rate  con- 
stants in  the  range  of  105  to  107  M~1s~1  for  eaq~ ,  and  >  109 
M-1s-1  for  H  and  OH.40 

C6H6  +  eaq"  —->    C6H6~  (+  H2°)  — >  C6HV  +  0H"     M 

C6H6  +  H  " >  C6HV  I_4 

C6H6  +  OH  >  C6H6OH.  1-5 

Products  of  these  addition  reactions  are  cyclohexadiene 
systems  containing  an  unpaired  electron  delocalized  about 
the  benzene  ring.   The  dienyl  radicals,  namely  cyclohexadi- 
ene, C6H7«,  and  hydroxycyclohexadiene ,  C6H6OH.,  have  been 
investigated  by  electron  spin  resonance  techniques  that 
indicate  a  positive  electron  density  assignment  to  the 

54 

ortho-  and  para-carbon  atoms. 

Extensive  pulse  radiolysis  work  has  shown  these  dienyls 
exhibit  a  reasonably  strong  absorption  {e   *  5  x  10  M  ■'cm   ) 
in  the  310-315  nm  range.55'56   Rate  constants  involving  the 
primary  radicals  with  benzene  as  well  as  the  dimerization 
and  disproportionation  of  the  resulting  cyclohexa-  and 
hydroxycyclohexadienyls  have  also  been  reported. 

Steady-state  7-radiolysis  studies  involving  liquid 
chromatography,62  radio-liquid  chromatography,63  and  gas 
chromatography,64'65  as  analyzing  methods,  indicate  a  com- 
plex mixture  of  products  formed  from  the  combination  and 
disproportionation  of  the  cyclohexa-  and  hydroxycyclohexa- 
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dienyl  radicals.   Phenol  and  biphenyl  are  the  major  products 
reported,  in  addition  to  mixed  dienyls  and  regenerated 
benzene.   All  of  the  previous  studies  were  carried  out  in 
N20-saturated  or  deoxygenated  aqueous  benzene  solutions.   In 
aerated  aqueous  benzene  solutions,  ea  ~  and  H  are  scavenged 
by  02 •   The  hydroxycyclohexadiene  radical  reacts  with  02 
forming  phenol  and  phenol-like  products.66'67 

1.2.4  Sodium  Tetraphenylborate  System 

Pulse  radiolysis  studies  of  the  aqueous  NaTPB  solutions 
have  been  carried  out  in  this  laboratory  to  help  establish  a 
mechanism  consistent  with  the  end  product  yields  obtained 
from  Co-60  7-radiolysis  on  aqueous  NaTPB  solutions.68  These 
steady-state  radiolysis  studies  were  done  in  conjunction 
with  the  United  States  Department  of  Energy  operation  at  the 
Savannah  River  Plant  in  South  Carolina  and  concerned  the 
pi  oblem  of  cesium-137  separation  and  removal  from  nuclear 
waste  obtained  from  spent  fuel  reprocessing.69   Precipita- 
tion of  the  radiocesium  as  its  insoluble  tetraphenylborate 
salt  is  a  method  of  separating  it  from  the  nonradioactive 
components  and  concentrating  it  to  reduce  the  volume  of  the 
liquid  radioactive  waste.   Other  decontamination  methods 
include  ion  exchange  and  evaporation  techniques.70'71   The 
radioactive  precipitate  will  be  stored  in  large  underground 
steel  tanks  until  it  is  sent  to  the  Defense  Waste  Processing 
Facility  (DWPF)  for  immobilization  in  glass.72'73 
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The  TPB  precipitate  will  undergo  intense  p-   and  ir- 
radiation from  the  cesium-137  during  processing  and  storage 
that  will  partially  decompose  the  tetraphenylborate,  produc- 
ing several  organic  compounds.   Understanding  the  reaction 
scheme  which  leads  to  these  various  organic  products  as  well 
as  their  impact  on  subsequent  processing  of  the  precipitate 
are  important  aspects  concerning  the  design  and  operation  of 

the  DWPF.68 

Since  the  tetraphenylborate  anion,  TPB" ,  is  rather 
different  structurally  from  benzene  or  simple  aromatics,  it 
is  not  obvious  that  previous  work  on  the  mechanism  of  ben- 
zene radiolysis  is  relevant.   However,  several  lines  of 
evidence  will  be  presented  in  this  dissertation  which  indi- 
cate that  there  is  a  substantial  connection.   The  prepara- 
tion and  properties  of  aqueous  metal  tetraphenylborate  salts 
have  been  described  by  Flaschka  and  Barnard.74   LiTPB  and 
NaTPB  are  found  to  be  the  only  TPB  salts  that  are  apprecia- 
bly soluble  in  water;  all  TPB  salts  tend  to  be  soluble  in 
polar  organic  solvents. 

Only  a  few  mechanistic  studies  involving  primary  radi- 
olysis radical  reactions  with  the  TPB"  anion  have  been  found 
in  the  literature.   Horii  and  Taniguchi ,  using  pulse  radiol- 
ysis, studied  the  oxidation  intermediates  of  the  TPB"  anion 
in  aqueous  solutions  using  the  oxidizing  radicals  N3 ■ ,  Br2 * , 
and  SCN2T.75  Kinetic  as  well  as  spectral  information  was 
reported  lor  the  resulting  intermediate,  B(C6H5)4',  which 
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was  produced  in  all  cases  by  an  electron  transfer  from  the 
TPB-  anion  to  the  oxidizing  radical.   Liu  et  al.  reported  an 
absorbing  species  at  A  =  300  nm  for  the  expected  radical 
derived  by  OH*  addition  to  the  phenyl  group  of  the  TPB- 
anion.76   Another  transient,  derived  from  H'  addition  to  the 
phenyl  group  was  also  discussed;  however,  no  kinetic  infor- 
mation was  supplied  regarding  the  fate  of  these  resulting 
transient  species. 

Other  pulse  radiolysis  investigations  involving  the 
NaTPB  salt  include  ion  association  studies,  in  which  the 
spectra  of  the  diamagnetic  species,  Na" ,  the  solvated  elec- 
tron, e  ",  and  the  ion  pair  species  (Na+,eaq")  are  reported 
for  alkali  metal  solutions  in  various  non-aqueous 
solvents.77-82 

The  photochemistry  literature  of  the  metal  tetraphenyl- 
borates  was  surveyed  to  find  possible  mechanisms  concerning 
the  reactions  of  the  TPB"  anion  and  subsequent  intermediates 
under  UV-irradiation.   Photolysis  of  aqueous  solutions  of 
NaTPB  using  UV  light  with  a  wavelength  of  2537  A  generates 
biphenyl  as  the  major  product  in  oxygenated  solutions  and  1- 
phenyl-l,4-cyclohexadiene  as  the  major  product  in  the  ab- 
sence of  02.83'84  Other  photochemical  studies  involve  the 
photolysis  of  tetraphenylborates  in  non-aqueous  organic 
solvents.85-87 

Additional  photochemical  studies  involving  TPB"  report 
on  photochemically-induced  electron  transfer  from  the  TPB" 
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anion  to  oxygen,88  and  various  intermolecular  charge-trans- 
fer transitions  of  donor-acceptor  ion  pairs  involving  the 
TPB"  anion. 89-91 

1.2.5  Agueous  Solution  Studies  of  Transition-metal  Complexes 

In  agueous  systems  containing  transition  metal  complex- 
es, the  predominant  reaction  occurring  between  the  solute 
and  the  primary  products  of  water  radiolysis  is  a  one-elec- 
tron transfer,  either  oxidation  or  reduction.   This  electron 
transfer  can  occur  by  an  inner-  or  outer-sphere 
mechanism,92'93  often  leaving  the  complex  ion  in  an  unstable 
oxidation  state.   Agueous  solutions  of  both  transition 
metals  and  their  complexes,  most  notably  bio-inorganic 
molecules,  have  been  irradiated  and  the  radiation  chemical 
processes  which  occur  have  been  extensively  reviewed.   ■ 

In  agueous  systems,  aguation  of  the  metal  complex  from 
ligand  exchange  with  water  molecules  can  occur,  and  has  been 
reported  in  both  pulse  radiolysis  studies96-99  and  photo- 
chemical studies.100'101   Changes  in  the  stereochemistry 
and  coordination  number  of  the  transition  metal  complex  upon 
pulse  irradiation  have  also  been  investigated.102'103 

1.2.6  Hexachloroiridate  System 

The  redox  chemistry  of  the  IrCl62-  and  IrCl63-  complex- 
es has  been  of  recent  interest  due  to  the  prospect  of  using 
the  system  to  produce  molecular  hydrogen  by  the  solar  photo- 
chemical reduction  of  protons  in  acidic  agueous  solu- 
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tions.104  The  solar  chemistry  of  metal  complexes  has  been 
reported  by  Gray  and  Maverick.105  Solar  chemistry  is  de- 
fined as  that  area  of  photochemistry  in  which  the  excitation 
energies  fall  within  the  spectrum  of  solar  irradiation  at 
the  earth's  surface.   The  proposed  scheme  involving  agueous 
HC1  solutions  of  IrCl62"/3-  is 

IrCl63-  +  H+  >  IrCl62"  +  1/2H2        1-6 

IrCl62"  +  CI"  >  IrCl63"  +  1/2C12      1-7 

The  net  result  is  the  conversion  of  HC1  to  hydrogen  and 
chlorine  using  UV  and  visible  light  as  the  only  energy 
sources.   While  Gray  reports  that  the  second  reaction  occurs 
with  low  efficiency  at  wavelengths  higher  than  430  nm,  the 
first  unfortunately  occurs  only  with  UV  light,  e.g.  at  254 
nm. 

In  this  study  we  attempt  to  further  elucidate  the 
radiation  chemistry  of  the  chloroiridates  under  various 
agueous  solution  conditions,  namely  in  alkaline,  N20- satu- 
rated solutions  and  in  neutral,  aerated  solutions. 

2—  3  — 

Dynamic  processes  occurring  in  the  IrClg   /  IrCl6 

system  have  been  previously  investigated  by  photochemical 
technigues.100'101'106'107   Dainton  and  Rumfeldt  made  perti- 
nent observations  concerning  the  chemistry  of  the  system, 
although  their  experimental  procedure  involved  Co-60  7- 
radiolysis  rather  than  direct  studies  of  the  transients.108 
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Mills  and  Henglein  used  steady-state  gamma  radiolysis 
of  aqueous  IrCl63~  to  explore  the  formation  of  colloidal 
iridium  and  the  catalytic  properties  of  the  colloid  with 
respect  to  H2  formation.109'110   Pulse  radiolysis  of  NaIrCl6 
aqueous  solutions,  carried  out  by  Mills  and  Henglein,  exam- 
ined the  reduction  of  IrCl63~  by  eaq"  and  the  (CH3)2COH- 
radical.   A  fundamental  pulse  radiolysis  study  on  the  aque- 
ous solution  redox  reactions  of  hexachloroiridate  complexes 
by  Broszkiewicz  addressed  both  the  reduction  and  oxidation 
of  aqueous  IrClg3-.111  From  his  work,  the  author  concluded 
that  octahedral  IrCl63~  is  oxidized  by  a  single  electron 
transfer  to  octahedral  IrCl62-  without  changing  its  struc- 
ture or  ligand  composition,  i.e.  aquation.   Edwards  also 
commented  that  the  iridium  chlorides  are  typical  examples  of 
complex  ions  that  are  relatively  inert  toward  substitution, 
but  which  undergo  rapid  electron  transfer  reactions.112 

Further  pulse  radiolysis  studies  on  transition  metal 
hexachloride  complexes  have  been  reported  by 
Broszkiewicz.113-115  Other  pulse  radiolysis  studies  involv- 
ing aqueous  IrCl63~  have  examined  its  oxidation  by  the 
azidyl  radical,  N3-,116  and  the  hydroxyl  radical,  OH'.93 


CHAPTER  II 
EXPERIMENTAL  APPARATUS  AND  PROCEDURES 

2.1   Pulse  Radiolysis  Instrumentation 


Reactive  intermediates  generated  in  radiation-induced 
chemical  processes  are  short-lived;  thus  the  pulse  radioly- 
sis instrumentation  must  be  able  to  monitor  the  transient 
behavior  with  a  high  degree  of  resolution  over  extremely 
short  intervals.   Since  the  concentration  of  such  intermedi- 
ates can  influence  both  the  mechanistic  pathway  and  the 
final  product  yields,  detection  techniques  that  are  sensi- 
tive to  very  low  concentrations  of  radicals  are  necessary 
for  a  complete  characterization  of  the  transient  behavior. 

The  University  of  Florida  pulse  radiolysis  system 
configuration  incorporating  computer  control  is  shown  in 
Figure  2.   The  accelerator  used  is  a  Febetron  (brand  name 
for  the  impulse  generator  produced  by  Field  Emission  Corp.) 
706,  which  operates  on  the  principle  of  the  Marx-bank  cir- 
cuit.  Other  accelerators  used  in  time-resolved  studies 
include  the  linear  accelerator  (linac),  and  the  Van  de 
Graaff  generator.   One  may  see  that  the  other  components 
shown  comprise  an  absorption  spectrophotometer.   Absorption 
spectroscopy  is  the  most  widely  used  technique  among  others 
mentioned  previously,  for  monitoring  transient  behavior  in 
pulse  radiolysis. 
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2.2  Absorption  Spectrophotometry 

Components  involved  in  an  absorption  spectrophotometer 
are  the  light  source,  an  optical  train  including  all  lenses, 
mirrors  and  absorption  cell,  the  monochromator ,  the  detector 
and  a  recorder.   In  flash  photolysis  and  pulse  radiolysis, 
the  absorption  cell  also  serves  as  the  reaction  region. 

2.2.1  Light  Source 

The  light  source  used  should  provide  high  monitoring 
light  output  to  optimize  the  signal-to-noise  ratio,  S/N, 
which  is  directly  related  to  the  intensity  of  light  striking 
the  photocathode  of  a  photomultiplier  by  the  eguation: 

S/N  =  (Is/  A  f)1/2         II-l 

where  Is  is  the  light  intensity  at  the  photocathode  and 
delta  f  is  the  bandwidth  of  the  recorder  system.117'118 
This  relationship  is  applicable  when  signal  shot  noise  is 
the  limiting  S/N  determining  factor  as  is  the  case  in  fast 
(  >  1  kHz)  kinetic  spectrophotometry.   Lasers  provide  high 
monitoring  light  intensities  but  are  limited  as  to  their 
wavelength  of  operation.   Thus  in  order  to  utilize  an  ab- 
sorption technigue  that  is  applicable  to  all  transients  or 
combination  of  transients  throughout  the  ultraviolet  and 
visible  spectrum,  a  high  power  spectral  lamp  with  a  continu- 
ous distribution  is  preferred.   Xenon-arc  lamps  filled  with 
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xenon  gas  at  above  atmospheric  pressures  provide  the  follow- 
ing features:  high  arc  stability  and  luminous  flux,  high 
reliability  and  long  operating  lifetimes  up  to  1000  hours. 
In  addition,  high-current  pulses  can  be  delivered  to  the  arc 
lamp,  which  enable  very  high  monitoring  lamp  outputs  during 
transient  measurement  times  of  microseconds  to  milliseconds. 

2.2.2  Absorption  Path  and  Monochromators 

The  flow  reaction  cell  previously  mentioned  allows  for 
the  production  of  transient  species  in  the  irradiated  volume 
which  can  be  directly  monitored  by  absorption  spectrophotom- 
etry.  For  example,  using  the  Beer-Lambert  relationship, 

A  =  €  b  c  H-2 

—  i   —l 
where     e  =   molar  extinction  coefficient  (M  cm   ) 

b    =  pathlength  (cm) 

c    =  concentration  (M) 
and  assuming  1)  a  transient  molar  absorption  coefficient,  e, 
of  approximately  1000  mol~1dm3cm~1,  2)  a  transient  popula- 
tion which  produces  a  10%  deflection  in  the  initial  light 
intensity,  and  3)  a  1  cm  absorption  pathlength,  one  calcu- 
lates the  initial  transient  concentration  formed  in  the 
irradiated  volume  immediately  following  the  pulse  to  be 

—  5 
approximately   4.5  x  10  ~"   M. 

Monochromators  function  to  isolate  a  small  wavelength 

band  from  a  polychromatic  source  and  consist  of  a  dispersive 

element,  either  a  prism  or  grating,  and  an  image  transfer 
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system  of  slits  and  mirrors.   By  using  a  monochromator  with 
an  adequate  light  flux  throughput  and  a  varying  wavelength 
range,  transient  signals  can  be  recorded  as  a  function  of 
wavelength  to  build  a  time-resolved  transient  spectrum. 
The  throughput  factors  which  determine  the  radiant 
power  passing  out  of  the  exit  slit  of  a  monochromator  are 

—  2    —1    —1 

the  source  spectral  radiance,  BA  in  units  of  Wmm  ■  sr   nm 
and  the  monochromator  throughput  factor,  Y  in  units  of  mm 
sr  nm.   The  product  of  these  two  factors  gives  the  output 
spectral  radiant  power,  Q,  in  watts.118   For  operating  a 
monochromator  with  a  broadband  or  continuum  source  the 
monochromator  throughput  factor,  Y,  is  defined  as 

(*/4) 

y  =  w2  h  T  Rd       II-3 

(F/n) 

where     w    =  slit  width  (mm) 
h    =  slit  height  (mm) 
F/n  =  f-number  or  effective  aperture 
T__  =  optical  transmission  factor  (unitless) 
Rd   =  reciprocal  linear  dispersion  (nm/mm) 
Since  the  light  image  is  usually  focused  down  on  the  en- 
trance slits  by  a  condensing  lens,  the  slit  width  is  the 
only  experimental  variable  parameter  in  determining  the 
light  throughput  at  a  particular  wavelength  for  a  monochro- 
mator with  a  fixed  optical  transmission  factor,  effective 
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aperture   and  reciprocal    linear  dispersion.    *' 
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2.2.3  Photomultipliers 

The  advantages  of  using  photomultiplier  detectors  to 
display  transient  absorption  signals  on  transient  recorders 
are  high  sensitivity  due  to  secondary  emission  multiplica- 
tion, fast  response  times  relative  to  the  half -life  of  the 
relaxation  processes  being  monitored,  and  the  ability  to 
handle  high  light  flux  levels  while  maintaining  a  linear 
signal  at  high  output  photocur rents.119  Two  types  of  photo- 
multipliers  were  used  in  this  work  which  have  different 
spectral  response  ranges.   An  EMI  Model  9250  QB  end-on  50  mm 
window  PMT  fitted  in  a  Pacific  Precision  Instruments  Model 
3262  housing  was  used  initially  with  the  geometry  shown  in 
Figure  3.   Because  of  its  extended  spectral  response  range 
above  600  nm,  a  Hamamatsu  Model  R928  side-on  30  mm  window 
PMT  fit  into  a  Pacific  Precision  Instruments  Model  3150 
housing  was  eventually  incorporated  into  the  geometry  shown 
in  Figure  4 . 

The  photomultipliers  were  powered  by  a  John  Fluke,  Inc. 
Model  412A  High  Voltage  DC  power  supply,  which  proved  to  be 
a  well  regulated  and  stable  dc  source.   The  overall  gain  on 
the  PMT  could  be  altered  by  changing  the  applied  voltage, 
typically  from  1000  Vdc  to  1250  Vdc  in  increments  of  10  or 
100  volts,  while  maintaining  maximum  gain  on  the  operational 
amplifier. 

The  EMI  end-on  tube  has  a  bialkali  photocathode ,  a 
Spectrosil  guartz  window,  and  an  optical  range  of  165-650  nm 
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with  a  peak  wavelength,  or  maximum  radiant  sensitivity 
(mA/W) ,  at  4  20  nm.   Figure  5  shows  a  diagram  of  the  base 
components  which  were  customized  to  achieve  optimum  rise 
times,  gains  and  linearities  for  use  in  the  pulse  radiolysis 
system.   The  EMI  Model  GB  1001A  fast  electronic  switching 
circuit  connected  through  an  optoisolator  and  switching 
transistor  gates  off  the  first  dynode  voltage  of  the  PMT  to 
eliminate  overloading  problems  caused  by  short  burst  expo- 
sure to  high  light  intensities.   The  gating  circuit  offers  a 
precise  way  to  regulate  the  PMT  from  a  +  5  volt  gate  signal 
and  has  the  specifications  of  an  on-to-off  fall-time  of  0.5 
microseconds  and  an  off-to-on  rise-time  of  1.5  microseconds. 
This  circuit  which  is  normally  applied  in  flashlamp  or  laser 
systems  to  protect  from  overloading  conditions  was  original- 
ly employed  to  eliminate  the  overloading  of  the  PMT  by 
fluorescence  of  OH  radicals  which  occurred  in  previous  gas 
phase  studies  on  the  argon  sensitized  pulse  radiolysis  of 
water  vapor.3   In  the  present  liguid  phase  work,  such  severe 
fluorescence  problems  derived  from  radiation  interactions 
with  either  the  agueous  chemical  system  or  the  guartz  cell 
have  not  been  encountered. 

Gating  of  the  PMT  was  applied  to  reduce  the  effect  of 
high  levels  of  electromagnetic  interferences  (EMI)  associat- 
ed with  the  high  beam  current  Febetron  706  machine.   Ac- 
counts of  these  electromagnetic  fields  and  their  effects  on 

"\  s~     ion 

detection  electronics  are  discussed  in  the  literature.^  ' 
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By  gating  off  the  initial  dynode  of  the  PMT  during  the 
radiation  pulse,  overloading  of  the  photocathode  by  the  X- 
ray  (X-radiation)  and  the  radio-frequency  (RF)  noise  coinci- 
dent with  the  electron  pulse  was  significantly  attenuated. 
Further  attempts  to  reduce  the  EMI  effects  included  the  use 
of  lead  shields  inserted  into  the  entrance  of  the  PMT  hous- 
ings, Figure  6,  and  extensive  lead  shielding  around  the 
outside  of  the  PMT  housing. 

The  Hamamatsu  side-on  PMT  was  obtained  to  expand  the 
wavelength  range  to  allow  transient  absorption  measurements 
above  600  nm,  and  was  used  in  all  subsequent  experiments 
after  installation.   This  PMT  has  a  multialkali  photocath- 
ode, a  UV-glass  window,  and  an  optical  range  of  195-900  nm 
with  a  peak  wavelength  at  400  nm.   Figure  7  shows  a  diagram 
of  the  customized  base  components  used  with  this  PMT  along 
with  a  gating  circuit  that  is  analogous  to  the  one  previous- 
ly described. 

2 . 3   Light  Source  Systems 

For  a  versatile  pulse  radiolysis  system  with  capabili- 
ties of  nanosecond  sampling  times  in  the  UV-VIS  region  of 
the  spectrum,  stable  and  high  analyzing  light  intensities 
are  needed  in  order  to  permit  the  observation  of  weak  or 
long-lived  absorption  signals,  to  reduce  shot  noise,  and  to 
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increase  the  signal  to  noise  ratio.  " 
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a)  2  1/16  inch  diameter  x  1/2  inch  thick  lead  shield 
for  end-on  PMT  housing,  Model  3262. 


/ 


! 


' 


\ 


b)  1  1/4  inch  diameter  x  1/2  inch  thick  lead  shield 
for  side  window  PMT  housing,  Model  31 50. 


Fig.  6.   Lead  shields  for  PMT  photocathodes  in  the  Pacific 
Precision  Instruments  PMT  housings. 
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2.3.1  Xenon-Arc  Lamp 

The  lamp  system  used  in  this  work  consisted  of  a  xenon- 
arc  high  pressure  discharge  lamp  contained  in  a  water-cooled 
housing  and  powered  by  a  water-cooled  ignition  device  /  DC 
power  supply.   Several  different  bulbs  including  an  Advance 
Radiation  Corporation  XSA  350  and  500  watt  and  a  Ushio  UXL 
302-0  300  watt  were  used  throughout  the  work.   The  ARC  bulbs 
were  of  the  deep  UV  design,  in  which  a  metal  vapor,  such  as 
mercury,  is  present  in  addition  to  the  xenon  gas  providing  a 
strong  distribution  in  the  deep  UV  region  (190  to  250  nm) . 
The  Ushio  bulb  was  not  seeded  with  such  a  metal  and  its 
continuous  spectrum  is  shown  in  Figure  8.   This  spectrum  was 
generated  by  feeding  the  PMT  output  across  a  load  resistor 
and  into  a  Linear  Instruments  Corporation  Model  591  chart 
recorder . 

The  bulb  was  contained  in  a  PRA-Canada  ALH  220  housing 
with  dual  water  cooling  loops  and  powered  by  a  PRA-Canada 
Model  M301  ignition-power  supply.   The  power  supply  uses  a 
continuous  feedback,  constant  current  mode  of  operation  to 
give  excellent  current  stability.   Current  programming  with 
sguare,  triangular,  or  sine  wave  modulation  can  be  used  in 
this  unit  to  vary  the  optical  output  of  the  lamp.   As  shown 
in  Figure  9,  a  2-3  L/minute  cooling  line  was  in  series  with 
the  power  supply  and  top  face  plate  portion  of  the  housing 
which  was  oriented  vertically.   A  separate  0.25  L/minute 
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Fig.  9.   Xenon-arc  lamp  schematic  diagram. 
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cooling  line  was  used  for  cooling  of  the  lamp  electrodes. 
The  approximate  dimensions  of  the  different  optical  arrange- 
ments are  also  shown  for  the  lamp  housing  which  contains  an 
f/2.4  exit  lens. 

For  conditions  of  1)  low  absorbance  monitoring  on  the 
order  of  a  1%  change  in  light  intensity  or  less,  or  2)  lower 
relative  light  output  (i.e.  the  lower  and  upper  regions  of 
the  lamp  output  spectrum  shown  in  Figure  8)  a  short  dura- 
tion, very  high  monitoring  light  output  is  desired.   This 
was  accomplished  by  using  a  conventional  signal  generator 
feeding  a  BNC  connector  in  the  rear  panel  of  the  M301  power 
supply.   The  timing  and  magnitude  of  the  positive  going 
sguare  pulse  was  determined  by  the  specifications  given  by 
the  manufacturer.121  The  pulse  width  is  reguired  to  be  <  1 
second  to  eliminate  damage  to  the  bulb.   The  pulse  magnitude 
can  be  varied  between  0  and  1  volt  according  to  the  follow- 
ing eguation 

lout  =  ^dle  +  tvmod/4)  x  5°1        IX-4 
where  Iout  =  M301  output  current 

Iidle  =  idle  current'  typically  3-5  amps 
V  d  =  output  voltage  of  the  signal  generator 
The  upper  limit  of  Vmod,  1  volt,  is  determined  so  that  Iout 
does  not  exceed  the  listed  nominal  operating  current  of  the 
lamp  which  for  the  Ushio  UXL  300  W  bulb  is  approximately  15 
amps.   To  attenuate  and  fine  control  the  5-volt  TTL  type 
signal  from  the  pulse  generator,  a  potentiometer  voltage 
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divider  circuit  was  placed  in  series  with  the  pulse  genera- 
tor and  M301  power  supply.   Preliminary  monitoring  of  the 
pulse  parameters  was  done  through  a  Tektronix  Model  35A 
oscilloscope  prior  to  actually  pulsing  the  lamp. 

2.3.2  He-Ne  Laser 

A  Spectra  Physics  Model  155  helium-neon  laser  was  used 
for  alignment  of  the  optical  systems,  calibration  of  the 
monochromator  gratings,  and  as  a  light  source  for  absorption 
spectrophotometry  involving  the  kinetics  of  the  hydrated 
electron,  ea  ~.   Characteristics  of  the  monochromatic  632.8 
laser  light  such  as  high  collimation,  low  divergence,  and 
high  intensity  make  this  source  ideally  suited  for  e  ~ 
studies.122   Figure  4  shows  the  optical  arrangement  used  for 
the  He-Ne  laser,  the  Jarrell-Ash  monochromator  and  side-on 
PMT.   Attempts  to  monitor  the  632.8  nm  light  with  the  end-on 
EMI  PMT  were  unsuccessful  due  to  this  tube's  lack  of  sensi- 
tivity above  600  nm. 

2.4   Reaction  Cells  and  Cellholder 

Two  types  of  reaction  cells,  both  consisting  of  modi- 
fied guartz  UV-VIS  spectrophotometric  cuvettes,  were  con- 
structed for  this  study.   Previous  accounts  of  similar 
reaction  cells  for  use  in  low  energy  electron  beam  radioly- 
sis  studies  involving  liguids  indicated  the  importance  of 
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using  very  thin  electron  entrance  windows  that  would  mini- 
mize the  attenuation  of  the  low  energy  electron  beam  in  the 
liguid.123  A  cuvette  fitted  with  a  reservoir  and  flow  exit 
for  flushing  purposes  was  used  in  all  the  on-line  absorption 
detection  work  and  a  non-flow  cell  was  used  exclusively  for 
the  Fricke  dosimetry  experiments.40   Both  cells  were  fitted 
into  an  aluminum  cell  holder  to  provide  reproducible  posi- 
tioning and  to  assure  the  delivery  of  equal  doses. 

2.4.1  Flow  Reaction  Cell 

The  flow  cell,  Figure  10,  was  constructed  from  a  60  ml 
separatory  funnel,  a  $  2  stopcock,  a  SB  5/20  glass  joint,  a 
glass/quartz  connector  seal  and  a  modified  quartz  UV-VIS 
cuvette.   A  thin  mylar  sheet  ~   200  microns  thick  was  held  by 
epoxy  glue  to  the  quartz  cell  and  served  as  the  electron 
entrance  window.  In  the  initial  construction  of  the  cell, 
FISHER  brand  microscopic  optical  glass  cover  slides  with  a 
thickness  of  0.13  to  0.17  mm  were  used.124   However,  the 
mylar  sheet  material  was  chosen  as  the  electron  entrance 
window  due  to  its  lesser  thickness  and  the  ease  by  which  it 
could  be  mounted  and  trimmed  on  the  cuvette.   The  optical 
pathlength  was  1  cm,  and  with  an  electron  penetration  depth 
of  ca.  1.5  mm  for  600  KeV  electrons  in  water,122  the  irradi- 
ated volume  was  1  x  1  x  0.15  cm.   This  electron  penetration 
depth,  which  is  less  than  the  extrapolated  range  of  1.9  mm, 
is  taken  as  the  depth  over  which  a  reasonably  uniform  dose 
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Fig.  10.   Flow-reaction  cell  with  modified  UV-VIS  cuvette. 
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deposition  can  be  assumed.125  For  studies  using  purging 
gases,  the  gas  was  continuously  bubbled  into  the  upper 
reservoir  for  the  duration  of  the  experiment. 


2.4.2  Non-Flow  Reaction  Cell 

This  cell,  Figure  11,  was  constructed  from  a  UV-VIS 
cuvette  to  provide  detection  of  304  nm  light  associated  with 
the  maximum  absorption  of  the  ferric  ion,  Fe3+,  the  oxida- 
tion product  derived  from  the  radiolysis  of  an  acid  solution 
of  ferrous  ion,  Fe2+,  in  the  Fricke  dosimetry  system  as 
discussed  in  Appendix  B.40  The  electron  entrance  window  was 
identical  to  that  of  the  flow  cell  and  a  total  volume  of  1.4 
ml  was  delivered  by  syringe  to  the  cell  prior  to  radiolysis. 
This  cell  could  be  readily  removed  from  the  aluminum  cell 
holder  in  the  pulse  radiolysis  system  and  transferred  to  a 
Beckman  DU  model  2400  spectrophotometer  with  Gilford  solid- 
state  electronics,  allowing  absorption  measurements  to  be 
made  as  a  function  of  the  number  of  radiation  pulses  re- 
ceived.  A  dilution  factor  involving  the  actual  irradiated 
volume  of  0.15  ml  and  the  total  volume  of  1.4  ml  was  used  in 
all  dose  calculations. 


2.4.3  Cell  Holder 

An  aluminum  cell  holder,  Figure  12,  was  used  to  enable 
reproducible  placement  of  the  solution  cell  with  respect  to 
both  the  electron  beam  and  analyzing  light.   The  holder 
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FRONT  VIEW 


1.4  cm 


1  cm 


mylar  sheet  e-beam 
entrance  window 


SIDE  VIEW 


ANALYZING  LIGHT 


solution  injection 
by  syringe 


i 


total  volume  =    1 .40  +/-  0.02  ml 


Fig.  11.  Non-flow-reaction  cell  formed  from  UV-VIS  cuvette. 
Reaction  volume  consists  of  1  x  1  x  0.15  cm  region 
immediately  behind  mylar  entrance  window. 
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6.0" 


wire  mesh  screens, 
60  -  80  % 
beam  reduction 


ANALYZING 
LIGHT 


adjustable  razor  blades, 
entrance  and  exit 
apertures 


Fig.  12.   Reaction  cellholder  with  adjustable  razor  blade 
mounts  and  wire-mesh  screens.   Electron  beam 
enters  through  rear  flange  and  traverses  8.0 
inches  in  air.   Radiation-induced  reaction  region 
occurs  in  cell  (not  shown)  located  at  intersection 
of  e-beara  and  analyzing  light. 
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consists  of  a  1  x  1  inch  open  rectangular  aluminum  column 
eight  inches  in  length  fitted  to  a  six-inch  diameter  alumi- 
num flange  that  is  attached  by  four  threaded  screws  to  the 
outer  edge  of  the  short-pulse  adapter  chamber  of  the  Febe- 
tron  706.   Adjustable  razor  blades  mounted  at  the  entrance 
and  exit  apertures  at  right  angles  to  the  electron  beam  axis 
defined  the  width  of  the  analyzing  light  beam.   The  use  of  a 
very  narrow  analyzing  light-beam  allowed  for  the  study  of 
the  solution  kinetics  by  monitoring  only  a  narrow  cross- 
section  for  which  the  dose  could  be  reasonably  defined. 

Using  this  arrangement,  the  electron  beam  traversed  ca. 
8  inches  of  air  between  the  field  emission  tube  exit  window 
and  the  thin  mylar  entrance  window  of  the  reaction  cell. 
This  distance  creates  a  slight  attenuation  of  the  electron 
beam  intensity  which  is  primarily  associated  with  a  spread- 
ing of  the  effective  beam  diameter  away  from  the  tube  face. 
For  a  distance  of  eight  inches  the  beam  will  have  an  effec- 
tive diameter  of  three  centimeters  as  compared  to  a  diameter 
of  1.5  centimeters  as  it  exits  the  tube  face.   The  loss  of 
beam  intensity  due  to  absorption  in  air  should  be  negligible 
due  to  the  maximum  penetration  range  of  80  inches  for  600 
KeV  electrons  in  air.125'126 

Experiments  with  FISHER  brand  optical  glass  microscopic 
coverslides  placed  in  the  cell  holder  at  the  solution  cell 
position  were  performed  to  determine  the  depth  distribution 
of  the  electron  beam.   Typical  depth-dose  curves  are  shown 


44 

in  Figure  13  which  includes  data  for  two  experimental  condi- 
tions.  Initially  a  stack  of  approximately  20  slides  was 
exposed  to  a  single  Febetron  shot  with  maximum  charging  of 
30KeV.   Absorption  measurements  were  then  made  on  each 
individual  slide  to  detect  the  amount  of  discoloration  in 
the  glass  due  to  the  e-beam  penetration.   Secondly,  a  simi- 
lar  stack  of  coverslides  was  exposed  to  10  cumulative  shots 
at  full  charge,  followed  by  individual  absorption  measure- 
ments.  Both  data  sets  indicate  an  extrapolated  range  for 
the  600  keV  electrons  in  the  glass  to  be  approximately  7 
slides.   Using  a  slide  thickness  of  0.15  mm,  the  extrapolat- 
ed depth  in  glass  is  found  to  be  1.0  mm.   Assuming  that 
glass  is  approximately  twice  as  dense  as  water,  these  data 
indicate  a  maximum  range  for  600  Kev  electrons  in  water  to 
be  approximately  2.0  mm,  which  agrees  with  the  earlier 
reported  value  of  1.9  mm.125 

Figure  14  shows  the  absorption  due  to  the  discoloration 
in  a  single  slide  as  a  function  of  the  number  of  pulses 
delivered  at  maximum  charging.   The  accumulated  discolora- 
tion resulting  from  the  e-beam  passage  through  the  glass 
slide  was  found  to  be  non-linear,  with  the  absorption  fall- 
ing off  as  a  function  of  pulse-number. 

Wire  mesh  screens  capable  of  60%  to  80%  electron  beam 
reductions  could  be  positioned  in  the  cell  holder  two  inches 
from  the  field  emission  tube  window  (Figure  12)  to  aid  in 
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Fig.  13.   Depth-dose  curves  for  glass  coverslide  experiment. 
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Fig.  14.   Absorbance  change  with  number  of  pulses  delivered 
for  glass  coverslide  experiment. 
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acquiring  transient  formation  signals  occurring  during  the 
first  few  microseconds  immediately  following  the  radiolysis 
pulse.   These  screens  restrict  the  flux  of  the  electron 
beam  which  was  maintained  at  a  maximum  charging  voltage  of 
30  KeV.   With  the  wire  mesh  screens  inserted,  attenuation  of 
the  beam  intensity  decreases  the  total  average  population  of 
initial  radicals  formed.   This  beam  flux  attenuation  helps 
to  prevent  radical-radical  interactions  and  provide  favora- 
ble radical-solute  interaction  conditions.   However,  due  to 
both  a  malfunction  in  the  fastest  data  acquisition  mode  of 
the  transient  digitizer,  and  electromagnetic  interferences 
associated  with  the  Febetron  706  system  that  were  discussed 
earlier,  attempts  to  capture  transient  formation  signals 
with  enough  resolution  for  kinetic  analysis  were  not  suc- 
cessful. 

2 . 5  Monochromators 

In  the  early  stages  of  this  work  a  0.3  meter  McPherson 
218  monochromator  was  used.   Later,  efforts  were  made  to 
reduce  electronic  interferences  associated  with  the  electron 
beam  pulse  by  removing  the  monochromator  and  PMT  couple  away 
from  and  to  the  rear  of  the  field  emission  tube.   Due  to  the 
rather  permanent  fixture  of  the  Mcpherson  218  and  its  asso- 
ciated vacuum  pump  assembly  relative  to  the  initial  experi- 
mental geometry,  a  second  0.25  meter  Jarrell-Ash  82-41550 
monochromator  was  employed. 
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2.5.1  0.3  Meter  Monochromator 

Initial  studies  were  carried  out  using  a  McPherson  218 
vacuum/atmospheric  monochromator  which  has  an  f/5.3,  a  2400 
G/mm  plane  grating  and  a  1.33  nm/mm  linear  dispersion.   The 
vacuum  capabilities  were  not  employed  for  any  of  this  work 
and  the  2400  G/mm  grating  having  a  0  -  5,000  angstrom  wave- 
length range  was  eventually  replaced  with  a  1200  G/mm  grat- 
ing capable  of  operating  in  the  wavelength  range  from  0  - 
10,000  angstroms.   A  synchronous  motor-driven  gear  box 
supplies  12  scanning  speeds  from  0.5  angstroms /minute  to  200 
angstroms/minute  and  slit  widths  were  adjustable  from  10 
microns  to  2  millimeters.   Slit  widths  of  1.5  millimeters 
and  a  fixed  slit  height  of  8  millimeters  were  used.   Figure 
3  shows  the  optical  configuration  used  with  the  McPherson 
monochromator.   Light  exiting  the  xenon  arc  lamp  was  passed 
through  a  condensing  lens  onto  a  first  surface  mirror  and 
reflected  through  the  solution  cell.   The  beam  was  then 
focused  on  the  entrance  slit  of  the  monochromator  and  after 
wavelength  dispersion  the  exiting  beam  was  reflected  onto 
the  photocathode  of  the  PMT  by  a  2.5  cm  focal  length  para- 
bolic mirror. 

2.5.2  0.25  Meter  Monochromator 

Initial  testing  of  the  Jarrell-Ash  monochromator  cou- 
pled to  the  end-on  PMT  with  both  components  located  approxi- 


49 


mately  2  meters  away  from  the  field  emission  tube  showed  a 
marked  improvement  in  reducing  the  electromagnetic  noise 
interferences.   Using  the  McPherson  monochromator  in  the 
previously  shown  geometry  in  conjunction  with  the  gating 
circuit  produced  a  net  dead-time,  caused  by  gross  overload- 
ing of  the  PMT,  of  >  5  microseconds  in  the  detection  system. 
When  the  detection  system  was  relocated,  as  shown  in  Figure 
4 ,  net  dead-times  of  2  to  3  microseconds  could  be  attained 
without  the  use  of  the  gating  circuit.   All  subseguent  work 
was  completed  with  this  arrangement  and  eventually  the  end- 
on  PMT  was  replaced  by  the  side-on  PMT.  To  attain  further 
reduction  of  the  EMI  effects,  the  entire  detection  system 
and  signal  lines  would  have  to  be  either  completely  removed 
from  the  room  containing  the  Febetron  706  or  extensively 
shielded.120'127 

In  this  modified  configuration  the  analyzing  light, 
after  passing  through  the  solution,  was  reflected  by  a  first 
surface  mirror  back  to  a  plano-convex  f/2.5  condensing  lens 
and  focused  onto  the  entrance  slit  of  the  Jarrell-Ash  mono- 
chromator.  This  monochromator  has  an  f/3.5,  a  1180  G/mm 
grating  and  a  3.3  nm/mm  linear  dispersion.   A  wavelength 
motor  drive  of  100  nm/minute  was  coupled  to  a  wavelength 
drive  knob  to  enable  scanning  of  the  light  source  spectrum. 

The  slit  arrangement  provided  with  this  monochromator 
was  altered  to  accommodate  the  direct  coupling  of  the  Jar- 
rell-Ash monochromator  to  both  a  previous  gas  phase  aluminum 
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reaction  cell,  and  the  Pacific  Precision  Instruments  PMT 
housing.128   In  this  work  only,  the  direct  couple  of  the 
monochromator  and  PMT  housing  was  utilized.   This  was  accom- 
plished by  screwing  the  PMT  housing  containing  the  side-on 
Hamamatsu  tube  directly  to  the  exit  slit  flange  adapter. 
Figure  15  shows  the  entrance  and  exit  system  flange  adapt- 
ers.  Brass  slit  bodies  of  fixed  width  and  height,  0.71  mm 
and  15  mm  respectively,  were  screwed  into  the  aluminum 
flange  by  turning  the  two  small  holes  in  the  slit  body  with 
a  sharp  tool . 

2.6   Febetron  706  System 

The  principle  of  pulse  radiolysis  involves  the  delivery 
of  ionizing  radiation  to  a  chemical  system  in  the  form  of  a 
short  pulse  that  creates  a  non-eguilibrium  system  with  an 
adeguate  amount  of  transient  species  that  are  monitored  in 
time.   An  accelerator  is  used  as  the  source  of  the  high 
kinetic  energy  ions  and  several  related  parameters  are 
important  in  pulse  radiolysis  measurements.   Primarily,  the 
radiolytic  pulse  should  be  very  short  relative  to  relaxation 
of  the  transient  system,  and  the  pulse  current  should  be 
large  enough  to  produce  an  observable  concentration  of 
transients.   Energy  deposition,  which  depends  on  1)  the 
particle's  mass,  charge,  and  velocity,  and  2)  the  medium 
density,  is  critical  in  determining  the  distribution  (either 
homogeneous  or  non-homogeneous)  of  radical  concentrations 
produced . 


slit  area  =  0.71  mm  x  15  mm 
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Fig.  15.   Entrance  and  exit  slit  system  flange  adaption  to 
the  Jarrell-Ash  monochromator . 
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The  Febetron  706  electron  accelerator  used  in  this  work 
is  fitted  with  a  Model  5515  high  vacuum  field  emission 
electron  tube  and  has  the  following  specifications:  12 
joules  maximum  pulse  energy,  8,000  amperes  of  peak  current, 
3  nanosecond  pulse  duration  (FWHM),  5%  peak  to  peak  repeat- 
ability, and  a  maximum  electron  energy  of  600  KeV.   The 
extremely  high  currents  make  this  model  particularly  ap- 
plicable to  gas  phase  studies. 

However,  it  was  suspected  that  the  high  vacuum  electron 
tube  had  become  gaseous  due  to  intermittent  use  in  previous 
years  and  the  extended  lifetime  of  the  electron  tube  beyond 
its  suggested  3500  pulses.  The  e-beam  tube  had  failed  after 
ca.  8000  pulses;  the  beam  had  burned  microscopic  holes  in 
the  window  which  allowed  the  pressure  in  the  tube  to  rise  to 
atmospheric.   It  was  expected  that  partial  operation  could 
be  restored  by  replacing  the  window  and  reevacuating  the 
tube.   This  attenuated  beam  feature  would  be  beneficial  to 
condensed  phase  studies  by  eliminating  radical-radical 
interactions  caused  by  the  absorbance  of  enormous  radiation 
doses.   Encouraged  by  other  somewhat  analogous  and  confirm- 
ing reports,122'129-132  we  pursued  the  application  of  the 
Febetron  706  accelerator  to  agueous  solution  work. 

Details  of  the  Febetron  706  system  can  be  found  else- 
where.125'126'133  Figure  16  shows  a  block  diagram  of  the 
Model  706/2677  system.   The  Febetron  706  system  was  main- 
tained by  keeping  the  various  chambers  with  their  associated 


53 


o 

+> 

(0 

! 


(0 

>i 

(0 
O 

C 

o 

M 

0) 


H 

• 

CP 

■H 

Cm 


54 


pressure  regulators  and  gauges  properly  pressurized.   A 
compressed  air  tank  is  used  to  pressurize  the  module  chamber 
at  a  maximum  of  80  psi  and  the  short  pulse  adapter  at  a 
maximum  of  240  psi.   Freon-12  was  maintained  at  13  psi  in 
the  pulser  case.   The  transformer  oil  in  the  short  pulse 
adapter  was  circulated  daily  by  an  electric  pump  connected 
to  a  large  oil  reservoir. 

The  accelerator  operates  on  the  principle  of  the  Marx- 
bank  circuit  which  consists  of  a  number  of  connected  high 
voltage  capacitors  that  are  charged  in  parallel  and  subse- 
guently  rapidly  switched  to  a  series  arrangement  through 
spark  gaps.   This  circuit,  located  in  the  Marx-surge  stage 
of  the  pulser,  produces  a  high  voltage,  extremely  high  peak 
power  pulse  that  is  used  to  pulse  charge  through  a  series 
spark  gap  in  the  central  electrode  of  the  short  pulse  adapt- 
er stage  of  the  pulser.   This  short  pulse  adapter  then 
discharges  an  extremely  short  3  ns  high  current  pulse  to  the 
field  emission  electron  tube.   This  pulse  generates  electron 
emission  from  the  cathode  and  the  electrons  are  then  accel- 
erated toward  an  exit  window  that  functions  as  the  anode. 

As  mentioned  above,  due  to  the  lack  of  an  inward  de- 
flection in  the  thin  exit  window  of  the  field  emission  tube, 
it  was  decided  that  gaseous  conditions  existed  inside  the 
tube.   Under  normal  conditions,  the  thin  metal  window  would 
be  noticeably  drawn  in  by  the  high  vacuum.   Experiments 
using  the  technigue  of  calorimetric  dosimetry  carried  out 


55 


directly  in  front  of  the  tube  exit  window  showed  negligible 
energy  outputs.   Therefore  a  modified  electron  beam  tube 
exit  system  was  fabricated. 

2.6.1  Refurbished  Beam  Tube 

A  modified  electron  beam  tube  flange  and  exit  window 
system132  was  constructed  to  renew  the  output  of  the  Febe- 
tron  706  system.   This  modification,  Figure  17,  was  formed 
from  a  stainless  steel  chamber  4.75  inches  in  diameter  and 
1.25  inches  thick,  and  fit  with  a  Model  202-100  Perkin-Elmer 
2  liter/second  differential  ion  pump  powered  by  a  Model  222- 
0200  Perkin-Elmer  IONPAK  200  power  supply.   The  high  vacuum 
current  limiting  power  supply  operated  at  a  voltage  between 
+4000  volts  DC  and  +6000  volts  DC  and  the  system  was  capable 
of  attaining  pressures  below  1  x  10~6  Torr. 

The  former  anodic  window  of  the  field  emission  tube  was 
cut  away  and  a  thin,  approximately  0.001  inch  stainless 
steel  window  and  o-ring  were  set  into  the  compartment.   The 
compartment  was  attached  by  three  threaded  screws  to  the 
field  emission  tube  window  face  plate  flange  of  the  Febetron 
706  with  a  second  o-ring  seal.   Initially,  the  thin  window 
was  placed  0 . 8  inches  away  from  the  original  anode  window 
location.   Although  an  adeguate  vacuum  was  attained,  further 
dosimetry  experiments  showed  no  improvement  with  this  geome- 
try so  the  window  was  set  further  into  the  compartment  0 . 4 
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inches  away  from  the  original  window  location.   After  rough 
pumping  the  chamber  down  to  a  few  microns  of  pressure  and 
attaining  the  minimum  +4000  volts  DC  reguired  for  proper 
operation  of  the  ion  pump,  calorimetric  dosimetry  measure- 
ments were  again  performed. 

2.6.2  Calorimetric  Dosimetry 

The  total-beam  calorimetric  technigue,  as  illustrated 
in  Figure  18,  involves  the  measurement  of  the  heat  deposi- 
tion in  an  aluminum  disk  upon  firing  the  Febetron  706  at 
maximum  charge.   The  thickness  and  diameter  of  the  disk, 
0.033  inches  and  1.385  inches  respectively,  were  chosen  to 
assure  complete  e-beam  deposition.   Aluminum  is  suitable  due 
to  its  high  thermal  conductivity,  specific  heat,  and  low 
backscatter  energy  coefficient.   The  temperature  of  the 
calorimeter  was  conveniently  measured  by  attaching  a  thermo- 
couple and  measuring  the  voltage  output  with  a  recording 
millivoltmeter.   The  thermocouple  was  attached  to  the  alumi- 
num disc  in  this  particular  calorimeter  by  drilling  a  small 
hole  in  the  center,  inserting  the  thermocouple  wires,  and 
swaging  the  aluminum  around  them  with  a  punch.   The  calibra- 
tion of  the  thermocouple/chart  recorder  was  accomplished 
using  room  temperature  of  ~   25  °C  as  the  base  line  and 
normal  body  temperature  of  ~   37  °C  as  the  maximum  scale 
value.   For  the  maximum  temperature  reading,  the  disk  was 
held  firmly  in  the  palm  of  the  hand  until  a  steady  eguili- 
brated  signal  was  obtained. 
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Fig.  18.   Total  beam  calorimetric  dosimetry  arrangement. 
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Using  the  equation 

Q  =  m  C_  A  T  II-5 

P 

with  Q  =  heat  required  to  raise  the  specific  heat  of  the 
calorimeter  by  A  J,  Joules 
m  =  2.3  +/-  0.05  grams  (mass  of  the  aluminum  disk) 
C_  =  0.22  cal/  °C  gram  (specific  heat  of  aluminum) 
AT  =  0.53  °C  (difference  between  pre-  and  post-  pulse 

read  from  chart  recorder) 
a  value  of  1.12  +/-  0.06  Joules  was  determined  as  the  beam 
energy  absorbed  per  pulse.   These  positive  results  indicated 
that  the  placement  of  the  new  anode  window  was  a  critical 
factor  with  possible  arcing  complications  arising  when  the 
new  window  was  not  located  close  enough  in  to  the  original 
window  position. 

2.7   Transient  Recorder  and  Computer  System 

While  early  pulse  radiolysis  studies  used  the  oscillo- 
scope-Polaroid photograph  technique  for  recording  transient 
signals,  development  of  sophisticated  transient  recorders 
followed  that  now  have  replaced  this  manual  data  analysis 
method.117   Digital  recorders  now  digitize  the  input  signal 
in  real  time  and  can  resolve  the  signal  into  separate  chan- 
nels thus  providing  a  stored  digitized  record  of  the  tran- 
sient trace  as  a  function  of  time.  In  addition,  interfacing 
of  the  recorder  with  a  computer  allows  for  integrated  con- 
trol of  the  operational  parameters  such  as  time  base,  sensi- 
tivity, delays,  signal  processing,  and  data  transfer. 


60 


The  data  acquisition  and  analysis  system  developed  for 
this  work  is  a  second  generation  system  having  been  preceed- 
ed  by  a  Biomation  Model  610B  transient  recorder  supported  by 
an  IMSAI  8080  S-100  bus  microcomputer.128   Current  instru- 
ments include  an  Analogic  Data  Precision  DATA  6000  waveform 
analyzer  fitted  with  a  Model  620  2-channel  plug-in  digitizer 
interfaced  through  a  serial  port  at  9600  baud  to  a  Tandy 
1200  IBM-PC/compatible  laboratory  microcomputer. 

The  DATA  6000  has  two  signal-input  ranges  of  +/-  1.2 
and  +/-  3.6  volts  full  scale,  and  a  maximum  sampling  rate  of 
100  MHz.   This  sampling  rate  gives  a  period  (time  between 
sampled  points)  of  10  nanoseconds  for  collection  of  512 
points  using  an  8-bit  recorder.   Although  the  maximum  period 
attainable  is  300  seconds,  a  maximum  period  of  1  microsecond 
was  used  in  this  work  giving  a  sampling  rate  of  1  MHz.  The 
DATA  6000  is  noted  for  its  versatile  capabilities  of  signal 
capture,  processing  and  data  transfer,  although  the  lack  of 
any  fine  control  on  the  input  voltage  parameter  was  found  to 
be  a  rather  severe  experimental  limitation  in  this  work. 
The  DATA  6000  is  controlled  by  a  single  internal  68000 
microprocessor  and  a  CRT  display  screen  enables  visual 
inspection  of  signals  before  data  transfer.134 

Data  recording  parameters  were  selected  either  manually 
through  keypads  on  the  front  panel  or  remotely  via  computer 
control.  Remote  control  of  the  DATA  6000  via  an  RS-23  2 
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serial  interface  port  utilized  three  basic  functions:  con- 
trol of  data  recording  parameters,  use  of  processing  and 
display,  and  transfer  of  data  to  and  from  the  digitizer.   To 
use  these  functions,  a  computer  program  written  in  BASIC 
language  provided  either  keystroke  emulation  by  key  code 
number  or  direct  control  using  English  commands.   These 
operations  were  sent  to  the  DATA  6000  as  string  commands  and 
data  transfers  to  the  computer  were  sent  in  ASCII  number 
format.   Various  timing  loops  set  up  in  the  BASIC  program 
assured  proper  hand-shaking.   During  operation  several 
parameters  such  as  trigger  source,  timebase,  screen  display 
and  marking  were  initially  set  manually  and  subseguently 
controlled  remotely  through  the  computer  program.   For 
short-lived  transient  signals,  the  multiple  trace  delayed 
mode  with  signal  averaging  was  used  which  provides  an  en- 
hancement to  the  S/N  ratio  in  proportion  to  the  sguare  root 
of  the  number,  n,  of  signals  averaged.118   Otherwise,  the 
single  trace  delay  mode  was  used  with  no  signal  averaging. 
As  shown  in  Figure  2,  the  triggering  and  delays  were 
originated  from  a  single  initial  manual  trigger  pendant 
delivering  a  90  volt  pulse  derived  from  three  30  volt  bat- 
teries in  series  with  a  standard  RC  trigger  circuit.   A 
typical  timing  seguence  is  presented  in  Figure  19.   Delay  of 
the  Febetron  firing  was  accomplished  using  an  Evans  Associ- 
ates programmable  time  delay  (PTD)  module,  Model  4145-2. 
Three  such  modules  are  present  which  could  be  used  in  series 
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for  a  maximum  delay  of  3  milliseconds.   The  voltage  divider 
and  amplifier  in  the  circuit  are  required  because  the  incom- 
ing trigger  for  the  Febetron  must  be  90  volts  in  magnitude 
while  the  PTD  uses  a  TTL  type  input  and  output.  The  lamp 
pulsing  and  shutter  control  was  provided  by  a  Vincent  Asso- 
ciates shutter  driver  and  timing  unit,  Model  SD  10,  which 
controlled  a  Model  225L  Vincent  Associates  UNIBLITZ  shutter. 
The  shutter  was  not  incorporated  into  the  timing  sequence 
but  was  manually  switched  on  and  off  during  a  calibration 
procedure  to  set  the  zero  and  full  light  levels,  VQ  and  Vm 
respectively.   The  lamp  pulse  duration  was  controlled  by  the 
shutter  driver  and  timing  unit  which  provides  a  positive 
going  5  volt  square  pulse  with  a  variable  8  to  8,000  milli- 
second duration. 

2  .  8  Computer  Program 

Full  interaction  between  the  operator  and  the 
computer/transient  recorder  was  provided  by  a  single  BASIC 
program.   This  program,  developed  by  J.  Brogdon  and  R.  J. 
Hanrahan  with  assistance  from  B.  Yarborough  and  Brent  Jones, 
was  locally  written  in  accordance  with  suggested  BASIC  rou- 
tines provided  by  the  DATA  6000  manufacturer.135   A  complete 
listing  is  provided  in  Appendix  C.   The  BASIC  source  code 
was  compiled  from  the  RUN  menu  of  Microsoft  Quick  BASIC 
version  4.0,  and  was  executed  directly  from  memory  or  from  a 
separate  executable  program  written  to  disk. 
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Initially,  the  opening  menu  would  offer  all  options 
available,  and  upon  full  completion  of  any  commands  the  main 
menu  would  again  be  presented.   Prior  to  data  collection  the 
CALIBRATION  command  was  used  to  set  all  data  recording 
parameters  and  to  provide  a  dark  or  no  light  trace,  VO ,  and 
a  full  light  trace,  VM.   The  DATA  command  was  entered  during 
collection,  storage,  and  averaging  of  actual  traces  by  the 
DATA  6000.   Upon  completion  of  this  loop  three  files  were 
stored  into  computer  memory:  all  digitized  voltage  values 
(QY(I)),  a  time  dependent  concentration  file  [QX(I), 
CON(I)],  and  a  time  dependent  optical  density  file  [QX(I), 
0D(I)].   After  transfer  of  the  data  from  the  digitizer  to 
the  computer,  the  SAVE  command  could  be  executed  to  save  a 
disk  copy  of  the  digitized  voltage  values.   Using  the  PLOT- 
TER command  provided  a  graphical  presentation  of  the  digi- 
tized voltage  optical  trace  on  the  computer  CRT  screen. 
Often  hard  copy  of  the  graphical  displays  would  be  obtained 
to  catalog  the  data. 

The  subroutine  at  line  1570  (LEFT  &  RIGHT  BORDERS) 
allowed  for  selection  of  any  portion  of  the  trace  followed 
by  storage  of  either  a  concentration  or  optical  density  file 
corresponding  to  the  section  of  the  trace  within  the  bor- 
ders.  The  subroutine  at  line  3980  (LINEAR  REGRESSIONS) 
provided  linear  least  sguares  analysis  of  the  concentration 
data  to  aid  in  determining  the  kinetics  involved  in  the 
transient  decays.   Operation  of  the  MANUAL  command  enables 
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adjustment  for  the  true  full  light  signal,  VM.  Even  though 
the  100%  light  signal  was  set  during  the  calibration  steps, 
any  fluctuation  in  the  light  during  the  actual  data  collec- 
tion mode  could  be  accounted  for  by  taking  an  average  of  the 
light  signal  immediately  preceding  the  pulse  as  an  accurate 
measure  of  the  true  100%  light  intensity.  Selection  of  the 
molar  absorption  coefficient,  E,  which  was  dependent  on  the 
specific  transient  being  observed  was  also  done  through  the 
MANUAL  command.   The  DIRECTORY  command  provided  directory 
inspection  to  insure  the  storage  of  all  files  to  disk  and  to 
access  available  memory. 


2.9  Reagents  and  Sample  Preparation 
2.9.1  Sodium  Tetraphenylborate  and  Sodium  Azide 


Kodak  reagent-grade  sodium  tetraphenylborate  and  sodium 
azide  were  used  as  received. 

2.9.2  IrCl63"( hydrate)  and  IrCl62" 


The  sodium  hexachloroiridate(III)  hydrate  and  sodium 
hexachloridate(IV)  complexes,  which  were  obtained  from  Al- 
drich  Chemical  Company,  were  stored  in  a  dessicator  and  used 
without  purification. 

2.9.3  Sodium  Carbonate 

Fisher  Scientific  Company-supplied  Na2C03  was  used  in 
all  the  carbonate  competition  experiments. 
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2.9.4  Sodium  Hydroxide  and  Ethanol 

Fisher  Scientific  Company-supplied  NaOH  and  USP  reagent 
quality  200  proof  absolute  CH3CH2OH  from  Florida  Distillers 
Company  were  used  in  the  hydrated  electron  dosimetry  stud- 
ies. 

2.9.5  Sulfuric  Acid  and  Ferrous  Ammonium  Sulfate 

Fisher  Scientific  Company  H2S04  and  Fe(NH4 )2 (S04 )2' 6H20 
were  used  as  received  for  the  Fricke  dosimetry  studies. 

2.9.6  Potassium  Thiocyanate 

Fisher  Scientific  Company  KSCN  was  used  as  received  for 
the  thiocyanate  dosimetry  studies. 


2.9.7  Tert-Butyl  Alcohol 


Fisher  Scientific  Company  reagent-grade  (CH3)3COH  was 
used  as  received  for  the  hydrated  electron  rate-constant 
studies. 

2.9.8  Potassium  Permanganate 

Fisher  Scientific  Company  reagent-grade  KMN04  and  NaOH 
were  used  as  received  for  the  water  distillation  procedure. 
KMN04  solid  was  dissolved  in  laboratory-distilled  water  that 
was  made  alkaline  by  addition  of  NaOH.   Upon  distillation, 
any  organic  impurities  were  removed  by  an  oxidation  reaction 
involving  the  Mn04~  ion. 
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2.9.9  Oxygen  and  Nitrogen 

Liquid  Air  Corporation-supplied  02  and  N2,  both  of 
which  were  >  98%  purity,  were  used  as  received  for  purging 
gases  in  the  'Super'  Fricke  studies  and  the  hydrated  elec- 
tron studies,  respectively. 

2.9.10  Nitrous  Oxide 

Matheson  CP-grade  N20  with  >  99.0%  purity  was  used  as 
received  for  generating  oxidizing  conditions  for  all  studies 
involving  the  OH*  radical  reactions. 

2.9.11  Sample  Preparation 

Aqueous  solutions  were  prepared  by  dissolving  known 
amounts  of  solid  reagents  in  doubly  distilled  water.   A 
Model  B5  Mettler  analytical  balance  with  200  gram  capacity 
and  +/-  0.05  milligram  accuracy  was  used  in  all  mass  meas- 
urements. Reference  spectra  of  the  Iridium(III)  chloride 
complex,  Iridium(IV)  chloride  complex,  and  the  Sodium  Tetra- 
phenylborate  solution  were  measured  using  a  Hewlett-Packard 
model  2400  UV/VIS  spectrophotometer. 

2.10   Sample  Irradiation 

Prior  to  radiolysis,  solutions  which  were  to  be  irradi- 
ated under  oxidizing  conditions  were  purged  for  approximate- 
ly thirty  minutes  with  nitrous  oxide,  N20.   For  work  in- 
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volving  reducing  conditions,  nitrogen,  N2 ,  was  used  as  the 
purge  gas.   Purging  of  the  solutions  was  done  by  introducing 
a  disposable  pipette  directly  into  the  60  milliliter  separa- 
tory  funnel  serving  as  the  upper  reservoir  of  the  solution 
flow  cell.   Gas  flow  was  maintained  throughout  the  pulse 
experiments.   Prior  to  calibrating  the  DATA  6000,  the  solu- 
tion was  introduced  by  gravity  flow  into  the  modified  cu- 
vette, which  was  then  inserted  into  the  aluminum  cell  hold- 
er.  The  cell  holder  was  designed  to  accommodate  the  solu- 
tion cell  to  give  reproducible  positioning  and  uniform 
irradiation  in  the  region  that  was  sampled  by  the  analyzing 
light  beam.   After  completing  the  calibration  procedure,  the 
data  collection  routine  was  executed  and  the  Febetron  was 
charged  and  triggered.   Typically,  due  to  product  buildup  or 
gas  bubble  evolution  on  the  inner  surface  of  the  electron 
beam  entrance  window,  a  fresh  aliguot  was  flushed  into  the 
cell  and  the  system  was  recalibrated  after  each  pulse. 


CHAPTER  III 

THE  PULSE  RADIOLYSIS  OF 

SODIUM  TETRAPHENYLBORATE  IN  AQUEOUS  SOLUTION 

3.1   Experimental  Results 

3.1.1  Oxidation  Intermediates  of  Tetraphenylborate  Ions 

The  pulse  radiolysis  of  aqueous  sodium  tetraphenylbo- 
rate (TBP)  solutions  under  oxidizing  conditions  was  studied 
with  two  different  methods  of  generating  the  oxidizing 
species.   As  mentioned  previously,  using  N20-saturated 
solutions  the  predominant  oxidant  produced  is  the  OH*  radi- 
cal generated  from  the  direct  radiolysis  of  water  and  the 
interconversion  of  eaa~  to  OH-  by  N20.   Alternatively,  a 
secondary  oxidant,  N3«,  is  prepared  under  conditions  of  N20- 
saturation  and  the  addition  of  excess  sodium  azide,  NaN3. 
The  azidyl  radical,  N3 • ,  is  formed  from  the  oxidation  of  the 
azide  anion  by  the  0H«  radical, 

OH*  +  N3"   — >   0H~  +  N3«         III-l 

Both  the  N3~  anion  and  N3«  radical  have  absorption  maxima  at 
wavelengths  <  280  nm.  36   Experiments  with  N3«  radical  were 
carried  out  to  utilize  the  selective  electron  transfer  path- 
way of  oxidation  as  compared  to  other  possible  reactions  of 
the  OH*  radical  with  aromatic  systems,  namely  abstraction  of 
an  H  atom  and  addition  of  the  OH*  radical  to  the  aromatic 
ring. 
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The  pulse  radiolysis  of  a  nitrous  oxide-saturated  (N20 
a  2.4  x  10-2  M)  solution  of  10  mM  NaTPB  was  carried  out  at 
room  temperature.   Using  the  Xenon  arc  lamp  in  the  pulsing 
mode,  the  transient  absorption  signals  were  recorded  at 
various  wavelengths  from  300  to  800  nm.   Figure  20  shows  the 
absorption  signal  recorded  at  3  30  nm  for  a  single  pulse 
experiment.   The  appearance  of  an  incipient  secondary  ab- 
sorption occurring  approximately  20  to  30  microseconds  after 
the  initiating  e-beam  pulse  was  characteristic  in  all  the 
time-concentration  traces  recorded  in  the  300-400  nm  range. 

A  transient  absorption  spectrum,  Figure  21,  was  gener- 
ated by  plotting  the  maximum  absorption  or  "end  of  pulse" 
values  as  a  function  of  the  wavelength  at  which  they  were 
recorded.   This  spectrum  shows  an  absorption  maximum  at  330 
nm;  the  only  other  transient  spectral  feature  is  a  small 
absorption  (ca.  10%  of  the  main  peak)  centered  at  about  750 
nm.   This  spectrum  is  to  be  compared  with  that  reported  by 
Liu  et  al.  for  the  radical  derived  by  OH*  addition  to  the 
phenyl  groups  of  TPB-  —  a  spectrum  with  a  large  absorption 
at  300  nm  that  was  reported  to  change  only  slightly  over  the 
first  few  microseconds.76 

Assuming  the  absorption  at  330  nm  is  caused  by  the 
radical  derived  by  the  addition  of  OH*  to  the  phenyl  groups 
of  TPB-,  its  extinction  coefficient  at  the  absorption  maxi- 
mum of  3  30  nm  can  be  calculated  from  the  transient  spectra 
shown  in  Figure  21  (a)  by  making  use  of  the  value  for  G(0H) 
in  N20-saturated  solutions  egual  to  5.4  (=  G(OH)  +  G(e   ~). 
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Fig.  20.   Graphical  displays  of  optical  signal  versus  time 
following  pulse  radiolysis  of  N2G— saturated  10  mM 
NaTPB  solution.   a)  recorded  at  330  nm,  200  nsec 
period,  100  usee  full  scale;   b)  recorded  at  345 
nm,  500  nsec  period,  250  usee  full  scale. 
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In  N20-saturated  TPB"  solutions,  this  total  OH*  yield  of  5.4 
will  have  been  completely  converted  to  the  TPBOH-*  radical 
by  the  time  optical  measurements  can  be  made.   Thus  the  end 
of  pulse  optical  density  can  be  used  to  measure  the  extinc- 
tion coefficient  of  the  TPBOH-*  adduct.   The  value  obtained 
from  Figure  21  is  approximately  (10,000  +/-  10%)  M-1cm-1  for 
the  maximum  peak  at  330  nm. 

The  optical  densities  were  calculated  from  the  initial 
pre-pulse  light  intensity  determined  during  calibration,  Vf , 
and  from  the  transmitted  optical  signal,  V^.,  using  the 
eguation 

OD  =  -log  (Vt  /  Vf)  -  -log[  (QY(I)  -  VQ)  /  Vf)  ]      III-2 

Vf  is  obtained  from  subtracting  the  dark  signal,  VQ,  from 
the  initial  100%  light  intensity,  Vm,  determined  under  the 
experimental  conditions  prior  to  irradiation  and  V+    is  ob- 
tained from  subtracting  VQ  from  the  digitized  optical  sig- 
nal, QY(I). 

The  pulse  radiolysis  of  a  N20-saturated  solution  of  10 
mM  NaTPB  ,  with  excess  NaN3  »  0.2  M,  was  carried  out  at 
various  wavelengths  yielding  typical  absorption  traces,  as 
shown  in  Figure  22.   The  accompanying  transient  absorption 
spectrum  is  plotted  in  Figure  21.   The  transient  absorption 
traces  resulting  from  the  NaTPB  with  added  NaN3  conditions 
also  illustrate  a  secondary  transient  production.   The 
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Fig.  22.   Graphical  display  of  optical  signal  versus  time 

following  pulse  radiolysis  of  N20- saturated  10  xnM 
NaTPB,  0.2  M  NaN3  solution.   a)  recorded  at  330 
nm,  200  nsec  period,  100  usee  full  scale;   b) 
recorded  at  330  nm,  500  nsec  period,  250  usee  full 
scale. 
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transient  absorption  spectrum  contains  two  absorption  bands, 
one  at  330  nm  and  a  broad  band  above  700  nm.   This  transient 
absorption  spectrum  as  well  as  the  associated  kinetics  of 
the  330  nm  decay  are  in  good  agreement  with  published  re- 
sults obtained  under  similar  solution  conditions.  5'76 

3.1.2  Transient  Kinetics 

Kinetic  analysis  of  the  optical  traces  was  carried  out 
using  a  curve  fitting  routine  from  the  BASIC  program,  and 
the  best  fit  was  judged  from  the  correlation  coefficient  and 
visual  inspection.   Figures  23  and  24  show  typical  curve 
fits  for  both  the  initial  and  secondary  transients  produced 
under  N20-saturated  millimolar  NaTPB  conditions.   Analogous 
curve  fitting  was  carried  out  on  the  experimental  absorption 
traces  obtained  from  N20-saturated  solutions  of  NaTPB  con- 
taining excess  NaN3.   The  resulting  decay  rate  constants  are 
presented  in  Table  2. 

3.1.3  Carbonate  Competition  Method 

The  reaction  of  0H«  radicals  with  the  TPB~  anion  occurs 
too  rapidly  to  be  observed  directly  using  the  experimental 
conditions  previously  described.   However,  the  formation 
rate  of  the  OH«  radical  reaction  with  TPB~  can  be  obtained 
by  the  carbonate  competition  method137  using  a  secondary 
scavenger,  C032~,  which  competes  with  the  TPB~  anion  for  the 
0H«  radical.   The  two  competing  reactions  are 
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Fig.    23.      Curve   fit   of   data;    a)    transient   trace   presented 

from  Figure    20    (a);    b)    First-order   fitting  routine 
applied  to  data   between  x  =   20   and   x   =   40   micro- 
seconds. 
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Fig.  24.   Curve  fit  of  data;  a)  transient  trace  presented 
from  Figure  20  (b);  b)  Second-order  fitting  rou- 
tine applied  to  data  between  x  =  80  and  x  =  250 
microseconds . 


78 


Table  2 
Transient  Decay  Rate  Constants 


N20-saturated    Initial  Transient 
Solutions       Decay  Rate  Constant 

(  x  104  sec-1  ) 


Secondary  Transient 
Decay  Rate  Constant 
(  x  108  M-1sec-1  ) 


10  mM  NaTPB 


5.88  ±  0.98 


3.76  ±  0.42 


10  mM  NaTPB 
0.2  M  NaN, 


8.25  ±  0.59 

a)  2 

b)  3.9  ±  0.2 


a)  reported  by  Horii  and  Taniguchi 

b)  reported  by  Liu  et  al . 


3.10  ±  0.91 

a)  4 

b)  — 
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OH*  +  TPB~    >    products  ,  kx  III-3 

OH«  +  C032-   >    OH"  +  C03~  (detected)  ,  k2    130 

Using  the  He-Ne  laser  experimental  setup,  the  C03~  radical 
anion  was  monitored  as  a  function  of  added  TPB-  anion.  The 
C03~  radical  anion  has  an  extinction  coefficient  on  the 
order  of  1000  M-1  cm-1,  a  known  formation  rate  constant  of 
4.1  x  108  M-1  s~  ,  and  a  broad  absorption  band  centered  on 
600-610  nm.138  The  equation  representing  the  competition  is 

k2[co32-] 

[C03  ]inf#  =      _ — -   x   [0H.]o   IH-4 

k2[C032  ]  +  k^TPB  ] 


with  [C03~]inf#  =  end  of  pulse  concentration  of  C03~  ion 
[OH' ]Q     =   initial  concentration  of  OH'  (includes 
primary  OH'  from  water  radiolysis  and 
from  interconversion  of  ea  ~  by  N20) 

Taking  [0H«]o  as  equivalent  to  the  [C03~]inf<  generated  with 
no  TPB"  present,  or  [C03~]o,  equation  III-4  can  be  rear- 
ranged to 

[C03-]o  k^TPB-] 

1     +      III-5 

[C03-]inf.  k2[C032  ] 

From  a  plot  of  [C03~]o  /  [C03"]inf#  versus  [TPB~]  /  [C032~], 
as  shown  in  Figure  25,  k1/k2  is  found  to  be  15.   Substitu- 
tion of  k2,  the  C03~  reference  formation  rate  constant  gives 
a  rate  constant  for  the  formation  (reaction  III-3  above)  of 
the  initial  transient  =  (6.2  ±  0.6)  x  109  M_1sec-1. 


30 


(TPB  )  /  (CO*  ) 


0.25 


Fig.  25.   Carbonate  competition  experiment  for  determining 
the  OH  +  TPB-  rate  constant. 
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Radical  formation  rate  constants  for  TPB*  equal  to  1.4  x  109 
M-1s-1  and  7.1  x  108  M-1  s-1  have  been  reported  for  the 
electron  transfer  reaction  involving  TPB"  and  N3«.75'76 

3.1.4  Reduction  Intermediates  of  Tetraphenylborate  Ions 

Pulse  radiolysis  experiments  were  performed  using 
reducing  conditions  to  examine  the  possibility  of  a  rapid 
reaction  between  the  TPB"  ion  and  ea  ~ .   Because  the  TPB- 
anion  is  a  strong  reducing  agent  and  negatively  charged, 
this  reaction  was  presumed  unlikely.   However,  by  monitoring 
the  ea  ~  transient  signal  as  a  function  of  added  TPB-  any 
such  reaction  would  be  evinced  by  either  a  decrease  in  the 
magnitude  of  the  end  of  pulse  ea  ~  transient  absorption  or  a 
change  in  the  ea  ~  transient  decay  kinetics.   A  typical 
transient  trace  of  the  pulse  experiments  conducted  on  ap- 
proximately 10  mM  NaTPB  solutions  under  reducing  conditions 
is  shown  in  Figure  26.   The  life-time  of  the  ea  ~  transient 
was  found  to  be  approximately  8-10  microseconds  with  a  best 
fit  analysis  giving  a  second  order  decay  rate  constant  of 
5.7  x  1010  M~1s"1. 

3.1.5  Computer  Simulation 

Modeling  calculations  were  utilized  to  help  establish 
the  detailed  chemical  mechanism  of  multi-step  reactions 
found  in  aqueous  pulse  radiolysis  systems.   Differential 
rate  equations  for  the  suspected  reaction  sequence  are 
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Fig.  26.   a)  Graphical  display  of  optical  signal  versus  time 
following  pulse  radiolysis  of  an  N2~saturated  10 
mM  NaTPB,  0.04  M  t-butanol  solution;  20  nsec 
period,  10  usee  full  scale;   b)  Curve  fit  of  data 
between  x  =  3  and  x  =  7  microseconds. 
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integrated  to  give  computer-simulated  time-dependent  concen- 
tration profiles  that  are  compared  with  those  measured 
experimentally.   These  comparisons  help  to  establish,  or 
validate,  proposed  mechanisms.   The  computer  modeling  was 
done  using  a  Gear  integrator  program.     This  program, 
written  by  R.  L.  Brown140  of  NIST,  was  modified  in  this 
laboratory  to  run  under  Microsoft  FORTRAN  77  on  an  MS-DOS- 
based  microcomputer. 

The  input  parameters  to  the  program  consist  of  a  list- 
ing of  all  mechanistic  steps  and  their  rate  constants  (see 
Appendix  B) ,  the  identification  of  all  reactants  with  their 
initial  concentrations,  and  the  time-step  which  controls  the 
50  computation  intervals. 

3.2   Discussion 

3.2.1  Electron  Transfer 

The  oxidation  intermediates  of  TPB"  ions  in  agueous 
solutions  have  been  examined  by  Horii  and  Taniguchi,75  and 
Liu  et  al. ,6   From  their  results  with  N20-saturated  solu- 
tions of  NaTPB  and  NaN3 ,  the  proposed  route  to  the  formation 
of  several  major  products  is  a  rapid  electron  transfer  from 
B(C6H5)4~  to  N3'  yielding  the  B(C6H5)4-  radical  and  N3~,  or 

B(C6H5)4~  +  N3*  >  B(C6H5)4-  +  N3_      III-6 

As  shown  earlier  the  azidyl  radical,  N3 * ,  is  produced  ini- 
tially from  primary  OH  radical  reaction  with  N3~,  reaction 
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III-l.   The  transient  absorbing  at  3  30  nm  and  800  nm  (Figure 
21)  is  attributed  to  the  B(C6H5)4'  radical.   Both  peaks 
decay  with  the  same  first-order  rate  constant,  indicating 
that  they  represent  different  absorbances  of  the  same  tran- 

•  •     7fi  ■         • 

sient  species.  "     Loss  of  the  330  nm  peak  is  attributed  to  a 
rapid  self -decomposition  first-order  decay,  producing  a 
second  transient  that  undergoes  a  relatively  slower  decay. 
Horii  and  Taniguchi  proposed  the  reaction,75 

B<C6H5V >  B(C6H5)2'  +  (C6H5)2      1J1-7 

The  secondarily-produced  transient  B(C6H5)2*  was  suggested 
to  disappear  following  second-order  kinetics,  although  no 
product  species  were  discussed.75 

The  pulse  results  from  our  oxidation  experiments  on  TPB 
solutions  containing  no  added  N3~  anion  could  be  interpreted 
as  above,  i.e.  a  rapid  electron  transfer  from  B(C6H5)4~  to 
OH*,  instead  of  N3«,  or 

B(C6H5)4"  +  OH'  >  B(C6H5)4*  +  OH"    III-8 

yielding  the  B(C6H5)4*  radical.   As  shown  in  Table  2,  the 
decay  kinetics  of  the  transients  produced  with  and  without 
added  N3~  are  similar,  that  is  a  rapid,  first-order  decay 
followed  by  a  secondary  transient  decay.   However,  compari- 
son of  the  transient  spectra,  Figure  21,  suggests  that  two 
different  transients  are  formed.   With  added  azide,  the 
transient  spectral  feature  at  750  -  800  nm  is  about  70%  of 


the  magnitude  of  the  main  peak  centered  at  3  30  nm.   With  no 
added  azide,  the  transient  spectral  feature  at  750  -  800  nm 
is  <  10%  of  the  magnitude  of  the  main  peak  centered  at  3  30 
nm. 

3.2.2  OH  Addition 

An  alternate  mechanism  to  the  electron  transfer  from 
TPB-  to  OH*,  as  described  above,  is  the  addition  of  the  OH" 
radical  to  the  aromatic  ring,  which  would  lead  to  an  OH' 
adduct  transient  species  different  from  the  TPB*  radical,  or 


B(C6H5)4   +  0H'  >  (C6H5)3BC6H5OH 
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The  following  self-decomposition  first-order  decay 
scheme  has  been  considered  for  the  observed  loss  of  this  OH' 
addition  adduct. 


(C6H5)3BC6H5OH'-  >  B(C6H5)3  +  C6H5OH* 

C6H50H-~  +  H20  >  C6H6OH"  +  OH" 


302 
303 


Triphenylborate  and  its  hydrolysis  products,  (C6H5)B(OH)2 
and  (C6H5)2BOH,  all  have  absorptions  at  wavelengths  below 
300  nm.74   The  secondarily-produced  transient  C6H5OH'~,  is 
suggested  to  undergo  a  rapid  hydrolysis,  yielding  the  hy- 
droxycyclohexadienyl  radical  C6H6OH« .   The  hydrolysis  of 
C6H50H  •  has  been  reported  in  a  previous  aqueous  phenol 
pulse  radiolysis  study.141   In  neutral  solutions  a  slow 
reaction  between  the  solvated  electron,  e  ~,  and  phenol 
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that  forms  C6H5OH~*  was  reported.   The  radical  anion, 
C6H5OH-',  subsequently  undergoes  rapid  protonation  by  H20  to 
yield  (C6H5) (OH)H* .   This  species  is  a  structural  isomer  of 
the  hydroxycyclohexadienyl  radical,  C6H6OH* .   A  molar  ex- 
tinction coefficient  of  (3800  +/-  800)  M"1  cm-1  at  a  wave- 
length maximum  of  330  nm  was  assigned  to  the  (C6H5)(OH)H* 
radical . 

Another  possible  competing  reaction  pathway  for  the 
C6H5OH*~  radical  anion  would  be  a  bimolecular  radical  anion 
reaction.   The  bimolecular  disappearance  of  negatively 
charged  radical  species  has  been  previously  suggested  in 
pulse  radiolysis  studies  involving  the  reaction  of  the 
hydrated  electron,  ea  ~,  with  benzene.58   In  this  aqueous 
benzene  study,  C6H6~,  the  addition  product  of  ea~  and 
benzene,  was  reported  to  decay  via  a  rapid  bimolecular  proc- 
ess before  protonation  in  solutions  of  pH  neutral  to  13. 

C6H6  +  eaq~  >  C6H6~  IIZ"9 

C6H6"  +  C6H6~  >  Products  111-10 

No  specific  products  of  the  radical  anion  reaction,  111-10, 
were  discussed. 

In  the  present  study,  concerning  the  combination  reac- 
tion of  C6H5OH*~,  a  possible  reaction  forming  a  charged 
dimer  species,  with  subseguent  formation  of  biphenyl  and 
hydroxide  ion,  is  proposed  below 

■ 
C6H5OH*"  +  C6H5OH'"  >  [HO(C6H5)-(C6H5)OH]2"    304 

[HO(C6H5)-(C6H5)OH]2"  >  C6H5-C6H5  +  2  OH"     305 
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The  scheme  involving  radical  combination  leading  to  an 
intermediate  dimer  is  prominent  in  the  mechanism  involving 
the  radiolysis  of  agueous  benzene,  discussed  below. 

The  hydroxycyclohexadienyl  radical,  C6H6OH' ,  mentioned 
above  (reaction  303),  is  an  important  intermediate  in  the 
formation  of  major  products  in  the  radiolysis  of  agueous 
benzene  solutions;  the  results  of  which  are  in  good  agree- 
ment with  the  end  products  derived  from  Co-60  steady-state 
r-radiolysis  of  NaTPB  solutions  carried  out  in  this  lab.86 
The  proposed  mechanism  concerning  the  benzene 
oxidation55'64'65  is  described  as 

C6H6  +  OH'  >  C6H6OH*  III-9 

2  C6H60H-  >  [HO(C6H6)-(C6H6)OH]        306 

[HO(C6H6)-(C6H6)OH]  >  C6H5-C6H5  +  2H20   307 

2  C6H6OH-  >  C6H6  +  C6H5OH  +  H20        308 

The  hydroxycyclohexadienyl  radical,  C6H6OH- ,  has  a  reported 
absorption  maxima  at  313  nm.58   The  combined  bimolecular 
rate  constant  for  its  disappearance  via  reactions  306  and 
308  has  been  determined  to  be  (7.05  ±  0.35)  x  108  M"1s-1.58 
The  disproportionation  to  combination  ratio  for  the  disap- 
pearance of  C6H6OH*  was  found  to  be  <  0.4,  indicating  that 
the  combination  reaction  306  is  predominant.63   The  major 
products  produced  are  biphenyl,  benzene  and  phenol.   In  the 
Co-60  r-radiolysis  of  aerated  0.05  M  sodium  tetraphenylbo- 
rate  solutions  the  products  and  their  G-values  were  found  to 
be  biphenyl,  1.9;  benzene,  1.3;  and  phenol,  0.41.68 
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It  is  also  reasonable  to  expect  the  H«  atom  to  add  to  TPB' 
ion,  giving  a  transient  species  that  decomposes  via  a  two- 
step  process  giving  the  cyclohexadienyl  radical 


B(C6H5)4"  +  H-  >  (C6H5)3BC6H6T  309 

(C6H5)3BC6H6T  >  B(C6H5)3  +  C6H67       310 

C6H6~  +  H2°  >    C6H7*  +  0H~  311 

The  (C6H5)3BC6H6«  transient  has  been  proposed  by  Liu  et  al. 
as  the  weakly  absorbing  species  at  higher  wavelengths 
(Figure  21)  in  the  pulse  radiolysis  of  TPB-  solutions  with- 
out added  azide.76   The  cyclohexadienyl  radical  could  also 
react  with  the  hydroxycyclohexadienyl  radical  in  a  bimolecu- 
lar  process  producing  benzene56 

C6H7*  +  C6H60H*  >   2  C6H6  +  H2°  312 

Concerning  the  reduction  of  TPB",  a  rapid  reaction 
between  eaq~  and  TPB"  appears  to  be  negligible  from  the 
results  of  the  NaTPB  pulse  radiolysis  under  reducing  condi- 
tions.  Comparison  of  Figure  26  with  e   ~  traces  found  in 
the  hydrated  electron  dosimetry  experiments,  Appendix  A, 
shows  little  difference  in  either  the  magnitude  or  kinetic 
behavior  of  the  eaq~  transient  signal.   The  e   ~  transient 
signal  decayed  with  a  second  order  rate  constant  of  5.7  x 


10 


10  „-U-l 


M"±s~1  with  added  TPB";  for  the  e   ~  dosimetry  experi- 
ments, the  eaq~  signal  exhibited  a  second-order  decay  =4.4 


x  1010  M^s"1 
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3.2.3  Computer  Simulations 

Computer  simulations  were  carried  out  with  a  Gear 
integrator  utilizing  the  data  set  consisting  of  the  elemen- 
tary reaction  steps  describing  the  transients  of  irradiated 
pure  water  (Appendix  B) .   The  mechanistic  scheme  also  in- 
cludes reactions  of  e  ~  with  N20,  and  OH*  with  TPB-. 

Figure  27  illustrates  typical  results  of  the  simulated 
calculations  for  conditions  involving  N20- saturated  solu- 
tions containing  10  mM  TPB-.   Time-concentration  profiles 
for  all  species  in  the  system  are  predicted  over  the  ca.  200 
microsecond  data  collection  time  period.   Experimental  data 
points  indicate  the  observed  decay  of  the  transient  trace 
produced  in  a  typical  experimental  run,  recorded  at  3  30  nm 
(Figure  20  (b)).  These  experimental  data  points  represent 
the  combined  absorbance  of  all  the  components  in  the  system 
contributing  to  the  absorption  trace  at  3  30  nm.   Dotted 
lines  show  predicted  yields  of  the  various  individual  tran- 
sients formed.   The  solid  line  represents  the  combined 
concentrations  of  the  following  intermediates;  the  initial 
OH*  adduct  (C6H5)3BC6H5OH"  with  a  molar  extinction  coefi- 
cient  of  10,000  M"1  cm"1,  the  radical  anion  C6H5OH~  and  the 
hydroxycyclohexadienyl  radical  C6H6OH*  with  molar  extinction 
coefficients  of  3800  M-1cm-1,  and  the  intermediate  dimers, 
[HO(C6H6)-(C6H6)OH]  and  [HO(C6H5 )-(C6H5 )OH] 2"  with  molar 
extinction  coef f iecients  of  7,500  M-1cm   . 
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Fig.  27.   Comparison  of  experimental  and  computer-simulated 
transient  concentrations  versus  time  plots  during 
the  pulse  radiolysis  of  a  10  mM  NaTPB  solution; 
( — )  summation  of  individual  transients; 
(•)  typical  experimental  results; 
( )  computed  individual  transient  traces. 
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3.3  Conclusions 

Results  obtained  in  this  work  lead  to  the  conclusion 
that  the  identity  of  the  transient  species,  absorbing  upon 
pulse  radiolysis  of  N20-saturated  NaTPB  solutions,  is  very 
much  dependent  on  the  nature  of  the  oxidizing  species  pro- 
duced in  the  system.   The  aqueous  azide  NaTPB  solution 
experiments  presented  here,  taken  together  with  the  results 
of  others,75,  6  indicate  that  the  one  electron  oxidation 
product  of  TPB-,  the  TPB*  radical,  is  the  species  that 
absorbs  at  3  30  and  775  nm. 

However,  experiments  in  aqueous  NaTPB  solutions  without 
added  azide  indicate  formation  of  a  different  transient 
species,  proposed  to  be  the  OH"  adduct  of  TPB~ ,  or 
(C6H5) 3B(C6H5)OHT.   This  OH*  addition  product  exhibits  a 
strong  absorption  at  330  nm.   A  weak  absorption  at  750  nm 
constitutes  approximately  10%  of  the  magnitude  of  the  main 
3  30  nm  peak.   As  mentioned  earlier,  this  weak  absorbance 
could  be  due  to  the  H*  addition  product  of  TPB~,  since  the 
addition  of  H*  to  aromatic  systems  with  resonating  7r-bond 
systems  has  been  observed.   Alternatively,  the  weak  750  nm 
absorption  could  be  indicative  of  the  TPB*  radical  product 
formation,  which  would  indicate  that  the  one  electron  trans- 
fer mechanism  occurs  to  a  small  extent  in  the  N20-saturated 
TPB  solutions  without  added  azide. 

The  present  study  provides  information  regarding  the 
fate  of  the  TPBOHT  adduct  immediately  following  its  forma- 
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tion,  whereas  the  previous  account  mentioned  only  alluded  to 
the  TPBOHT  adduct  transient  absorption  spectrum.76   Kinetic 
analysis  of  this  transient  species  indicates  that  it  under- 
goes a  rapid  first-order  disappearance  within  the  first  20 
microseconds  that  is  suggested  to  be  a  self -decomposition. 
It  is  presumed  that  this  initial  transient  loss  leads  to  the 
observed  secondary  absorption,  which  was  found  to  decay  on  a 
longer  timescale  by  second  order  decay  kinetics. 

A  primary  interest  in  this  work  was  to  establish  a 
plausible  reaction  mechanism  involving  the  radiolysis  of  TPB 
solutions.   Evidence  obtained  from  the  pulse  radiolysis 
experiments  has  lead  to  the  formulation  of  a  scheme  that  can 
explain  the  production  of  various  organic  products  such  as 
biphenyl,  benzene,  and  phenol,  reported  from  previous 
steady-state  radiolysis  studies  on  agueous  aerated  TPB 
solutions  sealed  from  the  atmosphere.68   For  further  valida- 
tion of  the  mechanism  presented  in  this  study  for  the  pulse 
radiolysis  conditions,  one  could  perform  numerical  integra- 
tion of  the  scheme  for  both  the  pulse  and  steady-state 
radiolysis  systems  under  eguivalent  solution  conditions. 
Initial  attempts  to  carry  out  this  work  have  been 
encouraging. 


CHAPTER  IV 

THE  PULSE  RADIOLYSIS  OF  IRIDIUM  HEXACHLOROIRIDATE 

IN  AQUEOUS  SOLUTION 

4.1   Results 


4.1.1  Absorption  Spectra 

Initial  experiments  performed  with  the  IrCl63~  and 
IrClg   ions  involved  measurement  of  the  aqueous  solution 
absorption  spectra  of  these  complexes.   These  absorption 
spectra,  Figure  28,  of  millimolar  aqueous  solutions  of  the 
IrCl6J   and  IrClg   complex  ions  were  in  agreement  with  the 

results  of  other  authors.111'142   Between  400  and  500  nm, 

P— 
IrClg    has  a  characteristic  spectrum  with  an  absorption 

maximum  at  488  nm  (extinction  coefficient  ca.  4000).   In 

contrast,  the  IrClg   complex  has  no  appreciable  absorption 

in  this  wavelength  range.   Thus  all  pulse  absorption  data 

were  recorded  at  488  nm. 

4.1.2  N20-Saturated  Iridium(III)  Solutions 

The  reaction  of  0H«  radicals  with  IrCl63~  occurs  too 
rapidly  to  be  observed  directly  under  our  experimental 
kconditions.   A  method  involving  competition  kinetics,  the 
carbonate  competition  method  described  earlier,  was  used  to 
determine  the  rate  constant  for  the  reaction 

OH*  +  IrCl63~   =  0H~  +  IrCl62~  201 
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A  series  of  experiments  was  performed  using  N20- 
saturated  solutions  containing,  typically,  1  x  10-3  M 
IrCl63-  and  various  amounts  of  C032-  from  2  to  14  x  10~3M. 
Under  these  conditions  the  hydrated  electron,  e  ~,  is 
converted  to  OH* .   The  eguations  representing  the  competi- 
tion are 

0H«  +  C032~   >  OH"  +  C03"  130 

0H«  +  IrCl63-   >  OH"  +  IrCl62~  201 

The  end-of-pulse  Ir(IV)  absorption,  Ir(IV)Q,  was  first 
monitored  for  a  millimolar  solution  of  Ir(III)  to  obtain  the 
maximum  Ir(IV)  produced,  i.e.  without  any  added  C032~.   The 
end-of-pulse  Ir(IV)  absorption,  Ir(IV)^nf  ,  was  then  ob- 
tained for  a  series  of  solutions  with  added  C032~.   A  plot, 
Figure  29,  of  [ Ir( IV)Q ]/[ Ir ( IV)inf _ ]  versus  the  added  C032~ 
concentration  yields  a  slope,  k13O/(k201 [Ir(III) ] ) ,  egual  to 
85  M-1.   Substitution  of  the  Ir(III)  concentration  of  1 
mM  and  k130  =  4.1  x  10°  into  the  slope  expression  above 
gives  a  rate  constant  for  reaction  130  above,  egual  to  (4.9 
±  0.6)  x  109  M~-*s"  •.   Broszkiewicz  measured  a  somewhat 
higher  value  of  8.9  x  109  M~1s~1  for  this  reaction.111 
However,  the  lower  value  gives  satisfactory  results  in  the 
computer  simulation  studies. 

In  all  the  experiments  on  the  pulse  radiolysis  of  C032~ 
/IrCl63~  solutions,  using  concentrations  of  IrCl63-  of  1  x 
10~3  M  and  concentrations  of  C032~  from  2  to  14  x  10~3  M, 
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Fig.    29.      Carbonate  competition  experiment   for  the  determi- 
nation of  the  OH  +   Ir(IV)    rate   constant. 
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the  product,  IrClg2-,  disappears  rather  rapidly  as  shown 
in  Figure  30  (a).   If  the  decay  is  fitted  as  a  simple  sec- 
ond-order process  in  IrCl62~,  an  apparent  rate  constant  of  4 
x  107  M-1s-1  is  obtained.   However,  loss  of  the  Ir(IV)  via  a 
bimolecular  disproportionation  reaction  is  not  reasonable 
mechanistically.   An  interpretation  of  the  cause  of  the 
decay  of  the  IrCl62~  optical  signal  in  these  experiments  is 
discussed  below.   Interestingly,  as  seen  in  Figure  30  (b) , 
the  product  iridium  (IV)  complex  is  almost  inert,  or  stable 
over  approximately  400  microseconds,  when  the  experiment  is 
performed  with  added  N20  using  a  neutral  medium,  i.e.  with 
no  added  NaOH  or  Na2C03 . 

4.1.3  Nitrogen  and  Air-saturated  Ir(III)  Solutions 

Experiments  were  performed  on  the  pulse  radiolysis  of 
both  nitrogen  and  air-saturated  Ir(III)  agueous  solutions  to 
further  investigate  the  reactivity  of  the  Ir(IV)  species. 
Using  N2-saturated  solutions,  roughly  egual  amounts  of  the 
reducing  and  oxidizing  species,  ea  ~  and  OH*  respectively, 
are  produced.   The  N2 -saturation  purges  all  dissolved  atmos- 
pheric oxygen,  an  efficient  electron  scavenger,  from  the 
solution,  thus  eliminating  the  reaction 

eag~  +  °2  >  °2~  113 

In  the  pulse  radiolysis  of  air-saturated  solutions,  which 
contain  02  at  0.25  mM,  the  superoxide  anion,  02~,  plays  an 
important  role40  due  to  its  formation  from  the  reaction 
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Fig.  30.   Graphical  display  of  IrCl62~  optical  signal  versus 
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above,  and  the  reaction 

H  +  02 >  H02  117 

With  a  pK  of  4.88  for  H02 ,  02~  is  the  predominant  species 
above  pH  =  5.   The  perhydroxyl  radical,  H02 • ,  or  02~  in 
basic  media,  can  act  as  an  oxidizing  or  reducing  agent 
depending  on  the  solute  present. 

The  transient  Ir(IV)  oxidation  product  was  found  to  be 
unstable  in  both  N2-  and  air-saturated  solutions,  as  shown 
in  Figure  31.   The  spectra  of  Ir(IV)  disappeared  rapidly  in 
the  N2 -saturated  solutions  with  a  second-order  rate  constant 
of  5  x  10°  M-1s-1  and  moderately  rapid  with  a  second-order 
rate  constant  of  1  x  107  in  the  case  of  air-saturated  solu- 
tions.  These  rate  constants  establish  the  time-scale  of  the 
phenomena.   Additionally,  since  concentrations  of  the  oxi- 
dizing and  reducing  reactants  (originating  from  OH  and  e   ~ 
respectively)  are  expected  to  be  nearly  eguivalent,  the  rate 
constants  should  correctly  describe  the  actual  reactions. 

4.1.4  Reduction  of  Irdll^ 

The  rate  constant  for  the  reduction  of  the  IrCl63~ 
complex 

IrCl63"  +  eaq~   =   IrCl64"  202 

was  obtained  by  applying  pseudo-first-order  kinetics. 
Formation  of  the  hydrated  electron  in  water  was  observed 
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Fig.  31.   Graphical  display  of  IrCl62~  optical  signal  versus 
time  following  pulse  radiolysis  of  neutral  1  mM 
IrCl63-  solutions  saturated  with  (a)  N2;  (b) 
air. 
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using  a  50  MHz  sampling  rate  with  our  transient  digitizer 
and  its  decay  was  followed  at  632.8  nm.   The  rate  expression 
for  reaction  202  above  is 

-d[eaq"]/dt  =  k202  [Ir(III)]  [eaq_]        IV-l 

By  making  the  concentration  of  Ir(III)  high  compared  to  the 
expected  concentration  of  ea~ ,  k2Q2  [Ir(Ii:c)]  i-s  treated 
as  a  constant  and  the  rate  expression  reduces  to 

-d[eaq"]/dt  =  k'  [eaq_]  IV-2 

where  k  =  k202  [Ir(III)]»   Expression  IV-2  is  a  first-order 
rate  expression.   The  first-order  decay  rate  constant  for 
e   ",  k  ,  was  determined  using  the  appropriate  fit  routine 

"Hi 

in  the  BASIC  data  analysis  program  found  in  Appendix  C. 
These  first-order  rate  constants  for  the  ea  ~  decay  were 
determined  and  plotted  versus  the  concentration  of  added 
Ir(III)  complex,  Figure  32.   The  nitrogen-saturated  solu- 
tions were  made  0.1  M  in  t-Butanol  to  remove  OH*  radicals. 
The  value  for  the  rate  constant  corresponding  to  reaction 
202  above  was  found  to  be  (6.1  ±  0.3)  x  109  M^s"1. 

4.2   Discussion 

The  results  presented  above  agree  in  general  respects 
with  those  of  Broszkiewicz;  OH*  radicals  react  rapidly  with 
IrCl63~  to  produce  the  (IV)  oxidation  state  of  iridium  via 
reaction  201  (Table  5).   Our  measured  rate  constant  is 


102 


1.6x10 


1.5x10 


O 
CD 

CO 


1.4  X  10 


CO 

n 

J£  6 

1.3X10 


1  2x10 


2  3 


4  5  6 

[lr(lll)]x105,M 


Fig.  32.   Change  in  first-order  decay  of  the  hydrated  elec- 
tron with  IrCl63_  concentration. 
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slower  than  that  seen  by  Broszkiewicz,  as  mentioned  earlier. 
Broszkiewicz  also  noted  that  IrCl62~  is  "stable  within  the 
period  of  ca.  0.5  ms,  though  its  relative  instability,  when 
compared  with  solutions  of  IrCl62~  as  prepared  from  solid 
salt  indicates  that  radiation  produces  either  an  active  form 
of  IrCl62~  of  higher  reactivity  towards  water,  or  a  system 
in  which  this  compound  is  more  readily  reduced."111   In 
addition,  Broszkiewicz  reported  that  IrCl62~  is  not  capable 
of  being  further  oxidized  by  any  primary  radiolysis  product, 
but  that  it  is  capable  of  being  reduced  within  a  few  micro- 
seconds. 

ill 
In  agreement  with  the  earlier  work  mentioned,  -LJ-  it 

was  found  that  the  IrCl62_  species  is  stable  when  formed  in 

neutral  solutions  saturated  with  N20,  which  converts  eaq~  to 

OH-  radical.   In  contrast,  however,  the  Ir(IV)  complex  is 

rather  short-lived  under  our  circumstances  in  either  of 

three  conditions:  1)  when  formed  in  a  neutral  N2-saturated 

solution,  2)  when  formed  in  an  air-saturated  solution  in  the 

absence  of  N20  (i.e.,  with  both  ea  ~  and  OH*  transients 

present,  and  3)  when  formed  in  an  alkaline,  N20-saturated 

solution. 

4.2.1  Loss  of  Iridium  (IV)    in  Nitrogen-Saturated  Solutions 

As  noted  in  the  earlier  work,111  iridium  (III)  is 
readily  reduced  to  iridium  (II)  by  either  hydrated  electrons 
via  reaction  202  or  by  hydrogen  atoms  via  reaction  203 
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IrCl<53-  +  eaq   -   IrCl64-  202 

IrCl63-   +   H-     =   IrCl64"   +    H+       203 

Our  measured  rate  constant  for  the  reduction  of  IrCl63~  by 
reaction  202  is  in  good  agreement  with  earlier  reported 
values.110'111 

Broszkiewicz  suggested  that  the  oxidation  state  (II) 
compound  disappears  rapidly  via  disproportionation  to 
Ir(III)  and  Ir(I).   While  no  evidence  was  found  contrary 
to  this  possibility,  an  alternate  suggestion  is  discussed 
below.   N9-saturated  solutions  pulse-irradiated  in  the 
absence  of  N20  contained  Ir(II)  formed  by  reactions  202  and 
203.   In  addition,  Ir(IV)  formed  from  reaction  201  was 
produced  under  these  conditions.   Thus  it  is  suggested  that 
under  these  circumstances,  the  most  likely  route  of  disap- 
pearance of  Ir(II)  is  an  electron  transfer  involving  Ir(IV) 

IrCl62-   +   IrCl64~  -   2  IrCl63"  206 

Interestingly,  this  reaction  was  proposed  by  Dainton108  in 
connection  with  steady-state  radiolysis  of  iridium  (III) 
systems  at  natural  pH,  as  studied  by  endproduct  analysis. 

4.2.2  Loss  of  Iridium  fIVt  in  Aerated  Solutions 

Air-saturated  solutions  showed  a  similar  decay  of  the 
IrCl64~  absorption,  although  to  a  lesser  extent  than  the  N2~ 
saturated  solutions.   However,  in  the  aerated  solutions 
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millimolar  concentrations  of  dissolved  02  (an  effective  ea  ~ 
scavenger)  were  in  competition  with  Ir(III)  for  ea  ~  and  H« 
atoms.   Thus  only  a  portion  of  the  H«  atom  and  ea  ~  yields 
reacted  via  processes  202  and  203.   The  superoxide  anion, 
02~,  formed  from  reactions  113  and  117,  can  act  as  either  an 
oxidant  or  reductant  depending  on  the  solute  present.   The 
IrCl62~  complex  has  a  standard  reduction  potential,  E°,  = 
+1.02  V.143   Thus  the  02~   anion,  with  a  standard  reduction 
potential  of  -0.33,143  is  identified  as  the  reductant  of 
IrCl62~  in  the  aerated  system. 

4.2.3  Loss  of  Iridium  (IV)  in  Basic,  N20-Saturated  Solutions 

As  indicated  previously,  our  experimental  conditions 
differ  somewhat  from  those  used  by  Broszkiewicz  in  that 
there  exists  a  much  higher  instantaneous  concentration  of 
OH*  radicals  leading  to  net  formation  of  H202.   As  suggested 
by  Broszkiewicz,  there  is  no  obvious  method  to  oxidize  the 
iridium  (IV)  compound,  but  its  disappearance  by  reduction  is 
a  clear  possibility.   For  this  disappearance  to  occur,  it  is 
necessary  that  some  other  component  in  the  solution  be 
oxidizable.   After  elimination  of  all  other  obvious  possi- 
bilities, it  was  concluded  that  the  only  plausible  reductant 
is  hydrogen  peroxide,  or  actually  the  H02~  anion  in  basic 
solution. 
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As  noted  by  several  workers,144-146  iridium  (IV)  chlo- 
ride will  ultimately  decompose  in  alkaline  aqueous  solution, 
with  formation  of  hydrated  species.   However,  this  process 
requires  hundreds  of  hours  and  is  far  too  slow  to  be  rele- 
vant to  our  pulse  radiolysis  conditions.  Fine147  reports  a 
rapid  and  quantitative  reduction  of  iridium  (IV)  chloride  in 
aqueous  solutions  at  pH  >  11  as  a  necessary  preliminary  step 
to  any  subsequent  hydroxylation  reactions.   In  concentrated 
alkaline  solution  and  in  the  absence  of  other  reducinq 
aqents,  a  reaction  with  hydroxide  ion  was  suqqested  by 
Fine147 

2  IrCl62~  +  2  0H~  =  2  IrClg3"  +  0.5  02  +  H20       IV-3 

However,  in  view  of  the  rate  constants  reported  by  Fine, 
this  reaction  is  too  slow  to  explain  our  results  obtained  at 
the  pH  values  of  our  solutions  (pH  ca.  10).   It  was  con- 
firmed by  simple  benchtop  experiments  that  the  decay  of  the 
characteristic  IrCl62-  spectrum  in  millimolar  alkaline 
solution  is  very  slow. 

Dainton108  proposed  that  iridium  (IV)  is  capable  of 
oxidizinq  hydroqen  peroxide  accordinq  to  the  net  reaction 

H202   +   2  IrCl62~  =  2    IrCl63-  +   02  +   2  H+       IV-4 

This  reaction  has  also  been  mentioned  by  Waltz  and 
Adamson.101   However,  it  was  found  in  the  present  study  that 
the  reaction  of  H202  with  IrCl62~  in  neutral  solution  is 
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only  moderately  rapid,  requiring  typically  ten  to  fifteen 
minutes  to  produce  a  substantial  color  change  in  neutral, 
millimolar  solutions  of  each  reagent. 

In  contrast  with  the  above  observations ,  it  was  found 
that  if  millimolar  solutions  of  H202  and  IrCl62"  were  sud- 
denly also  made  ca.  1  millimolar  in  0H~  by  adding  three 
drops  of  0.1  molar  NaOH  to  20  ml.  of  solution,  the  color 
change  indicating  transformation  of  Ir  (IV)  to  Ir  (III)  was 
essentially  instantaneous;  that  is,  the  color  changed  as 
rapidly  as  the  solution  could  be  mixed  by  swirling.   These 
results  suggest  that  the  redox  process  that  causes  a  de- 
struction of  the  IrCl62"  optical  signal  in  the  experiments 
mentioned  above,  involves  the  peroxide  anion 

IrCl62~    +   H02~   =  IrCl63_   +  H02-  208 

The  peroxide  anion,  H02~,  can  act  as  either  an  oxidant  or 
reductant  depending  on  the  solute  present.   As  mentioned 
earlier,  the  IrCl62-  complex  has  a  standard  reduction  poten- 
tial, E°,  =  +1.02  V.142   Thus  the  H02~  anion,  with  a  stand- 
ard reduction  potential  of  -0.15  V,143  is  attributed  as  the 
reductant  of  IrCl62_  for  the  conditions  described  above. 
The  hydroperoxy  free  radical  is  consumed  by  several  subse- 
quent processes,  in  particular 


H02«   +  H02«    =   H202   +  02  111 
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Also,  in  solutions  of  pH  >  5,  H02 •  is  converted  into  the 
reducing  species  02~  which  could  also  contribute  to  the  loss 
of  the  IrCl62~  optical  signal 

IrCl62-  +  02"   =  IrCl63"  +  02  207 

Reaction  208  was  confirmed  directly  by  stopped-flow  experi- 
ments using  a  Durrum  Model  110  instrument  with  a  Nicolet 
digital  oscilloscope.   A  rate  constant  of  approximately  1  x 
106  M-1s_1  was  found  for  this  process.   Since  the  time 
response  of  the  stopped-flow  apparatus  is  limited  by  a 
finite  mixing  time,  this  value  may  represent  a  lower  limit 
on  the  rate  constant.  Computer  simulations  of  our  pulse 
radiolysis  data  indicate  a  rate  constant  of   4.7  x  108  M-1 
s_1  for  this  reaction. 

4.2.4  Computer  Simulations 

The  consistency  of  the  above  postulates  has  been  exam- 
ined by  carrying  out  computer  simulations  with  a  Gear  inte- 
grator utilizing  the  data  set  consisting  of  the  elementary 
reaction  steps  describing  transients  of  irradiated  pure 
water  (Appendix  B).   The  mechanistic  scheme  also  includes 
reactions  of  ea  ~  with  N20  and  with  IrCl63~,  and  of  OH* 
radical  with  IrCl63~,  or  with  C032-  (if  present). 

Figure  3  3  shows  a  typical  result  of  the  simulation 
calculations  for  conditions  involving  5  x  10~3  M  carbonate, 
2.5  x  10~2  M  N20,  and  1  x  10"3  M  IrCl63-.   Experimental  data 
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Fig.  33.   Comparison  of  experimental  and  computer-simulated 
IrCl62~  concentration  versus  time  plots  during  the 
pulse  radiolysis  of  a  N^O-saturated  solution  con- 
taining 1  mM  IrCl63-  and  5  mM  Na2C03 .  ( — )  comput- 
ed; (•)  typical  experimental  results. 
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points  showing  formation  and  decay  of  the  IrCl62~  transient 
in  a  typical  run  are  indicated.   The  Gear  integrator  pre- 
dicts time-concentration  behavior  of  all  species  in  the 
system;   the  figure  also  illustrates  predicted  yields  of 
C03T,  C042~,  IrCl62~,  IrCl64~  and  H02~- 

Further  simulation  calculations  for  conditions  involv- 
ing 1  x  10~3  M  IrCl63-  solutions  saturated  with  air,  and 
nitrogen  are  presented  in  Figures  34  and  35  respectively. 
It  will  be  seen  that  the  agreement  between  measured  and 
predicted  behavior  for  IrCl62_  is  reasonable  but  not  perfect 
for  the  various  solution  conditions. 

The  results  of  simulated  calculations  performed  for  all 
solution  conditions  mentioned  previously  are  presented  in 
Figure  36.   These  concentration-versus-time  profiles  are  to 
be  compared  with  the  experimental  digitized  optical  traces 
of  the  IrCl62~  species  shown  in  Figures  30  and  31. 
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Fig.  34.   Comparison  of  experimental  and  computer-simulated 
IrCl62~  concentration  versus  time  plots  during  the 
pulse  radiolysis  of  a  nitrogen-saturated  solution 


containing  1  mM  IrCl6 
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Fig.  35.   Comparison  of  experimental  and  computer-simulated 
IrCl62~  concentration  versus  time  plots  during  the 
pulse  radiolysis  of  an  air-saturated  solution 


containing  1  mM  IrCl6 
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4.3  Conclusions 

Aqueous  solutions  of  IrClg-'  and  the  one-electron 
oxidation  or  reduction  products,  IrCl62~  and  IrCl64~  respec- 
tively, form  a  relatively  complex  reaction  system.   However, 
the  method  of  pulse  radiolysis  allows  for  the  collection  of 
some  experimental  evidence  on  the  early  stages  of  the  redox 
reactions  involved.   It  is  entirely  possible  there  are 
further  nuances  in  the  chemistry  of  the  iridium  complexes 
that  have  not  been  included  at  this  time.   For  instance, 
several  authors  have  suggested  the  participation  of  partial- 
ly hydrolyzed  iridium  chloride  complexes,  such  as 
IrCl5(OH2)-.148-150 

Results  from  this  work  lead  to  the  following  conclu- 
sions.  The  reported  relative  stability,  i.e.  no  evidence  of 
decay  for  the  ca.  500-microsecond  period  of  observation,  of 
the  IrCl62-  under  neutral  N20-saturated  conditions  was 
confirmed.   The  transient  IrCl62-  complex  formed  during  the 
pulse  radiolysis  of  aqueous  millimolar  IrCl63~  solutions  was 
found  to  be  unstable  under  conditions  of  neutral  air  satura- 
tion, neutral  nitrogen  saturation,  and  alkaline  nitrous- 
oxide  saturation.   For  each  of  the  solution  conditions,  an 
appropriate  species  has  been  identified  and  discussed  as  the 
most  likely  reductant  for  the  one-electron  reduction  of  the 
IrCl62~  complex.   In  the  case  of  alkaline  N20-saturated 
IrCl63-  solutions,  pertinent  benchtop  experiments,  as  well 
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as  stopped-flow  techniques,  were  utilized  to  identify  the 
peroxide  anion,  H02~,  as  the  reducing  species. 

Further  work  might  confirm  the  proposed  reaction  be- 
tween  the  Ir(IV)  and  Ir(II)  species  produced  from  the  pulse 
radiolysis  of  nitrogen-saturated  Ir(III)  solutions,  reaction 
206.   The  Ir(II)  complex  could  be  exclusively  formed  by 
pulse  irradiating,  under  totally  reducing  conditions, 
Ir(III)  solutions  containing  trace  amounts  of  Ir(IV)  added 
as  a  minor  solute.   Subsequent  analysis  and  comparison  of 
the  resulting  Ir(IV)  decay  to  the  Ir(IV)  decay  presented 
earlier  in  this  study  could  substantiate  an  electron  trans- 
fer  involving  Ir(II)  and  Ir(IV). 

Lastly,  as  mentioned  previously  the  reaction  rate 
constant  of  4.9  x  109  M-1s_:L  for  the  oxidation  of  Ir(III)  by 
OH*  was  attained  by  the  carbonate  competition  method.   This 

q 

value  is  somewhat  less  than  the  reported  value  of  8.9  x  10J 
M~1s~1.111   This  referenced  value  was  attained  by  applying 
the  kinetic  competition  method  using  t-butyl  alcohol  as  the 
competing  species,  rather  than  the  carbonate  anion.   It  is 
suspected  that  the  lower  value  reported  under  the  experimen- 
ta  conditions  in  this  study  could  be  due  to  complications 
arising  from  the  data  analysis  of  the  end-of -pulse  Ir(IV) 
concentrations .   In  the  presence  of  added  carbonate  under 
alkaline  conditions,  the  Ir(IV)  complex  undergoes  a  rapid 
decay  immediately  upon  its  formation  from  OH*  reaction  with 
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Ir(III).   However,  these  end-of -pulse  Ir(IV)  concentrations 
must  be  referenced  to  the  end-of -pulse  Ir(IV)  concentrations 
derived  from  stable  Ir(IV)  traces  that  are  obtained  when  no 
competing  species,  i.e.  carbonate,  is  present.   The  discrep- 
ancy in  the  rate  constants  could  be  explored  by  attempting 
the  t-butyl  alcohol  competition  method  for  the  reaction  of 
OH*  with  Ir(IV)  under  the  experimental  conditions  described 
in  this  work. 


APPENDICES 


APPENDIX  A 
CHEMICAL  DOSIMETRY  FOR  A  PULSED  RADIATION  SOURCE 


Chemical  dosimetry  involves  the  determination  of  the 
average  absorbed  dose  in  a  material  from  a  radiation- 
induced  chemical  change  produced  in  a  suitable  substrate. 
The  fraction  of  the  total  dose  or  flux  from  the  radiation 
source  that  is  absorbed  will  induce  the  chemical  change  in 
the  system,  and  the  most  common  unit  of  absorbed  dose  is 
the  rad,  which  is  defined  as 

1  rad  =  100  erg/g  =  6.24  x  1013  ev/g  =0.01  J/kg  A-l 

These  systems  consist  of  a  bulk  component,  which  absorbs 
the  deposited  radiation,  and  a  minor  component,  which 
reacts  with  any  radiation-induced  species  of  the  bulk  to 
produce  the  observed  chemical  change  that  is  guantitative- 
ly  measured.   Calculation  of  the  absorbed  dose  reguires 
the  knowledge  of  an  accurate  G-value  or  radiation  chemical 
yield  which  is  defined  as 

G(x)  =  n(x)/E  =  #  molecules/  100  eV    A-2 

with  n(x)  =  the  amount  of  x  produced,  destroyed,  or  changed 
by  an  absorbed  energy,  E.  The  dosimetry  system  should  have 
the  same  atomic  composition  and  density  as  the  sample  to 
be  irradiated;  thus  all  systems  presented  in  this  work  are 
agueous  solutions  that  contain  various  solutes  dissolved 
in  the  bulk  water  solvent.   Optimal  response  characteris- 
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tics  of  the  chemical  dosimeter  are  as  follows:  the  dosime- 
ter should  be  1)  proportional  to  the  radiation  dose  over 
a  large  range  that  covers  from  10  rads  up  to  100  Mrads  for 
radiation  chemistry  experiments,  2)  independent  of  both 
the  dose  rate  and  the  radiation's  energy  and  LET,  3) 
independent  of  temperature  and  4)  highly  precise,  i.e.  to 
<  5%.40 

Chemical  dosimetry  procedures  for  use  with  pulsed 
radiation  sources  are  different  from  standard  systems  used 
with  steady-state  radiolysis.   Larger  dose  rates  occur 
with  pulsed  sources  (108  to  1012  rad  sec-1)  compared  to 
continuous  dose  rates  (  10  to  10 3  rad  sec"1)  from  Co-60  r- 
ray  sources.   The  possibility  of  primary  radical-radical 
interactions,  as  opposed  to  primary  radical-solute  inter- 
actions, lead  to  a  net  loss  of  measured  product.   This 
effect  causes  the  G-value  determined  at  lower  dose  rates 
to  be  inaccurate  for  pulse  experiments.   Also,  with  the 
pulse  radiolysis  experimental  setup,  one  can  measure  tran- 
sient species  on  the  microsecond-to-millisecond  time 
scales  as  opposed  to  chemical  analysis  carried  out  after 
the  completion  of  a  series  of  radiation-induced 
reactions.40 

Multiple  chemical  dosimetry  systems  were  utilized  in 
this  work,  and  Table  3  shows  the  various  properties  in- 
volved in  the  chemical  analysis  for  each  respective  sys- 
tem.  The  first  system  applied,  using  the  1.4  ml  non-flow 
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reaction  cell,  was  a  modified  version  of  the  standard 
Ferrous  sulfate  or  Fricke  system.151   This  system  ( ImM 
Fe(NH4)2(S04) 3«6H20,  0.4  M  H2S04,  air-saturated)  involves 
the  radiation-induced  oxidation  of  ferrous  ions  to  the 
ferric  state  (Fe2+  — >  Fe3+)  in  the  presence  of  oxygen  at 
low  pH.   In  the  modified  version  (lOmM 

Fe(NH4) 2(S04) 3«6H20,  0.4  M  H2S04,  02-saturated)  recommend- 
ed for  pulse  experiments,  the  solution  is  oxygen-saturated 
and  a  relatively  higher  initial  [Fe2+]  is  used.  Detection 
of  the  product  species,  Fe3+,  was  carried  out  after  the 
completion  of  all  relevant  reactions,  approximately  thirty 
minutes,  as  is  the  practice  with  the  standard  Fricke 
system. 


Table  3 

Physical  Constants  of  the  Dosimetry  Systems 

Super  Fricke        Thiocyanate        Hydrated 

Electron 

Species 
detected 

Fe3+           (CNS)2-             eaq" 

G( product) 
#molecules 
100  eV 

16.1              5.8                  2.63 

A,  monitored 
(nm) 

304               472                632.8 
(Beckman)       (xenon-arc)          (laser) 

e,    (M-1cm_1) 

2205              7580                14860 

p,  (g/ml) 

1.024             1.001                1.012 

*   See  references  40  and  151 
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The  other  systems  used  —  thiocyanate  and  hydrated 
electron  —  involved  transient  species  detection  methods 
using  the  pulse  radiolysis  detection  instrumentation.   The 
thiocyanate  system  (10  mM  KCNS,  N20-saturated)  involves  a 
potassium  thiocyanate  aqueous  solution  that,  when  irradi- 
ated under  oxidizing  conditions,  produces  a  transient 
species,  (CNS)2~,  formed  in  a  two-step  sequence  from  the 
OH*  radical 

OH*  +  CNS"   — >  OH"  +  CNS       A-3 
CNS  +  CNS"  — >  (CNS)2~  A-4 

The  purpose  of  the  N20  saturation  is  to  convert  the  hy- 
drated electron,  ea  ~,  to  OH  that  eliminates  any  optical 
interferences  from  ea  ~  and  produces  a  larger  yield  of  the 
product  species,  (CNS)2~,  which  decays  via  a  second-order 
process  with  a  rate  constant  of  3  x  109  M_1sec-1.151 

In  the  hydrated  electron  system  (0.1  M  NaOH,  0.1  M 
ethanol,  N2 -saturated)  the  hydrated  electron,  ea  ~,  a 
primary  radical  species  produced  by  the  direct  interaction 
of  the  radiolysis  of  water,  is  monitored.   Alkaline  condi- 
tions, alcohol  addition,  and  nitrogen-purging  help  to 
minimize  the  following  reactions  of  the  primary  e. 


-aq 


eaq~  +  H+   — >  H  103 

ea  "  +  OH  — >  OH"  +  H20        102 

eaq~  +  °2  — >  °2~  113 
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The  hydrated  electron  decays  via  a  bimolecular  reaction 

eaq"  +  eaq"   ~>   H2  +  20H"    101 

with  a  rate  constant  of  2k  =  1.2  x  1010  M-1  sec-1.151 

Figure  35  shows  the  results  of  the  accumulated  ab- 
sorbance  detected  in  a  1.4  ml  volume  of  a  modified  Fricke 
solution  contained  in  the  non-flow  cell  described  earlier. 
The  absorbance  is  plotted  as  a  function  of  the  number  of 
pulses  delivered  with  the  resulting  slope  of  0.1  equal  to 
the  change  in  optical  density  per  pulse. 

Figures  36  and  37  show  the  transient  traces  of  the 
(CNS)2~  and  hydrated  electron,  ea  ~,  respectively.   The 
life-time  of  the  (CNS)2~  transient  is  shown  to  be  approxi- 
mately 30-40  microseconds  and  a  best  fit  analysis  gives  a 
second  order  decay  with  k  =  3.3  x  109  M"1sec"1  while  the  life- 
time of  the  ea  ~  transient  is  approximately  8-10  microsec- 
onds with  a  best  fit  curve  giving  a  second  order  decay 
with  k  =  4.3  x  1010  M~1sec-1. 

Using  the  experimental  AA  values  determined  along 
with  the  various  constants  collected  in  Table  3,  the 
absorbed  dose  (DD),  presented  in  Table  4,  can  be  calcu- 
lated as40 
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Dn  =  moles  product  formed  per  kg  — 


(kg) 

(molecules)        100         (eV) 
x  6.022  x  1023   x  


(mol)         G( product)   (molecule) 


x  1.602  x  10  19  x   100   A-5 


(J)                (kg  rad) 
x   100  

(eV)  (J) 


moles  product  formed  per  kg 

=   9.647  x  108  x  rads 

G ( product ) 

or, 

moles  product  formed  per  liter 

=  9.647  x  108  x   rads 

p   G( product) 

with, 

AA 
moles  product  formed  per  liter  =  

A   €   1 

then, 

A  A 

=  9.647  x  108  x    rads 

A  €  1  p  G( product) 

with, 

AA  =  difference  in  absorbance  between  irradiated  and 

non-irradiated  solution 

Ae  =  difference  in  molar  extinction  coefficient 

of  reactant  and  product  at  monitored  wavelength 

1  =  optical  pathlength,  cm 

p  =  density  of  solution,  g/ml 

G( product)  =   number  of  moles  of  product  formed 

per  100  eV  of  energy  absorbed 


124 


1.5    - 


CD 
O 

C     1 
03 

n 
o 

CO 

< 


0.5    - 


10 

Pulse  # 


15 


Fig.  37.   Change  in  absorbance  with  accumulated  radiation 
pulses  for  the  modified  Fricke  dosimetry  system, 
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Fig.  38.   Thiocyanate  dosimetry  system  (a)  graphical 

display  of  the  optical  signal  following  pulse 
radiolysis  of  N20-saturated  solution  containing 
10  mM  KCNS;  data  recorded  at  472  nra,  100  nsec 
period,  50  usee  full  scale;  (b)  curve  fit  of 
data  between  x  =  10  to  x  =  40  usee. 
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Fig.  39.   Hydrated  electron  dosimetry  system  (a)  graphical 
display  of  optical  signal  following  pulse  radi- 
olysis  of  N2_saturated  solution  containing  0.1  M 
NaOH  and  0.1  M  ethanol ;  data  recorded  at  63  2  nm, 
20  nsec  period,  10  usee  full  scale;  (b)  curve 
fit  of  data  between  x  =  3  and  x  =  8  usee. 
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Table  4 


Absorbance  Change  and  Absorbed  Dose  for 
Chemical  Dosimetry  Systems 


System  AA  /  pulse  DD,  rads 

Super  Fricke   0.101  ±  0.005  (2.5027  ±  .2741)  x  10 4 

Thiocyanate    a)  0.418  ±  0.076  (8.709  ±  1.739)  x  10 3 

b)  1.129  ±  0.332  (2.3542  ±  .5056)  X  104 

Hydrated       a)  0.378  ±  0.137  (9.220  ±  2.001)  x  10 3 


Electron 


b)  1.103  ±  0.267      (2.6894  ±  .6521)  X  10 4 


a)  Maximum  end  of  pulse  absorption  from  transient  trace 

b)  Results  from  extrapolation  of  best  fit 


APPENDIX  B 
RELEVANT  REACTION  SYSTEMS 

The  reactions  used  for  the  calculations  of  product 
concentrations  in  the  aqueous  tetraphenylborate  and  hexa- 
chloroiridate  systems  are  listed  in  Table  5.   The  100  series 
involves  the  well-known  reactions  describing  transients  of 
irradiated  pure  water152'153  as  well  as  the  reactions  rele- 
vant to  N20-  and  air-saturated  solutions.   Reactions  101  - 
107  and  118  describe  the  recombination  of  the  primary  radi- 
cal  species.   Reactions  108-117  describe  the  reactions 
between  the  radicals  and  the  molecular  products. 

Reactions  118-123  give  the  dissociation  equilibria  of 
H20,  H02,  and  H202.   Reactions  124-128  should  be  important 
at  pH  10  and  above.   Reactions  130-134  describe  the  carbon- 
ate anion  and  the  carbonate  radical  anion,  C03",  produced 
in  the  pulse  radiolysis  of  carbonate  solutions  containing 


N20. 


The  200  series  involves  all  the  reactions  in  the  mecha- 


nistic scheme  describing  the  IrCl63~  complex  under  the 
various  conditions  presented  in  the  text.   The  300  series 
involves  all  the  reactions  in  the  mechanistic  scheme  for  the 
pulse  radiolysis  of  N20-saturated  solutions  containing  the 
TBP~  anion. 
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Table  5 
Chemical  Reaction  Systems  for  Computer  Simulation  Studies 


No. 

Reaction 

Rate  Constant 
(M~1s"1) 

Ref . 

101 

e.,  ~  +  e„  ~ 
aq     aq 

— > 

H2  +  20H~ 

5.4 

X 

109 

(a) 

102 

eaq~  +  OH 

— > 

OH"  +  H20 

3.0 

X 

1010 

(a) 

103 

eaq"  +  H+ 

— > 

H 

2.4 

X 

1010 

(a) 

104 

eaq"  +  H 

— > 

H2  +  OH" 

2.5 

X 

1010 

(a) 

105 

H  +  H 

—  > 

H2 

7.8 

X 

109 

(a) 

106 

OH  +  OH 

— > 

«2°2 

5.3 

X 

109 

(a) 

107 

H  +  OH 

— > 

H20 

7 

X 

109 

(a) 

108 

eaq"  +  H2°2 

— > 

OH  +  OH" 

1.2 

X 

1010 

(a) 

109 

OH  +  H2 

— > 

H  +  H20 

4.9 

X 

107 

(a) 

110 

OH  +  H202 

— > 

H02  +  H20 

2.7 

X 

107 

(a) 

111 

H02  +  H02 

— > 

H202  +  02 

2.7 

X 

106 

(a) 

112 

02"  +  H02 

— > 

H02"  +  02 

4.4 

X 

107 

(a) 

113 

eaq"  +  °2 

— > 

°2~ 

1.9 

X 

1010 

(a) 

114 

H  +  H202 

— > 

OH  +  H20 

1 

X 

1010 

(a) 

115 

H02  +  H 

— > 

H202 

1 

X 

1010 

(a) 

116 

H02  +  OH 

— > 

H20  +  02 

6.6 

X 

109 

(a) 

117 

H  +  02 

— > 

H02 

1.9 

X 

1010 

(a) 

118 

H+  +  OH" 

— > 

H20 

1.4 

X 

1011 

(b) 

119 

H20 

— > 

H+  +  OH" 

2.6 

X 

10"5 

(b) 

120 

02"  +  H+ 

— > 

H02 

4.5 

X 

1010 

(c) 

121 

H02 

— > 

H+  +  02" 

8 

X 

105 

(d) 

130 


122 
123 
124 
125 
126 
127 
128 
129 
130 
131 
132 
133 
134 


H+  +  HO- 


H2°2 

H202  +  OH" 


H02   +  H20 
OH  +  OH~ 


O   +  H20 


02   +  OH 
N20  +  eaq- 


OH  +  CO 


2- 


C03   +  C03 


co3  +  o2 

CO-,2-  +  Ho0 


HCO3   +  OH 


■>  H202 


— >  W     +    HO. 


•>  H02   +  H20 


•>  H202  +  OH 


— >  H20  +  O 
— >  OH  +  OH" 


— >  OH  +  02 
— >  OH"  +  OH 
— >  OH~  +  CO- 


X  10 


10 


— >  CO,  +  CO 


2- 


4 


— >  CO-,2   +  o 


— >  HCO3   +  OH" 
— >  CO-,2-  +  HoO 


3 
1 


X   10 


X   10 


X  10' 


-2 


10 


1.2   X  10 


10 


9.3   X  10' 


1.0   X  10 


8.7   X  10 


4.1   X  10< 


10 


X  10 


4  X  108 
1.78  X  104 
1     X  108 


c 

e] 
c' 
f 

g: 

[h] 
a' 
i 
i 
J 

K 
c' 


201 
202 


IrCl63   +  OH 
IrCl63"  +  ea„~ 


— >  IrCl62"  +  OH" 

202  IrCl63"  +  e   ~  — >  IrClg4" 

203  IrCl63"  +  H     — >  IrCl64~  +  H+ 

204  IrCl62"  +  eaq"  ->  IrClg3" 

205  IrCl62~  +  H     — >  IrCl63~  +  H+ 

206  IrCl64~  +  IrCl62~ 

— >  2IrCl63" 

207  IrCl62~  +  02"   — >  IrCl63~  +  0: 

208  IrCl62" 


2 
+  H02"  — >  IrCl63- 

209 


HO- 


IrCl64"  +  IrCl64" 


IrCl63   +  IrCl65~ 


4  X  109 

6.8  X  109 
6  X  109 

9.9  X  109 
9.2  X  109 

6  X  108 
1  X  107 
4.7   X  108 

5  X  108 


1 
1 

[n] 
n' 

1 

1 
[o] 

n' 
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301  B(C6H5)4"  +  OH-  -->  (C6H5)3BC6H5OH~       6  X  109  (1) 

302  (C6H5)3BC6H5OH"   — >  B(C6H5)3  +  C6H5OH~   5.5  X  104  (1) 

303  C6H5OH~  +  H20  — >  CgHgOH'  +  OH"  1  X  105  (p) 

304  C6H5OH"'  +  C6H5OH"- 

~>  [HO(C6H5)-(C6H5)OH]2~     2  X  108  (c) 

305  [HO(C6H5)-(C6H5)OH]2"  ~>  C6H5-C6H5  +  20H"   .5  X  104  (c) 

306  C6H6OH  +  C6H6OH   ~>  [HO(C6H6 ) - ( C6H5 )0H]   6  X  108  (q) 

307  [HO(C6H6)-(C6H5)OH]  — >  C6H5-C6H5  +  2H20   .3  X  104  (o) 

308  C6H6OH  +  C6H6OH   ~>  C6H6  +  C6H5OH  +  H20   3  X  108  (q) 

309  B(C6H5)4~  +  H-   — >    (C6H5)3BC6H6T       2  X  107  (c) 

310  (C6H5)3BC6H6T  — >  B(C6H5)3  +  CgHgT       1  X  104  (c) 

311  C6H6"  +  H20  — >  C6Hy  +  0H~  1  X  105  (c) 

312  C6Hy-  +  C6H6OH-  ~>  2C6H6  +  H20  3  X  107  (q) 


aref.  153  nref .  111. 

bref.  37.  °to  qive  best  fit. 

cassuraed.  Pref .  141 

dto  give  pK  =  4.8  with  k120-  qrefs.  58,  60. 

eto  give  pK  =  11.8  with  k122. 

to  give  pK  =  2.3  with  k124. 
^ref.  43. 
^ef.  154. 
iref.  138. 
3ref.  155. 
^o  give  pK  =  3.74  with  k134. 

this  work. 
manalogous  to  reaction  205. 
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APPENDIX  C 
INPUT  DATA  FILES 


The  data  files  used  in  the  various  computer  integration 
simulations  for  the  aqueous  TPB  and  Iridium  Chloride  systems 
are  presented  below.   Each  input  file  consist  of  1)  title 
line,  2)  reaction  steps  and  appropriate  rate  constants,  3) 
initial  reactant  molar  concentrations,  4)  integration  step 
size  and  time  period  parameters,  5)  1-50  time  increments. 


Simulation  Tetraphenylborate  System;  N20  saturated 

1110 

1 , EAQ+EAQ=H2+H202 — , 5 . 4E9 

2 , EAQ+OH=OH-+H20 , 3 . 0E10 

3 , EAQ+HPLUS=H , 2 . 4E10 

4 , EAQ+H=H2+OH- , 2 . 5E10 

5,H+H=H2,1.0E10 

6,OH+OH=H202,5.3E9 

7,H+OH=H2O,2.4E10 

8 , HPLUS+OH-=H20 , 3E10 

9 , H20=HPLUS+OH- , 5 . 5E-6 

10 , EAQ+H202=OH+OH- , 1 . 2E10 

11 , OH+H2=H+H20 , 4 . 9E7 

12 , OH+H202=H02 . +H20 , 2 . 7E7 

13 , H02 . +H02 . =H202+02 , 2 . 7E6 

14 , 02-+H02 . =H02-+02 , 4 . 4E7 

15 , EAQ+02=02- , 1 . 9E10 

16 , H+H202=OH+H20 , 1E10 

17 ,H02 . +H=H202 , 1E10 

18 , H02 . +OH=H20+02 , 1E10 

19,H+02=H02. ,1E10 

20,O2-+HPLUS=HO2. ,3E10 

21 ,HPLUS+H02-=H202 , 3E10 

2  2 , H02 . =HPLUS+02- , 1E6 

23 ,H202=HPLUS+H02- , 3E-2 

24,H202+OH-=H02-,1E10 

25,H02-=H202+OH-,3E7 

26 , N20+EAQ=OH-+OH , 8 . 7E9 

27 , OH+OH-=0-+H20 , 1 . 2E10 

28 , 0-+H20=OH+OH- , 9 . 3E7 

29 , BP4-+OH=BP40H- , 6E9 

30 , BP40H-=BP3+POH- , 5 . 5E4 

31 , POH-=PHOH+OH- , 1E5 

32 , PHOH+PHOH=X , 6E8 

33,X=C6H5-C6H5, . 3E4 

34 ,PHOH+PHOH=POH+C6H6 , 3E8 
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35 , P0H-+POH-=Y , 2E8 

36 , Y=C6H5-C6H5 , . 5E4 

37 , BP4-+H=BP4H- , 1E5 

38 , BP4H-=BP3+PH- , 1E4 

3  9 , PH-+H20=PH+0H- , 1 . 0E7 

40 , PHOH+PH=C6H6+C6H6 , 3 . 5E7 

END 

EAQ,1.12E-4 

HPLUS,1.12E-4 

OH,1.12E-4 

H,2.91E-5 

H2, .975E-5 

H202,1.5E-5 

H20,55 

OH-,lE-7 

N20,2.5E-2 

BP4-,lE-2 

END 

1.0E-30,1.0E-33,4E-6, 0.00345 

1,2,3,4,5,6,7,8,9,10, 

11,12,13,14,15,16,17,18,19,20, 

21,22,23,24,25,26,27,28,29,30, 

31,32,3  3,34,35,36,37,38,39,40, 

41,42,43,44,45,46,47,48,49,50 

Simulation  Hexachloroiridate  System;  N20  saturated,  alkaline 

1110 

1 , EAQ+EAQ=H2+H202 — , 5 . 4E9 

2 , EAQ+OH=OH-+H20 , 3 . 0E10 

3 , EAQ+HPLUS=H , 2 . 4E10 

4 , EAQ+H=H2+OH- , 2 . 5E10 

5,H+H=H2,7.8E9 

6 , 0H+OH=H2O2 , 5 . 3E9 

7,H+OH=H20,7E9 

8 , HPLUS+OH-=H20 , 1 . 43E11 

9 , H20=HPLUS+OH- , 2 . 6E-5 

10 , EAQ+H202=OH+OH- , 1 . 2E10 

11 , OH+H2=H+H20 , 4 . 9E7 

12 , OH+H202=H02 . +H20 , 2 . 7E7 

13 , H02 . +H02 . =H202+02 , 2 . 7E6 

14 , 02-+H02 . =H02-+02 , 4 . 4E7 

15 , EAQ+02=02- , 1 . 9E10 

16 , H+H202=OH+H20 , 1E10 

17 , H02 . +H=H202 , 1E10 

18 ,H02 . +OH=H20+02 , 6 . 6E9 

19,H+02=H02. ,1.9E10 

20,O2-+HPLUS=HO2. ,4.5E10 

21 , HPLUS+H02-=H202 , 3E10 

22 , H02 . =HPLUS+02- , 8E5 

23,H202=HPLUS+H02-,3E-2 

24 ,H202+OH-=H02- , 1E10 

25 , H02-=H202+OH- , 5E7 

26 ,OH+OH-=0-+H20 , 1 . 2E10 
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27 , 0-+H20=0H+0H- , 9 . 3E7 

28 , N20+EAQ=0H-+0H , 8 . 7E9 

29 , 0H+C03 — =0H-+C03- , 4 . 1E8 

30 , C03-+C03-=C02+C04 — , 2E7 

31 ,C03 — =HC03-+OH- , 1 . 78E4 

32 , HC03-+OH-=C03 — , 1E8 

33 ,C03-+02-=C03 — +02 , 4E8 

34 ,02-+OH=OH-+02 , 1 . 01E10 

35 , IR3+OH=IR4+OH- , 4E9 

36 , IR3+EAQ=IR2 , 6 . 8E9 

37 , IR2+IR4=IR3+IR3 , 6E8 

38 , IR4+EAQ=IR3 , 9 . 9E9 

39 , IR4+H=IR3+HPLUS , 9 . 2E9 

40 , IR3+H=IR2+HPLUS , 6E9 

41 , IR4+H02-=IR3+H02 . , 4 . 7E8 

42 , IR4+O2-=IR3+02 , 1E7 

43 , IR2+IR2=IR3+IR1 , 5E8 

END 

EAQ,1.5E-4 

HPLUS,1.5E-4 

OH,1.55E-4 

H,3.14E-5 

H2,2.57E-5 

H202,3.88E-5 

H02,1.48E-6 

H02-, .53E-5 

H20,55 

0H-,9E-4 

IR3,lE-3 

N20,2.5E-2 

C03 — ,5E-3 

HC03-,9E-4 

02-,2E-7 

END 

1.0E-30,1.0E-33, 9. 25E-6, 0.00345 

1,2,3,4,5,6,7,8,9,10, 

11,12,13,14,15,16,17,18,19,20, 

21,22,23,24,25,26,27,28,29,30, 

31,32,33,34,35,36,37,38,39,40, 

41,42,43,44,45,46,47,48,49,50 

Simulation  Hexachloroiridate  System;  N20  saturated,  neutral 

1110 

1 , EAQ+EAQ=H2+H202 — , 5 . 4E9 

2 , EAQ+OH=OH-+H20 , 3 . 0E10 

3 , EAQ+HPLUS=H , 2 . 4E10 

4 , EAQ+H=H2+OH- , 2 . 5E10 

5,H+H=H2,7.8E9 

6 , OH+OH=H202 , 5 . 3E9 

7,H+OH=H20,7E9 

8 , HPLUS+OH-=H20 , 1 . 43E11 

9 , H20=HPLUS+OH- , 2 . 6E-5 

10 , EAQ+H202=OH+OH- , 1 . 2E10 


135 


1 1 , OH+H2=H+H20 , 4 . 9E7 

12 ,OH+H202=H02 . +H20 , 2 . 7E7 

13 ,H02 . +H02 . =H202+02 , 2 . 7E6 

14 ,02-+H02 . =H02-+02 , 4 . 4E7 

15 , EAQ+02=02- , 1 . 9E10 

16,H+H2O2=OH+H2O,lE10 

17 , H02 . +H=H202 , 1E10 

18 , H02 . +OH=H20+02 , 6 . 6E9 

19,H+02=H02. ,1.9E10 

20 , 02-+HPLUS=H02 . , 4 . 5E10 

21 ,HPLUS+H02-=H202 , 3E10 

22 , H02 . =HPLUS+02- , 8E5 

23,H202=HPLUS+H02-,3E-2 

24 , N20+EAQ=0H-+0H , 8 . 7E9 

25 ,02-+OH=OH-+02 , 1 . 01E10 

26 , IR3+OH=IR4+OH- , 4E9 

27 , IR3+EAQ=IR2 , 6 . 8E9 

28 , IR2+IR4=IR3+IR3 , 6E8 

29 , IR4+EAQ=IR3 , 9 . 9E9 

30 , IR4+H=IR3+HPLUS , 9 . 2E9 

31 , IR3+H=IR2+HPLUS , 6E9 

32 , IR4+H02-=IR3+H02 . , 4 . 7E8 

33 , IR4+02-=IR3+02 , 1E7 

34 , IR2+IR2=IR3+IR1 , 5E8 

END 

EAQ,1.5E-4 

HPLUS,1.5E-4 

OH,1.55E-4 

H,3.14E-5 

H2,2.57E-5 

H202,3.88E-5 

H02,1.48E-6 

H20,55 

0H-,lE-7 

IR3,lE-3 

N20,2.5E-2 

END 

1.0E-30,1.0E-33, 9. 25E-6, 0.00345 

1,2,3,4,5,6,7,8,9,10, 

11,12,13,14,15,16,17,18,19,20, 

21,22,23,24,25,26,27,28,29,30, 

31,32,33,34,35,36,37,38,39,40, 

41,42,43,44,45,46,47,48,49,50 

Simulation  Hexachloroiridate  System;  Nitrogen  saturated 

1110 

1 , EAQ+EAQ=H2+H202 — , 5 . 4E9 

2 , EAQ+OH=OH-+H20 , 3 . 0E10 

3 , EAQ+HPLUS=H , 2 . 4E10 

4 , EAQ+H=H2+OH- , 2 . 5E10 

5,H+H=H2,7.8E9 

6 , OH+OH=H202 , 5 . 3E9 

7,H+OH=H20,7E9 
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8 , HPLUS+0H-=H20 , 1 . 43E11 

9 , H20=HPLUS+0H- , 2 . 6E-5 

10 , EAQ+H202=OH+OH- , 1 . 2E10 

11 , OH+H2=H+H20 , 4 . 9E7 

12 ,OH+H202=H02 . +H20 , 2 . 7E7 

13 , H02 . +H02 . =H202+02 , 2 . 7E6 

14 ,02-+H02 . =H02-+02 , 4 . 4E7 

15 , EAQ+02=02- , 1 . 9E10 

16 ,H+H202=OH+H20 , 1E10 

17 , H02 . +H=H202 , 1E10 

18 , H02 . +OH=H20+02 , 6 . 6E9 

19,H+02=H02. ,1.9E10 

20 ,02-+HPLUS=H02 . , 4 . 5E10 

21 ,HPLUS+H02-=H202 , 3E10 

22 , H02 . =HPLUS+02- , 8E5 

23,H2O2=HPLUS+H02-,3E-2 

24 ,02-+OH=OH-+02 , 1 . 01E10 

25 , IR3+OH=IR4+OH- , 4E9 

26 , IR3+EAQ=IR2 , 6 . 8E9 

27 , IR2+IR4=IR3+IR3 , 6E8 

28 , IR4+EAQ=IR3 , 9 . 9E9 

29 , IR4+H=IR3+HPLUS , 9 . 2E9 

30 , IR3+H=IR2+HPLUS , 6E9 

31 , IR4+H02-=IR3+H02 . , 4 . 7E8 

32 , IR4+02-=IR3+02 , 1E7 

33  , IR2+IR2=IR3+IR1 , 5E8 

END 

EAQ,1.5E-4 

HPLUS,1.5E-4 

OH,1.55E-4 

H,3.14E-5 

H2,2.57E-5 

H202,3.88E-5 

H02,1.48E-6 

H20,55 

0H-,lE-7 

IR3,lE-3 

END 

1.0E-30,1.0E-3  3,9.25E-6,0.00345 

1,2,3,4,5,6,7,8,9,10, 

11,12,13,14,15,16,17,18,19,20, 

21,22,23,24,25,26,27,28,29,30, 

31,32,33,34,35,36,37,38,39,40, 

41,42,43,44,45,46,47,48,49,50 

Simulation  Hexachloroiridate  System;  Aerated 

1110 

1 , EAQ+EAQ=H2+H202 — , 5 . 4E9 

2 , EAQ+OH=OH-+H20 , 3 . 0E10 

3 , EAQ+HPLUS=H , 2 . 4E10 

4 , EAQ+H=H2+OH- , 2 . 5E10 

5,H+H=H2,7.8E9 

6 , OH+OH=H202 , 5 . 3  E9 
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7,H+OH=H20,7E9 
8 ,HPLUS+0H-=H20 , 1 . 43E11 
9 , H20=HPLUS+0H- , 2 . 6E-5 
10 , EAQ+H202=OH+OH- , 1 . 2E10 
, 0H+H2=H+H2O , 4 . 9E7 
, OH+H202=H02 . +H20 , 2 . 7E7 
, H02 . +H02 . =H202+02 , 2 . 7E6 
, 02-+H02 . =H02-+02 , 4 . 4E7 
, EAQ+02=02- , 1 . 9E10 
, H+H202=OH+H20 , 1E10 
,HO2.+H=H2O2,lE10 
, H02 . +OH=H20+02 , 6 . 6E9 
,H+02=H02. ,1.9E10 
,02-+HPLUS=H02. ,4.5E10 
, HPLUS+H02-=H202 , 3E10 
, H02 . =HPLUS+02- , 8E5 
, H202=HPLUS+H02- , 3E-2 
,02-+OH=OH-+02 , 1 . 01E10 
, IR3+OH=IR4+OH- , 4E9 
, IR3+EAQ=IR2 , 6 . 8E9 
, IR2+IR4=IR3+IR3 , 6E8 
,IR4+EAQ=IR3,9.9E9 
, IR4+H=IR3+HPLUS , 9 . 2E9 
, IR3+H=IR2+HPLUS , 6E9 
, IR4+H02-=IR3+H02 . , 4 . 7E8 
, IR4+02-=IR3+02 , 1E7 
, IR2+IR2=IR3+IR1 , 5E8 


11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

END 

EAQ,1.5E-4 

HPLUS,1.5E-4 

0H,1.55E-4 

H,3.14E-5 

H2,2.57E-5 

H202,3.88E-5 

H02,1.48E-6 

H20,55 

0H-,lE-7 

IR3,lE-3 

02,2.5E-4 

END 

1. 0E-30,1.0E-33, 9. 4E-6, 0.00345 

1,2,3,4,5,6,7,8,9,10, 

11,12,13,14,15,16,17,18,19,20, 

21,22,23,24,25,26,27,28,29,30, 

31,32,33,34,35,36,37,38,39,40, 

41,42,43,44,45,46,47,48,49,50 


APPENDIX  D 
PULSE  RADIOLYSIS  COMPUTER  PROGRAM 


1 
2 
3 
4 
5 
6 
7 

8 

9 

10 

11 

12 

15 

18 
19 
20 
21 
23 

24 
25 
26 

28 

33 

35 

40 

43 

45 

48 
50 
52 

72 

73 


**********    DATA  GRABBER    **************** 

*****  File  DEBEST8.BAS,  OCT.  12,  1989  ***** 

PULSE  RADIOLYSIS  DATA  ACQUISITION  PROGRAM  FOR  DATA 

PRECISION  6000 

TRANSIENT  DIGITIZER  INTERFACED  THROUGH  A  SERIAL  PORT  AT 

9600  BAUD 

WITH  AN  IBM-PC  COMPATIBLE  MICROCOMPUTER. 

FOR  USE  WITH  MICROSOFT  QUICK  BASIC  COMPILER  VERSION  4.0 

WITH  8087 

COPYRIGHT  1988,  1989  BY  JOHN  BROGDON  AND  ROBERT 

HANRAHAN 

RADIATION  CHEMISTRY  LAB,  406  NSC,  U  OF  FLA  GAINESVILLE 

32611 

TEL  904-392-1442  OR  376-7754 

WORK  DONE  UNDER  DOE  CONTRACT  DE-ASO5-76ERO-3106 

COMPILE  FROM  THE  QBASIC  DIRECTORY  WITH  THE  WORKFILE  ON 


B: 


AS  FOLLOWS 


BC  B:  DEBEST8/O/X/E/C:6000;  THE  LAST  SETS  UP  THE 

COMMUNICATION  BUFFER;  THEN  LINK.   IN  TURBO  BASIC  THE 

METACOMMAND  $C0M1  6000  IS  SUPPOSED  TO  WORK  BUT  WE  GET 

ERRORS 

ABOUT  EVERY  NINTH  DATA  POINT  WHEN  WE  TAKE  DATASET  FROM 

DP  6000. 

IN  BASICA,  THE  BUFFER  MAY  BE  CHANGED  WITH  THE  COMMAND 

"BASIC/C:6000" 

WHEN  ENTERING  INTO  THE  BASIC  MODE,  BUT  INTERPRETER  IS 

VERY  SLOW 

AND  SEVERAL  TIMING  LOOPS  MUST  BE  CHANGED.  SOME  TIMING 

LOOPS 

PRESENTLY  MAY  BE  LONGER  THAN  NECESSARY  ON  A  STANDARD 

SPEED  4.77  MHz 

COMPUTER  ****  BUT  MIGHT  NEED  TO  BE  LONGER  ON  A  FASTER 

MACHINE 


****  NOTE,  THE  UNDERLINE  IS  A  LINE  CONTINUATION 

IN  QBASIC  **** 
WARNING  —  NOT  ALL  FEATURES  HAVE  BEEN  FULLY  TESTED.  IN 
PARTICULAR, 

THE  COMMAND  TO  COMPARE  SPECTRA   HAS  NOT  BEEN  USED 
EXTENSIVELY. 


138 


139 


75 
80 

90 

100 
110 
120 
130 
140 
145 
150 
170 
180 
190 
200 
210 
220 
230 
240 
250 
260 
270 
280 
290 
300 
310 
320 

330 
340 
350 
360 
370 
380 
390 
400 
410 
420 
430 
432 
433 
435 
436 
437 
438 
440 
450 
455 
460 
470 
480 
490 
500 


DIM  AGI$(100) ,  BGI$(100),  DSPM$(10),  TRGSRC$(10),  OD(512) 
DIMQY(520),  QX(512),  VERLN(147),  QLY(512),  BTARY(3400), 

CON(512) 
OPEN  "coml:9600,n,8,l,cs,ds,cd"  FOR  RANDOM  AS  #1 
ON  ERROR  GOTO  5220 

SCREEN  2:  KEY  ON:  CLS  :  GET  (1,  80)- (639,  99),  BTARY 
DLY$  =  "-3.0US":  PER$  =  "200nS":  SWP$  =  "3" 
DSPM  =  2:  YSC$  =  "1.0":  TRGSRC  =  5 

PRINT  "FIRST  TURN  ON  THE  DP6000,  THEN  PRESS  CONT" 
FOR  KKQ  =  1  TO  520:  QY(KKQ)  =  0:  NEXT 
L  =  1:  EE  =  3000:  FLORD  =  0:  FLCAL  =  0:  FLDATA  =  0 
KEY  OFF:  B$  =  "PROG;KEYSRQ=1;KEY=1017;KEY=1033;REMOTE" 
GOSUB  1090 

INPUT  "IS  MODE  ON  RUN  OR  EDIT";  B$ 
B$  =  LEFT$(B$,  1) 

IF  B$  <>  "R"  AND  B$  <>  "E"  GOTO  190 
IF  B$  =  "R"  GOTO  250 
B$  =  "KEY=1017" 
GOSUB  1090 

INPUT  "IS  EXECUTE  ON  RUN/STOP  (YES  or  NO)";  B$ 
B$  =  LEFT$(B$,  1) 

IF  B$  <>  "Y"  AND  B$  <>  "N"  GOTO  250 
IF  B$  =  "Y"  GOTO  320 
B$  =  "KEY=1033" 
GOSUB  1090 
GOTO  250 

DEF  FNQSC  (QS,  QXS ,  QNS ,  QXO ,  QNO)  = 
/  (QXS  -  QNS)  *  (QXS  -  QS) 
PRINT  :  GOTO  350 
'***************  MENU  ************************************ 


INT (QXO  -  (QXO  -  QNO) 
+  .5) 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
FOR  I 
SOUND  RND 


:  PRINT  "RETURNING 

"CALIBRATION 

"DATA 

"PLOTTER 

"LEAST  SQUARES 

"LOAD 

"SAVE 

"ORDER 

"COMPARE 

"MAN 

"DIR 

"END 

=  1  TO  25 

*  1000  +  37 


FOR  EXPERIMENT" 


TO  MENU":  PRINT 
CALIBRATE  DP6000 
AQUIRE  DATA  FROM  DP6000" 
PLOT  RAW  DATA  ON  SCREEN" 
LINEAR  REGRESION  FOR  CURRENT  DATA" 
LOAD  DATA  FROM  DISK" 
SAVE  DATA  ON  DISK" 
CHOOSE  ORDER  OF  REACTION" 
COMPARE  DATA  SETS  ON  GRAPH" 
MANUAL  CHANGE  V0,VM,L,OR  E" 
DISK  DIRECTORY" 
END  SESSION  ■ 


NEXT  I 

PRINT  :  INPUT  "COMMAND";  B$ 

B$  =  LEFT$(B$,  3) 

IF  B$  =  "END"  THEN  END 

RESTORE 

FOR  I  =  1  TO  11 

IF  I  =  11  THEN  PRINT  "BAD  COMMAND": 

READ  A$ 

IF  B$  =  A$  THEN  GOTO  520 


PRINT 


GOTO  440 
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510   NEXT  I 

520   ON  I  GOSUB  570,  1140,  1480,  3950,  4940,  4780,  5140,  5330, 

8000,  9800 
530   PRINT  :  GOTO  350 
550   DATA  "CAL" , "DAT" , "PLO" , "LEA" , "LOA" , "SAV" , "ORD" , "COM" , 

"MAN","DIR" 
555   GOTO  9999 

560  '************  CALIBRATION  ******************************** 
570   PRINT  :  PRINT  TAB(25);  :  PRINT  "CALIBRATION" 
580   DSPM$(1)  =  "SINGLE":  DSPM$(2)  =  "2  SEPR" :  DSPM$(3)  = 

"2  OVLY" 

585  TRGSRC$(1)  =  "NONE":  TRGSRC$(2)  =  "CH  1":  TRGSRC$(3)  = 

"CH  2" 

586  TRGSRC$(4)  =  "LINE":  TRGSRC$(5)  =  "EXT  TRIG" 
590   PRINT 

592  PRINT  "1)TIME  DELAY="  +  DLY$ ,  "2)PERIOD="  +  PER$ ,  "3)# 

SWEEPS="  +  SWP$, 

593  PRINT  "4)DSPL  MODE="  +  DSPM$(DSPM),  "5)   Y  SCALE="  +  YSC$ , 

"6) NO  CHANGE" 
595   PRINT  "7)TRIGGER="  +  TRGSRC$ (TRGSRC) 
600   FLCAL  =  1 
610   PRINT 

620   PRINT  "WHICH  NUMBER":  B$  -  INPUT$(1) 
630   B  =  VAL(B$) 

640   IF  B  <  1  OR  B  >  7  GOTO  620 

650   ON  B  GOTO  660,  670,  680,  690,  700,  710,  702 
660   INPUT  "1)TIME  DELAY=  ",  DLY$:  GOTO  590 
670   INPUT  "2)PERIOD=",  PER$ :  GOTO  590 
680   INPUT  "3)#  SWEEPS=" ,  SWP$ :  GOTO  590 
690   PRINT  :  PRINT  "CHOICES:",  "1)  SINGLE",  "2)  2  SEPR", 

"3)  2  OVLY" 
695   PRINT  :  INPUT  "4)DSPL  MODE=" ,  DSPM:  GOTO  590 
700   INPUT  "5)Y  SCALE=" ,  YSC$ :  GOTO  590 
702   PRINT  :  PRINT  "CHOICES:",  "l)NONE",  "2)CH  1",  "3)CH  2", 

"4)LINE", 
705   PRINT  "5) EXT  TRIG":  INPUT  "7 )TRIGGER=" ,  TRGSR:  GOTO  590 
710   PRINT  "INIALIZING  DP6000" 
715   TRG$  =  STR$( TRGSR) 
720   B$  =  "DARM;PROMPT=l;TRIG;TRGM=l;TMB;PERIOD="  +  PER$  + 

";TRGSRC="  +  TRG$ 
730   GOSUB  1090 
732   B$  =  "TMB;KEY=1008;PERIOD="  +  PER$  +  " ;KEY=1004 ; 

PERIOD="  +  PER$  +  "" 
734   GOSUB  1090 

740   B$  =  "DELAY="  +  DLY$  +  " ;DISP;DSPM="  +  STR$(DSPM) 
750   GOSUB  1090 
760   B$  =  "TRACE=1;X;TRCSRC=BUF.A2;XFILL=2;Y;YSCL="  +  YSC$  +" 

;TRACE=2" 
770   GOSUB  1090 

780   B$  =  "TRCSRC=BUF.A2;YSCL="  +YSC$  +  " ;X;XFILL=2 ;GRID=2 ; ARM" 
790   GOSUB  1090 
792   B$  -  "TMB" 
794   GOSUB  1090 
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PRINT  :  PRINT  "V0=";  VO ,  "VM=" ;  VM  +  VO :  PRINT 
INPUT  "DO  YOU  WANT  TO  CHANGE  THE  PARAMETERS  OF  THE 

EXPERIMENT";  B$ 
B$  =  LEFT$(B$,  1) 

IF  B$  <>  "Y"  AND  B$  <>  "N"  GOTO  810 
B$  =  "KEY=4001" 
GOSUB  1090 
B$  =  "KEY=1008" 
GOSUB  1090 
IF  B$  =  "N"  AND  VM 
990 

IF  VM  =  0  AND  B$  = 
V0=0)" 
PRINT 
B$ 
B$ 
NO 


800 
810 

820 
830 
832 
834 
836 
838 
840 

850 

860 
870 
880 
882 
884 
886 
888 
910 
912 
914 
916 
918 
920 
930 
932 
934 
936 
938 
939 
940 
950 
960 
970 
980 
990 
1000 


1010  PRINT 


<>  0  THEN  B$  =  "DARM":  GOSUB  1090:  GOTO 
"N"  THEN  PRINT  "YOU  HAVE  NO  CHOICE  (VM- 


INPUT  "HAVE  YOU  TRIGGERED  FOR  NO  LIGHT";  B$ 


IF  B$  <>  "Y" 
GOSUB  7000 
"BUF.A2(150)H 
"BUF.A2(250)M 
"BUF.A2(350)" 


THEN  860 


GOSUB 
GOSUB 
GOSUB 


7000 
7000 
7000 


LEFT$(B$,  1) 

"BUF.A2(2)": 

MYINFO:  B$  = 
Yl  =  MYINFO:  B$  = 
Y2  =  MYINFO:  B$  = 
Y3  =  MYINFO 

VO  =  (Yl  +  Y2  +  Y3)  /  3 
B$  =  "KEY=4001" 
GOSUB  1090 
B$  =  "KEY=1008" 
GOSUB  1090 
PRINT  :  INPUT  "HAVE  YOU  TRIGGERED  FOR  FULL  LIGHT"; 


B$ 


1020 
1030 
1040 
1050 
1060 
1070 
1080 
1090 
1100 
1105 
1110 


B$  =  LEFT$(B$,  1):  IF  B$  <>  "Y"  THEN  920 
B$  =  "BUF.A2(2)":  GOSUB  7000 

NO  =  MYINFO:  B$  =  "BUF. A2 ( 150 ) " :  GOSUB  7000 
Yl  =  MYINFO:  B$  =  "BUF.A2(250)M:  GOSUB  7000 
Y2  =  MYINFO:  B$  =  "BUF. A2 ( 350) " :  GOSUB  7000 
Y3  =  MYINFO 

REM  B$="BUF.A2(2) ;BUF.A2(400) ;BUF.A2(300) ;BUF.A2(500) ;DARM" 
REM  GOSUB  1090 
REM  INPUT  #1,N0,Y1,Y2,Y3 
VM  =  (Yl  +  Y2  +  Y3)  /  3 
VFULL  =  VM  -  VO 
REM  L  =  L  /  EE 
PRINT  :  PRINT  "VMAX 
»E="; 
"L-CHANGE  L", 
INPUTS ( 1 ) 
IF  B$  <>  "L"  AND  B$ 
IF  B$  =  "L"  THEN  INPUT 
IF  B$  -  "E"  THEN  INPUT 
LEFF  =  L  *  EE  *  2.303 


=  ";  VM;  ' 
EE:  PRINT 
"E-CHANGE 


VO  = 


VO;  "L=! 


<> 


"E"  AND 
"L="  , 
"E="  , 


E",  "N-NO  CHANGE   " : 

B$  <>  "N"  GOTO  1010 
L:  GOTO  1000 
EE:  GOTO  1000 


B$  = 


QY(516)  = 

RETURN 

' ******** 

PRINT  #1, 
FOR  1=1 
QQWER  =  LOG(I) 
NEXT  I 


VO:  QY(517)  =  VM:  QY(518)  =  L:  QY(519)  =  EE 

OUTPUT  B$  TO  DP6000  *************************** 

B$ 

TO  8000 
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1120  RETURN 

1130  '***********  DATA  *************************************** 

1140  PRINT  TAB(40);  "DATA" 

1150  IF  FLCAL  =  0  THEN  PRINT  "CALIBRATION  NEEDED  FIRST":  GOSUB 

570 
1160  B$  =  "DARM;X;TRACE=1;TRCSRC=BUF.A2" 
1170  GOSUB  1090 

1172  B$  =  "X;TRACE=2;TRCSRC=AVEGA2" 

1173  GOSUB  1090 

1174  FOR  BADDOG  =  1  TO  500:  AHUNK1  =  LOG(BADDOG):  NEXT 

1175  REM  B$  =  "PROG;SAVG;AVGCLR;AVGCNT=0;NAVG="  +  SWP$  + 
";KEY=1067;DISP;" 

1176  B$  =  "PROG;SAVG;AVGCLR;AVGCNT=0;" 

1177  GOSUB  1090 

1178  B$  =  "NAVG="  +  SWP$  +  " ;KEY=1067 ;DISP; " 

1179  GOSUB  1090 

1182  B$  =  "AVGCLR;AVGCNT=0;CLRSUM" 

1184  GOSUB  1090 

1192  B$  =  "X;TRACE=l;TRCSRC=BUF.A2;ARM" 

1194  GOSUB  1090 

1200  SWP  =  VAL(SWP$) 

1210  FOR  K  =  1  TO  SWP 

1212  B$  =  "KEY=4001" 

1214  GOSUB  1090 

1216  B$  =  "KEY=1008" 

1218  GOSUB  1090 

1220  PRINT  "TYPE  T  FOR  TRIGGERED  OR  S  FOR  STOP.   "; 

1230  B$  =  INPUT$(1):  PRINT  B$ 

1240  IF  B$  <>  "T"  AND  B$  <>  "S"  THEN  GOTO  1220 

1250  IF  B$  =  "S"  GOTO  1300 

1260  B$  =  "KEY=2001" 

1270  GOSUB  1090 

1280  PRINT  "#  SWEEPS=" ;  K 

1290  NEXT  K 

1300  B$  =  "DARM;X;TRACE=2;TRCSRC=AVEGA2" 

1309  GOSUB  1090 

1310  REM  GOSUB  7000 

1311  B$  =  "PERIOD":  GOSUB  7000 

1312  NO  =  MYINFO 

1313  B$  =  "PERIOD":  GOSUB  7000 

1314  PER  =  MYINFO 

1315  B$  =  "MIN":  GOSUB  7000 

1316  QMINY  =  MYINFO 

1317  B$  =  "MAX":  GOSUB  7000 

1318  QMAXY  =  MYINFO 

1319  B$  =  "DELAY":  GOSUB  7000 

1320  N  =  MYINFO 

1329  REM   INPUT  #1 , NO , PER, QMINY , QMAXY ,N 

1330  IF  QMAXY  <>  QMINY  THEN  GOTO  1370 

1340  PRINT  "BAD  DATA":  PRINT  NO,  PER,  QMINY,  QMAXY,  N 
1350  REM  IF  LOC(l)<3  GOTO  1300 
1360  REM  INPUT  #l,PER:GOTO  1300 
1365  GOTO  350 
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1370  PRINT  "creating  file  from  avega2" 

1371  PRINT  #1,  "AVEGA2!" 

1372  FOR  AAA  =  1  TO  200:  SAM  =  SAM  +  3:  NEXT 
1378  PRINT  "finist  sam  stall" 

1380  FOR  I  -    1  TO  512 
1383  REM  QWQWQE  =  2+3 
1390  INPUT  #1,  QY(I) 
1400  NEXT  I 

1405  PRINT  "FINISH  W  INPUT  INTO  #1" 
1410  FOR  I  =  1  TO  512 
1420  QX(I)  =  I  *  PER 

1425  IF  VM  =  0  THEN  VM  =  -3:  IF  L  =  0  THEN  L  =  39:  IF  V0  =  0 
THEN  V0  =  -1 

1430  CON(I)  =  -LOG(ABS((QY(I)  -  VO )  /  VFULL) )  /  LEFF 

1431  OD(I)  =  -LOG(ABS( (QY(I)  -  VO )  /  VFULL))  /  2.303 
1440  NEXT  I 

1450  QY(513)  =  PER:  QY(514)  =  QMAXY:  QY(515)  =  QMINY 

1460  FLDATA  =  1 

1465  PRINT  "going  to  main  menu" 

1470  RETURN 

1480  '****************  PLOTTER  ****************************** 

1490  N  =  1:  QMINX  =  QX(1):  QMAXX  =  QX(512) 

1491  RHGRFLAG  =  0 

1500  IF  FLDATA  >  0  THEN  GOTO  1520 

1510  PRINT  "NO  DATA  SOURCE.  RETURNING  TO  MENU":  RETURN 

1520  INPUT  "WHAT  RESOLUTION  WOULD  YOU  LIKE  TO  GRAPH  (1  TO  20)"; 

QRES 
1525  IF  QRES  <  1  OR  QRES  >  20  THEN  QRES  =  1 
1530  QN  =  512:  QFLG  =  0:  GOSUB  2160 
1540  GOSUB  1570 
1550  RETURN 

1560  '****************  LEFT  &  RIGHT  BORDERS  ****************** 

1561  RHGRFLAG  =  0 

1570  PUT  (0,  QLNB  +11),  BTARY ,  PSET:  LOCATE  INT( (QLNB  /  200)  * 
25  +  1)  +  2,  1 

1571  IF  RHGRFLAG  =  1  THEN  FOR  ASEC  =  1  TO  12000:  AMIN  = 
LOG(ASEC):  NEXT 

1575  PRINT  "WOULD  YOU  LIKE  TO  LEAST  SQUARES  FIT  DATA  ON  GRAPH" 

1590  INPUT  B$:  B$  =  LEFT$(B$,  1) 

1600  IF  B$  =  "N"  THEN  RETURN 

1610  IF  B$  <>  "Y"  THEN  GOTO  1590 

1620  GET  (QLNL,  QLNT)-(QLNL,  QLNB),  VERLN 

1630  PUT  (0,  QLNB  +11),  BTARY,  PSET:  LOCATE  INT((QLNB  /  200)  * 

25  +  1)  +  2,  1 
1640  INPUT  "DO  YOU  WISH  TO  MOVE  THE  LEFT  BORDER";  B$ 
1650  B$  =  LEFT$(B$,  1) 
1660  LNPOS  =  QLNL 
1670  IF  B$  =  "N"  THEN  GOTO  1830 
1680  IF  B$  <>  "Y"  THEN  GOTO  1630 
1690  LNPOS  =  QLNL  +  25 
1700  C$  =  "R" 
1710  PUT  (LNPOS,  QLNT),  VERLN 
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1720 

1730 

1740 
1750 
1760 
1770 
1780 
1790 
1800 
1810 
1820 
1830 

1840 
1850 

1860 
1870 
1880 
1890 
1900 
1910 
1920 
1930 
1940 

1950 

1960 
1970 
1980 
1990 
2000 
2010 
2020 
2030 
2040 
2060 
2070 
2075 
2076 
2077 
2078 
2079 
2080 
2081 
2082 

2083 
2084 
2093 


BTARY,  PSET:  LOCATE  INT((QLNB  /  200)  ■ 

;  :  B$ 


PUT  (0,  QLNB  +  11) 

25  +  1)  +  2,  1 

PRINT  "L-LEFT,R-RIGHT,S-STOP  OR  NUMBER  OF  SPACES 

=  INPUT$(1) 

D  =  VAL(B$)  A  2 

IF  D  =  0  AND  (B$  =  ML"  OR  B$  =  "R"  OR  B$  =  "S")THEN  C$=  B$ 

IF  C$  =  "S"  THEN  GOTO  1830 

PUT  ( LNPOS ,  QLNT ) ,  VERLN 

IF  C$  =  "R"  THEN  LNPOS  =  LNPOS  +  D 

IF  C$  =  "L"  THEN  LNPOS  =  LNPOS  -  D 

IF  LNPOS  <  QLNL  THEN  LNPOS  =  QLNL: 

IF  LNPOS  >  QLNR  THEN  LNPOS  =  QLNR: 

PUT  (LNPOS,  QLNT),  VERLN:  GOTO  1720 

N  -  ((QMAXPX  -  QMINPX)  *  QSCLX  *  (LNPOS  -  QLNL))  /  (QLNR  - 

QLNL) 

N  =  INT((N  - 

PUT  (0,  QLNB 

25  +  1)  +  2, 

INPUT  "DO  YOU 


C$ 
C$ 


"R" 

"L" 


QX(i)) 

+  11), 

1 

WISH  TO 


/  QY(513)) 

BTARY,  PSET:  LOCATE  INT ( (QLNB  /  200)  * 

MOVE  THE  RIGHT  BORDER";  B$ 


LOCATE  INT((QLNB  /  200)  * 

;  :  B$ 


B$  =  LEFT$(B$,  1) 

LNPOS  =  QLNR 

IF  B$  =  "N"  THEN  GOTO  2060 

IF  B$  <>  "Y"  THEN  GOTO  1850 

LNPOS  =  QLNR  -  100 

C$  =  "L" 

PUT  ( LNPOS ,  QLNT ) ,  VERLN 

PUT  (0,  QLNB  +  11),  BTARY,  PSET: 

25  +  1)  +  2,  1 

PRINT  "L-LEFT,R-RIGHT,S-STOP  OR  NUMBER  OF  SPACES 

=  INPUT$(1) 

D  =  VAL(B$)  A  2 

IF  D  =  0  AND  (B$  =  "L"  OR  B$  =  "R"  OR  B$  =  "S")THEN  C$=  B$ 

IF  C$  =  "S"  THEN  GOTO  2060 

PUT  ( LNPOS ,  QLNT ) ,  VERLN 

IF  C$  =  "L"  THEN  LNPOS  =  LNPOS  -  D 

IF  C$  =  "R"  THEN  LNPOS  =  LNPOS  +  D 

IF  LNPOS  <  QLNL  THEN  LNPOS  =  QLNL:  C$ 

IF  LNPOS  >  QLNR  THEN  LNPOS  =  QLNR:  C$ 

PUT  (LNPOS,  QLNT),  VERLN:  GOTO  1940 

J  =  ((QMAXPX  -  QMINPX) *QSCLX  *( LNPOS  -  QLNL)) /(QLNR  -  QLNL) 

J  =  INT((J  -  QX(1))  /  QY(513)) 

IF  J  <  N  THEN  CC  =  J:  J=N:  N=CC 

IF  J  >  512  THEN  J  =  512 

IF  N  <  1  THEN  N  =  1 

CLOSE  #1 

INPUT  "CREATE  A  NEW  FILE  FOR  COMP  SIMU 

C$  =  (LEFT$(RNBANS$,  1)) 

IF  C$  =  "N"  OR  C$  =  ""  THEN  GOTO  2106 

INPUT  "CREAT  'OD'  FILE  (YES)  ELSE  CON  FILE  CREATED?" 

CLCANS$ 

C$  =  ( LEFTS (CLCANS$,  1)) 

IF  C$  =  "Y"  THEN  GOTO  9892 

INPUT  "WHICH  DRIVE";  DRV$ :  IF  DRV$  = 


"R" 
»L" 


■?•! 


RNBANS$ 


THEN  DRV$  =  "C" 
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2094 
2095 
2096 
2097 
2098 
2099 
2100 
2101 
2102 
2103 
2104 
2105 
2106 

2110 
2120 
2130 

2140 
2150 
2160 
2170 
2180 

2190 
2200 
2210 

2220 
2230 
2240 
2250 
2260 
2270 
2280 
2290 
2300 
2310 
2320 
2330 
2340 
2350 
2360 
2370 
2380 
2390 
2400 
2410 
2420 
2430 
2440 
2450 
2460 


DRV$  -  LEFT$(DRV$,  1) 

IF  DRV$  =  "A"  OR  DRV$  ■  "B" 

IF  DRV$  =  "C"  OR  DRV$  =  "D" 

GOTO  2093 

DRV$  =  DRV$  +  ":" 

INPUT  "FILE  NAME";  FILENAME$ 

B$  =  DRV$  +  LEFT$(FILENAME$, 

OPEN  "O",  #1,  B$ 

TO  J  STEP  5 

QX(I),  »,»,  CON(I), 


FOR  I  =  N 
PRINT  #1, 
NEXT  I 
CLOSE  #1 
PUT  (0,  QLNB 
25  +  1)  +  2, 


GOTO 
GOTO 


8) 


ii  ii 


2098 
2098 


11),  BTARY,  PSET:  LOCATE  INT((QLNB  /  200)  * 


ORDER 


HERE  COMES  PLOTTER********************** 


=  1  -  grid  only,  =  2  -  plot 


INPUT  "ORDER  OF  REACTION  (0,1,2  or  3 ) " ; 

IF  ORDER  <  0  OR  ORDER  >  3  GOTO  2100 

PUT  (0,  QLNB  +  11),  BTARY,  PSET:  LOCATE  INT ( (QLNB  /  200) 

25  +  1)  +  2,  1 

GOSUB  3980 

RETURN 

QQYY  =  95:  QQXX  -  317 
REM  qflg=0  -  grid  and  plot 
only 

IF  QFLG  >  1  THEN  2910 
Q$  =  "##.#" 

REM  approximate  number  of  vertical  and  horizontal 
intervals 

QNINX  =  5:  QNINY  =4.8 
REM  left  and  right  borders 
QLNL  =  40:  QLNR  =  639 
REM  bottom  and  top  borders 


164:  QLNT  =  9:  QLOCX1 
2:  CLS  :  KEY  OFF 


-  1 


INT (QLNT 
QLOCY1  + 
INT (QLNB 
QLOCY3  + 
QLOCY4  + 
INT (QLNL 
QLOCX2  + 
QLOCX3  + 


8  + 


8  + 


.5) 
.625) 


/ 

1 

/ 

1 
1 

/ 
1 
INT((QLNR  -  QLNL)  /  8  +  .625)  -  1 


8-4) 


QLNB  = 

SCREEN 

QLOCY1 

QLOCY2 

QLOCY3 

QLOCY4 

QLOCY5 

QLOCX2 

QLOCX3 

QLOCX4 

REM  find 

QMIN  =  QMINX: 

GOSUB  3500 

QSCLX  =  QSCL:  QMINPX  =  QMINP: 

QNINPX  =  QNINP:  QWDTX  =  QWDT 

REM  find  "nice"  y-scale  interals 

QMIN  =  QMINY:  QMAX  =  QMAXY:  QNIN  =  QNINY 

GOSUB  3  500 

QSCLY  =  QSCL:  QMINPY  =  QMINP:  QMAXPY  =  QMAXP 

QNINPY  =  QNINP:  QWDTY  =  QWDT 

REM  draw  title 


'nice"  x-scale  intervals 

QMAX  =  QMAXX:  QNIN  =  QNINX 

QMAXPX  =  QMAXP 
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2470  LOCATE  1,  9 

2480  IF  FLDATA  =  1  THEN  DSRC$  =  "DP6000" 
2490  IF  FLDATA  =  2  THEN  DSRC$  =  FILENAME$ 
2492  RH33  =  QSCLX:  RH33  =  RH33  *  (1.0001) 

2500  REM  PRINT  "MULTIPLY  Y-AXIS  BY" ;QSCLY ; "Units 

2501  PRINT  "MULT.  X-AXIS  BY  "; 

2502  PRINT  USING  "#.##AAAA";  RH33;  :  PRINT  "  Sec.   Y-AXIS  BY  "; 

2503  PRINT  USING  "# . ##A A A A"  ;  QSCLY ;  :  PRINT  "  Units.  SOURCE:"; 
DSRC$ 

2510  IF  FLORD  =  0  GOTO  2620 

2520  LOCATE  2,  35 

2530  IF  ORDER  =  1  THEN  B$  =  "ln(C)=ln(Co) -K*t" 

2540  IF  ORDER  =  2  THEN  B$  =  " ( 1/C)=( l/Co)+K*t" 

2550  IF  ORDER  =  3  THEN  B$  =  " ( 1/C) A2=( 1/Co) A2+2*K*t" 

2560  IF  ORDER  =  4  THEN  B$  =  "C=Co- K*t" 

2570  PRINT  B$ 

2580  IF  B  >  0  THEN  LOCATE  3,  8 

2590  IF  B  <  0  THEN  LOCATE  3,  INT(QLNR  *  80  /  640  -  .5)  -  15 

2600  PRINT  "K=";  K 

2610  FLORD  =  0 

2620  QLOCWY  =  (QLOCY3  -  QLOCY2  +  1)  /  QNINPY 

2630  QLNWY  =  (QLNB  -  QLNT)  /  QNINPY 

2640  REM  draw  vertical  grid  and  labels 

2650  FOR  QI  =  0  TO  QNINPY 

2660  IF  QI  =  0  THEN  2690 

2670  LOCATE  INT(QLOCY3  -  QI  *  QLOCWY  +  1),  QLOCX2 

2680  PRINT  USING  Q$ ;  QMINPY  +  QI  *  QWDTY; 

2690  QNLY1  =  QNLB  -  QI  *  QLNWY 

2700  NEXT  QI 

2710  GOSUB  3840 

2720  LINE  (QLNL,  QLNB)-(QLNR,  QLNB) 

2730  QLOCWX  =  (QLOCX4  -  QLOCX3  +  1)  /  QNINPX 

2740  QLNWX  =  (QLNR  -  QLNL)  /  QNINPX 

2750  REM  draw  horizontal  grid  and  labels 

2760  FOR  QI  =  0  TO  QNINPX 

2770  QLOC  =  QLOCY4 

2780  LOCATE  QLOC,  INT(QLOCX3  +  QI  *  QLOCWX  +  .5) 

2790  PRINT  USING  Q$ ;  QMINPX  +  QI  *  QWDTX; 

2800  QNLX1  =  QLNL  +  QI  *  QLNWX 

2810  NEXT  QI 

2820  LINE  (QLNL,  QLNB) -(QLNL,  QLNT) 

2830  GOSUB  3730 

2840  IF  QFLG  =  1  THEN  RETURN 

2850  REM  draw  gx-gy  array 

2860  FOR  QI  =  N  TO  QN  STEP  QRES 

2870  GOSUB  3440 

2880  LINE  (QX2,  QY2)-(QX2,  QY2),  ,  B 

2890  NEXT  QI 

2900  RETURN 

2910  IF  QSCL  >  0  THEN  2930 

2920  PRINT  "NO  GRID  DRAWN,  CHECK  QFLG  PARAMETER":  RETURN 

2930  FOR  QI  =  N  TO  QN  STEP  QRES 

2940  GOSUB  3440 
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2950 
2970 

2980 
3110 
3430 
3440 

3460 
3480 
3490 
3500 
3510 

3520 
3530 
3540 
3550 
3560 
3570 
3580 
3590 
3600 
3610 
3620 
3630 
3640 
3650 
3660 
3670 
3680 
3690 
3700 
3710 

3720 
3730 
3740 
3750 
3760 
3770 
3780 
3790 
3800 
3810 
3820 
3830 
3840 
3850 
3860 
3870 
3880 
3890 


QLNB  OR  QY2  <  QLNT 


QY2)-(QX2,  QY2) 


B 


QLNL) 
QLNB) 


REM  chek  for  points  outside  the  grid 

IF  QX2  <  QLNL  OR  QX2  >  QLNR  OR  QY2  > 

THEN  QFL2  =  1  ELSE  QFL2  =  0 

IF  QFL2  =  0  THEN  LINE  (QX2, 

NEXT  QI 

RETURN 

REM  routine  for  linear  transformation  from  subject  to 

object  space 

QX2  =  FNQSC(QX(QI)  /  QSCLX,  QMAXPX,  QMINPX,  QLNR, 

QY2  =  FNQSC(QY(QI)  /  QSCLY ,  QMAXPY ,  QMINPY ,  QLNT, 

RETURN 

REM  routine  to  find  "nice"  scale  intervals 

IF  QMIN  >=  QMAX  THEN  BEEP:  PRINT  "Array  can  not  be 

plotted":  STOP 

QEPS  =  .025:  QA  =  ABS(QMIN) 

IF  ABS(QMIN)  <  ABS(QMAX)  THEN  QA  -  ABS(QMAX) 

QSCL  =  10  A  INT(LOG(ABS(QA))  /  LOG(10)) 

QMINA  =  QMIN  /  QSCL:  QMAXA  =  QMAX  /  QSCL 

IF  QMIN  =  0  THEN  QMIN  =  -1 

QD  =  (QMAXA  -  QMINA)  /  QNIN:  QJ  =  QD  *  QEPS 

QE  =  INT(LOG(ABS(QD))  /  LOG(10)) 

QF  =  QD  /  10  A  QE:  QV  =  10 

IF  QF  <  SQR(50)  THEN  QV  =  5 

IF  QF  <  SQR(10)  THEN  QV  =  2 

IF  QF  <  SQR(2)  THEN  QV  =  1 

QWDT  -  QV  *  10  A  QE 

QG  =  INT (QMINA  /  QWDT) 


■  QMINA  /  QWDT)  <  QJ  THEN  QG  =  QG  +  1 
QG 
/  QWDT) 


QWDT 
QH 


+  1 
+  1  -  QH)  <  QJ  THEN  QH  =  QH  -  1 


IF  ABS(QG  +  1  - 

QMINP  =  QWDT  * 

QH  =  INT (QMAXA 

IF  ABS( QMAXA  / 

QMAXP  =  QWDT  * 

QNINP  =  QH  -  QG 

IF  ABS(QMAXP)  >=  10  OR  ABS(QMINP)  >=  10  THEN  QSCL  =  QSCL 

10:  GOTO  3550 

RETURN 

REM  routine  to  draw  x  tic  marks 

QLNXB1  =  QLNR 

QDTXP  -  INT ((QLNR  -  QLNL)  /  QNINPX) 

QLNXB1  =  QLNXB1  +  QDTXP 

FOR  QQ  =  1  TO  QNINPX 

QLNXB1  =  QLNXB1  -  QDTXP 

IF  QLNXB1  <  QLNL  THEN  RETURN 

QLNB3  =  QLNB  -  3 

LINE  (QLNXB1,  QLNB)-(QLNXB1 , 

NEXT  QQ 

RETURN 

REM  routine  to  draw  y  tic  marks 

QYLNB1  -  QLNB 

QDTYP  =  INT ((QLNB  -  QLNT) 

FOR  QQ  =  1  TO  QNINPY 

QYLNB1  =  QYLNB1  -  QDTYP 

QLNL4  =  QLNL  +  4 


QLNB 3 ) 


/  QNINPY) 
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3900    LINE    (QLNL,    QYLNB1 ) -(QLNL4 ,    QYLNB1) 

3920    IF   QQ   =   QNINPY   AND   QDTYP    =    1    THEN   QQ   =   QQ    -    34 

3925    NEXT   QQ 

3930    RETURN 

3940    '**********    LINEAR   REGRESSION*************************** 

3950    INPUT    "LEFT   MOST   BORDER    (1-512)";    N 

3960    INPUT    "RIGHT   MOST   BORDER    (1-512)";    J 

3970    QLNB   =    164 

3980    DEFDBL   P,    S 

3990   M   =   J   -    N   +    1 

4000    SX   =    0:    SY   =   0:    PX  =   0 :    PY   =    0 :    PC   =   0 

4010    IF   ORDER   <    1    OR   ORDER    >    4    THEN   ORDER   =    4 

4020    FOR    I    =   N   TO   J 

4030    IF   ORDER   =    4    THEN   QLY(I)    =    CON(I) 

4040    IF   ORDER   -    1    THEN   QLY(I)    =    LOG( ABS(CON( I ) ) ) 

4050    IF   ORDER   -    2    THEN   QLY(I)    =    1    /    CON(I) 

4060  IF   ORDER   =    3    THEN   QLY(I)    =    (1    /    CON(I))     A    2 

4061  RHGRFLAG   =    1 

4070    SX   =   SX   +   QX(I):    SY   =   SY    +   QLY(I) 

4080    PX   =   PX   +   QX(I)    *    QX(I):    PY   =   PY   +   QLY(I)    *    QLY(I) 

4090    PC   =    PC    +   QX(I)    *    QLY(I) 

4100    NEXT    I 

4110    D#    =   M    *    PX    -    SX    *    SX 

4120   A   =    (SY    *    PX   -    PC    *    SX)    /    D# 

4130    B   =    (M    *    PC    -    SX    *    SY)    /    D# 

4140    VX#    =    (PX   -    SX    *    SX    /   M)    /    (M   -    1) 

4150    VY#    =    (PY    -    SY    *    SY    /   M)    /    (M   -    1) 

4160    RR#    =   B    *    B    *    VX#    /   VY# 

4170   R   =   SQR(RR#) 

4180    E   =    SQR((1    -   RR#)    /    (M    -    2))    /   R 

4190   RE   =    (M   -    1)    *    VY#    *    (1    -   RR# ) 

4200    GB   =   ABS(E    *    B) 

4210    GA   =   GB    *    SQR(PX   /   M) 

4220    GP   =    SQR(RE    /    (M   -    1)) 

4230    K   =   ABS(B):    IF   ORDER   =    3    THEN   K   =   K   / 

4240    PRINT    "INTCPT=";     :    PRINT   USING    "##.## 

4250    PRINT    "    (STD   DEV=" ;     :    PRINT   USING    "## 

4260    PRINT    ")       SLOPE=";     :    PRINT   USING    '•##. 

4270    PRINT    "(STD   DEV=" ;     :    PRINT   USING    "##, 

4280    PRINT    "R=";    R;    "         E=" ;    EE;     "  STD 

4290    IF   R   >    .8   GOTO    4320 

4300    PUT    (0,    QLNB   -    64),    BTARY ,    PSET 

4305    LOCATE    INT( ( (QLNB    -    64)    /    200)    *    25    +    1)    +    1,     1 

4310    PRINT   TAB(30);     :    PRINT    "THIS    IS   NOT   A  WELL   FIT   CURVE" 

4  315    REM   FOR   LLKK   =    1    TO    2000 : MY BADDOG=LOG ( LLKK ): NEXT 

4320    IF    INKEY$    =    ""    THEN    4320 

4330    PUT    (0,    QLNB   +    11),    BTARY,    PSET:    LOCATE    INT((QLNB    /    200) 

25    +    1)    +    2,    1 
4340    INPUT    "WOULD   YOU   LIKE   A   PLOT   OF   THE   RESULTS";    B$ 
4350    B$    =   LEFT$(B$,    1) 
4360    IF    B$    =    "N"    THEN   RETURN 
4370    PUT    (0,    QLNB   +    11),    BTARY,    PSET:    LOCATE    INT((QLNB   /    200) 

25    +    1)    +    2,    1 


2 

AAA/* 

II   •        A    • 

-##' 

\    A    A    A   If     . 

/ 

GA; 

##A' 

.  A  A  II    • 
/ 

B; 

##A' 

^  A   A  ||    • 
t 

GB ; : PRINT 

") 

DEV 

PTS=r 

1 ;   gp; 
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4380 
4390 
4400 

4405 
4410 
4420 
4430 
4440 
4450 
4460 
4465 
4470 
4480 
4490 
4500 
4510 
4520 
4530 
4540 
4550 
4560 
4570 
4580 
4590 
4600 
4610 
4620 
4630 
4640 
4650 
4660 
4665 
4670 
4680 
4690 
4695 
4700 
4710 
4720 
4730 
4740 
4745 
4750 
4760 
4765 
4770 
4780 
4790 
4800 
4810 
4820 
4821 


OR 

J 


IF   B$    <>    "Y 
FLORD   =    1 
INPUT    "WHAT 
QRES 

IF   QRES    <    1 
FOR    I    =   N   TO 
CC   -   QY(I) 
QY(I)    =   QLY(I) 
QLY(I)    =   CC 
NEXT    I 

QFLG  =  0:  QMINX  = 
QMAXY  =  QY(J):  QN 
QYN  =  B  *  QX(N)  + 
QYJ  =  B  *  QX(J)  + 
GOTO  4520 
QY(N)  THEN 
QY(J)  THEN 
GOTO  4550 
QY(N)  THEN 
QY(J)  THEN 
2160 

FNQSC(QYN  / 
FNQSC(QX(N) 
FNQSC(QYJ  / 
FNQSC(QX(J) 


THEN   GOTO    4  3  30 
RESOLUTION   WOULD    YOU   LIKE   TO   GRAPH    (1    TO    10)" 


QRES    >    20    THEN   QRES    =    1 


QX(N) 
=   J 
A 
A 


QMAXX   =   QX(J):    QMINY   =   QY(N) 


IF 
IF 
IF 
IF 
IF 
IF 


B  < 
QYN 
QYJ 
B  > 
QYN 
QYJ 


0 

< 
> 
0 

> 
< 


QMINY 
QMAXY 

QMAXY 
QMINY 


=   QYN 
=   QYJ 

-  QYN 

-  QYJ 


ELSE   QMINY 
ELSE   QMAXY 


ELSE   QMAXY 
ELSE   QMINY 


QY(N) 
QY(J) 

QY(N) 
QY(J) 


GOSUB 

QYN   = 

QXN   = 

QY2    = 

QX2    = 

LINE    (QXN,    QYN)-(QX2, 

FOR    I    =   N   TO   J 

CC   =   QY(I) 

QY(I)    =   QLY(I) 

QLY(I)    =   CC 

NEXT    I 

QMINX   =   QX(1):    QMAXX   =   QX(512) 

QMAXY   =   QY(514):    QN   =    512 

FLLOOP   =    1    +    FLLOOP 

IF    FLLOOP    >    1    THEN   GOTO    4765    ELSE 

PUT    (0,    QLNB   +    11) 

LOCATE    I NT ((QLNB   / 

PRINT    "INTCPT=";     : 

"    (STD   DEV=" 

"(    SLOPE="; 

"    (STD   DEV=" 

"R=" •   R ■   " 

1   TO 

THEN 


QSCLY,    QMAXPY,    QMINPY,    QLNT,    QLNB) 
/   QSCLX,    QMAXPX,    QMINPX,    QLNR,    QLNL) 
QSCLY,    QMAXPY,    QMINPY,    QLNT,    QLNB) 
/   QSCLX,    QMAXPX,    QMINPX,    QLNR,    QLNL) 
QY2) 


QMINY   =   QY(515) 


GOSUB    1570 


+ 


2,    1 


BTARY,    PSET 
200)    *    25    +    1) 
PRINT   USING    "# 

:    PRINT   USING 

PRINT   USING    »##.##*      *"; 

:PRINT   USING    "##.##AAAA" 


"## 


A; 

ff AAAAH.   GA; 

B; 
GB ; : PRINT 


E=" ;    EE;    "  STD   DEV   PTS=" ;    GP; 

2000:MYBADDOG=LOG(LLKK) :NEXT 
GOTO    4750 


PRINT 

PRINT 

PRINT    "    (STD   DEV=";     : PRINT   USING    "##.##"      A";    GB;: PRINT    ■)■ 

PRINT 

REM   FOR   LLKK   = 

IF    INKEY$    =    "" 

FLLOOP   =    0 

REM   GOSUB    1570 

RETURN 

'*************    SAVE    ************************************** 

CLOSE    #1 

IF    FLDATA   >    0    THEN   GOTO    4820 

PRINT    "NO   DATA   SOURCE.    RETURNING   TO   MENU":    GOTO    4910 

IF   FLDATA   =    2    THEN   PRINT    "DATA   PREVIOUSLY   SAVED" 

INPUT    "SAVE   AGAIN   UNDER   ANOTHER   NAME?    "?    RJHANS$ 
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4822  IF  LEFT$(RJHANS$,  1)  =  "N"  THEN  GOTO  4900 

4823  INPUT  "WHICH  DRIVE  ■ ;  DRV$ :  IF  DRV$  =  ■■  THEN  DRV$  =  "C" 

4824  DRV$  =  LEFT$(DRV$,  1) 

4825  IF  DRV$  =  "A"  OR  DRV$  -  "B"  GOTO  4829 

4826  IF  DRV$  =  "C"  OR  DRV$  =  "D"  GOTO  4829 

4827  GOTO  4823 

4829  DRV$  =  DRV$  +  "S" 

4830  INPUT  "FILENAME  " ;  FILENAME$ 
4840  B$  =  DRV$  +  LEFT$ ( FILENAME$ ,  8) 
4850  OPEN  "O",  #1,  B$ 

4860  SIZE  =  518 

4870  FOR  I  =  1  TO  SIZE  STEP  5 

4880  PRINT  §1,    QY(I),  QY(I  +  1),  QY(I  +  2),QY(I  +  3),QY(I  +  4) 

4890  NEXT  I 

4900  CLOSE  #1 

4910  OPEN  "coml:9600,n,8,l,cs,ds,cd"  FOR  RANDOM  AS  #1 

4920  FLDATA  =  2 

4930  RETURN 

4940  '***************  LOAD  *********************************** 

4950  CLOSE  #1 

4951  INPUT  "WHICH  DRIVE  ";  DRV$ :  IF  DRV$  =  ""  THEN  DRV$  =  "CM 

4952  DRV$  =  LEFT$(DRV$,  1) 

4953  IF  DRV$  =  "A"  OR  DRV$  =  "B"  GOTO  4959 

4954  IF  DRV$  =  "C"  OR  DRV$  =  "D"  GOTO  4959 

4955  GOTO  4951 

4959  DRV$  =  DRV$  +  " : " 

4960  INPUT  "FILENAME";  FILENAME$ 
4970  B$  =  DRV$  +  LEFT$ ( FILENAME$ ,  8) 
4980  OPEN  "I",  #1,  B$ 

4982  ON  ERROR  GOTO  5230 

4990  SIZE  =  518 

5000  FOR  I  =  1  TO  SIZE  STEP  5 

5010  INPUT  #1,  QY(I),  QY(I  +  1),  QY(I  +  2),QY(I  +  3),QY(I  +  4) 

5020  NEXT  I 

5030  PER  =  QY(513):  QMAXY  =  QY(514):  QMINY  =  QY(515) 

5035  VO  =  QY(516):  VM  =  QY(517):  L  =  QY(518):  EE  =  QY(519) 

5040  PRINT  "FILE  ";  B$ ;  "  FOUND." 

5050  CLOSE  #1 

5060  OPEN  "coml:9600,n,8,l,cs,ds,cd"  FOR  RANDOM  AS  #1 

5061  LEFF  =  2.303  *  L  *  EE:  VFULL  =  VM  -  VO 

5062  PRINT  "VO,VM/VFULL/L/LEFF,E":  PRINT  VO ,  VM,  VFULL,  L;  " 
";  LEFF,  EE 

5063  INPUT  "WANT  TO  CHANGE  V0,VM,L,OR  E  " ;  ANANS$ 

5064  IF  ANANS$  =  "Y"  THEN  GOSUB  8000 

5068  LEFF  =  2.303  *  L  *  EE:  VFULL  =  VM  -  VO 

5070  REM  FOR  ZSAM  =  1  TO  5000:  ZSAM  =  LOG(ZSAM):  NEXT 

5085  FOR  I  =  1  TO  512 

5088  QX(I)  =  I  *  PER 

5090  CON(I)  =  -LOG(ABS((QY(I)  -  VO )  /  VFULL))  /  LEFF 

5092  OD(I)  =  -LOG(ABS((QY(I)  -  VO )  /  VFULL))  /  2.303 

5095  REM  PRINT  CON(I) 

5100  NEXT  I 

5110  FLDATA  =  2 
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5115 
5120 
5130 
5140 
5150 
5160 
5170 
5180 
5190 
5200 
5210 
5220 
5230 
5240 
5250 
5255 
5256 
5260 
5270 
5280 
5290 
5300 
5310 
5320 
5330 
5340 
5350 
5360 
5370 
5380 
5390 
5400 
5410 
5420 
5430 
5440 
5450 
5460 
5470 
5480 
5490 
5500 
5510 
5520 
5530 
5540 
5550 
5560 
5570 
5580 
5585 
5590 
5600 


RHGRFLAG  =  1 

RETURN 

'********•*****•  N'th 


PRINT  :  PRINT  "CHOOSE 

PRINT  "1)  FIRST  ORDER 

PRINT  "2)  SECOND  ORDER 

PRINT  "3)  THIRD  ORDER 

PRINT  "4)  ABSORBANCE  ONLY 

INPUT  "WHICH  NUMBER" ;  ORDER 

IF  ORDER  <  1  OR  ORDER  >  4  GOTO 


ORDER  ***************************** 
THE  ORDER  OF  THE  REACTION":  PRINT 

ln(C)=ln(Co)-K*t" 

(l/C)=(l/Co)+K*t" 

(l/C)A2=(l/Co)A2+2*K*t" 

C=Co-K*t":  PRINT 


5140 


'LOAD 

' PLOTTER 


RETURN 

***********************  ERROR  *************************** 

'IF(57=ERR)AND(180  >  ERL)  THEN  RESUME  160 

IF  57  =  ERR  THEN  PRINT"DP600  OR  DISK  I/O  ERROR" :RESUME  140 

IF  53  =  ERR  GOTO  5290 

IF  4970  =  ERL  THEN  GOTO  5290 

IF  4980  =  ERL  THEN  GOTO  5290 

REM  PRINT  "STOP":  PRINT  ERR,  ERL 

REM  IF  ERR=6  THEN  BAD=ERL:GOTO  ERL 

RESUME  NEXT 

PRINT  "FILE  ";  B$ ;  "  NOT  FOUND" 

FILENAMES  =  "■ 

CLOSE  #1 

RESUME  4950 

*******************  COMPARE  ***************************** 

INPUT  "HOW  MANY  PLOTS  DO  YOU  WISH  TO  MAKE  (1  TO  3)";  NOPLO 

IF  NOPLO  <  1  OR  NOPLO  >  3  GOTO  5340 

PRINT  :  PRINT  "LOAD  FIRST  DATA  SET  FROM  DISK":  PRINT 

GOSUB  4940 

GOSUB  1480 

Nl  =  N:  Jl  =  J 

DIM  QY1(515) 

FOR  I  =  N  TO  J 

QY1(I)  =  QLY(I) 

NEXT  I 

IF  B  <  0  GOTO  5460 

QY1(515)  =  QY1(N):  QY1(514) 

QY1(515)  =  QY1(J):  QY1(514) 

IF  NOPLO  =  1  GOTO  5750 

PRINT  :  PRINT  "LOAD  SECOND  DATA 

GOSUB  4940 

GOSUB  1480 

N2  =  N:  J2  =  J 

DIM  QY2(J) 

FOR  I  =  N  TO  J 

QY2(I)  =  QLY(I) 

NEXT  I 

IF  B  <  0  GOTO  5590 

IF  QY2(N)  <  QY1(515)  THEN 

IF  QY2(J)  >  QY1(514)  THEN 

GOTO  5610 

IF  QY2(N)  >  QY1(514)  THEN 

IF  QY2(J)  <  QY1(515)  THEN 


=  QY1(J) :  GOTO  5470 
-  QY1(N) 


SET  FROM  DISK":  PRINT 


'LOAD 

' PLOTTER 


QY1(515) 
QY1(514) 

QY1(514) 
QY1(515) 


QY2(N) 
QY2(J) 

QY2(N) 
QY2(J) 
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5610  IF  NOPLO  =  2  GOTO  5750 

5620  PRINT  :  PRINT  "LOAD  THIRD  DATA  SET  FROM  DISK":  PRINT 

5630  GOSUB  4940  'LOAD 

5640  GOSUB  1480  'PLOTTER 

5650  N3  =  N:  J3  =  J 

5660  DIM  QY3(J) 

5670  FOR  I  =  N  TO  J 

5680  QY3(I)  =  QLY(I) 

5690  NEXT  I 

5700  IF  B  <  0  GOTO  5730 

5710  IF  QY3(N)  <  QY1(515)  THEN  QY1(515)  =  QY3(N) 

5720  IF  QY3(J)  >  QY1(514)  THEN  QY1(514)  =  QY3(J) 

5725  GOTO  5750 

5730  IF  QY3(N)  >  QY1(514)  THEN  QY1(514)  -  QY3(N) 

5740  IF  QY3(J)  <  QY1(515)  THEN  QY1(515)  =  QY3(J) 

5750  N  =  Nl:  J  =  Jl 

5760  IF  N2  <  Nl  AND  NOPLO  >  1  THEN  N  =  N2 

5770  IF  J2  >  Jl  AND  NOPLO  >  1  THEN  J  -  J2 

5780  IF  N3  <  N  AND  NOPLO  >  2  THEN  N  =  N3 

5790  IF  J3  >  J  AND  NOPLO  >  2  THEN  J  =  J3 

5800  QFLG  =  0:  FLORD  =  0:  QMINX  =  QX(N):  QMAXX  =  QX(J) 

5805  QMINY  =  QY1(515):  QMAXY  =  QY1(514):  QN  -  Jl :  N  =  Nl 

5810  PRINT  :  INPUT  "WHAT  RESOLUTION  WOULD  YOU  LIKE  TO  GRAPH  (1 

TO  10)";  QRES 
5820  FOR  I  =  Nl  TO  Jl 
5830  QY(I)  =  QY1(I) 
5840  NEXT  I 
5850  ERASE  QY1 
5860  FLORD  =  0 
5870  GOSUB  2170 
5880  LOCATE  1,  59 

5890  PRINT  ■  " 

5900  IF  NOPLO  =  1  GOTO  6040 
5910  QFLG  =2:  QN  =  J2 :  N  =  N2 
5920  FOR  I  =  N2  TO  J2 
5930  QY(I)  =  QY2(I) 
5940  NEXT  I 
5950  ERASE  QY2 
5960  GOSUB  2170 
5970  IF  NOPLO  =  2  GOTO  6040 
5980  QN  =  J3:  N  =  N3 
5990  FOR  I  =  N3  TO  J3 
6000  QY(I)  =  QY3(I) 
6010  NEXT  I 
6020  GOSUB  2170 
6030  ERASE  QY3 
6040  PUT  (0,  QLNB  +11),  BTARY ,  PSET:  LOCATE  INT((QLNB  /  200)  < 

25  +  1)  +  2,  1 
6050  PRINT  "YOU  WILL  NEED  A  NEW  DATA  SOURCE  AFTER  RETURNING  TO 

MENU" 
6060  IF  INKEY$  =  ""  THEN  GOTO  6060 
6070  PUT  (0,  QLNB  +  11),  BTARY,  PSET:  LOCATE  INT ( (QLNB  /  200)  ' 

25  +  1 |  +  2,  1 
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6080  PRINT  "PRESS  ANY  KEY  ONCE  TO  GET  A  SCREEN  DUMP" 

6085  PRINT  "TWICE  TO  RETURN  TO  MENU"; 

6090  IF  INKEY$  =  ""  THEN  GOTO  6090 

6100  PUT  (0,  QLNB  +  11),  BTARY ,  PSET:  LOCATE  INT( (QLNB  /  200) 

25  +  1)  +  2,  1 
6110  IF  INKEY$  =  ""  THEN  GOTO  6110 
6120  FLDATA  =  0 
6130  RETURN 

7000  REM  SUBROUTINE  GET  INFO 
7040  ' 
7070  CLS 

7080  REM  OPEN  "coml : 9600 ,n, 8 , 1 ,cs ,ds ,cd"  AS  #1 
7100  BOUT$  =  B$ 
7110  PRINT  #1,  BOUT$ 

7120  FOR  I  =  1  TO  1500:  QWERT1  =  QWERT1  +  1 
7140  NEXT  I 
7150  J  =  1 

7160  V  -  LOC(l):  IF  V  <  1  THEN  PRINT  "NO  MESSAGE":  GOTO  7350 
7170  FOR  ZQQZ  =  1  TO  500:  FF  =  ZQQZ  *  ZQQZ :  NEXT 
7180  AGI$(J)  =  INPUT$(V,  #1) 
7190  ON  ERROR  GOTO  5220 

7200  V  =  LOC(l):  IF  V  <  1  THEN  GOTO  7230 
7210  J  =  J  +  1 
7220  GOTO  7160 
7230  FOR  HH  =  1  TO  J 
7240  BGI$(HH)  =  "" 

7250  REM  PRINT  AGI$(HH),  VAL( AGI$ (HH) ) 
7260  FOR  UU  =  1  TO  LEN( AGI$ (HH) ) 
7270  QAQ$  =  MID$(AGI$(HH),  1) 
7280  ZORP$  =  BGI$(HH) 

7285  IF  ASC(QAQ$)  =  45  THEN  GOTO  7300 

7290  IF  ASC(QAQ$)  <  48  OR  ASC(QAQ$)  >  91  THEN  GOTO  7330 
7300  BGI$(HH)  =  BGI$(BB)  +  QAQ$ 
7310  NEXT  UU 
7320  NEXT  HH 
7330  MYINFO  =  VAL(ZORP$) 
7350  RETURN 

8000  REM  CHANGE  V0,VMAX  ********************************** 
8010  INPUT  "DO  YOU  WANT  TO  CHECK  ACTUAL  VOLTAGE  VALUES  OF 

DATASET";  ANANS$ 
8020  IF  LEFT$(ANANS$,  1)  =  "N"  THEN  GOTO  8090 
8030  INPUT  "FIRST  ,  LAST  DATA  TO  USE  ";  POINT1 ,  POINT2 
8035  TEMPSUM  =  0 

8040  FOR  KZ  -  POINT1  TO  POINT2 
8050  TEMPSUM  =  TEMPSUM  +  QY(KZ) 
8060  NEXT  KZ 

8070  VNOW  =  TEMPSUM  /  (POINT2  -  POINT1  +  1) 
8080  PRINT  "AVERAGE  VOLTAGE  FOR  THAT  RANGE  IS  " ;  VNOW 
8090  INPUT  "CHANGE  VO   "j  ANANS$ 
8100  IF  ANANS$  =  "N"  THEN  GOTO  8200 
8110  INPUT  "NEW  VO";  VO 

8200  INPUT  "CHECK  VOLTAGES  AGAIN  ";  ANANS$ 
8210  IF  LEFT$(ANANS$,  1)  =  "Y"  THEN  GOTO  8030 
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8240  INPUT  "CHANGE  VM  " ;  ANANS$ 

8250  IF  ANANS$  =  "N"  THEN  GOTO  8300 

8260  INPUT  "NEW  VMAX  ";  VM 

8270  VFULL  =  VM  -  V0 

8300  PRINT  "L  IS  NOW  ";  L:  INPUT  "CHANGE  L   " ;  ANANS$ 

8320  IF  ANANS$  =  "N"  THEN  GOTO  8400 

8330  INPUT  "NEW  L" ;  L 

8400  PRINT  "E  IS  NOW  « ;  EE:  INPUT  "CHANGE  E   ";  ANANS$ 

8420  IF  ANANS$  =  "N"  THEN  GOTO  8500 

8430  INPUT  "NEW  E  ";  EE 

8500  LEFF  =  2.303  *  L  *  EE:  VFULL  =  VM  -  VO 

8510  QY(516)  =  VO:  QY(517)  =  VM:  QY(518)  =  L:  QY(519)  -  EE 

8520  PRINT  "VO, VM, VFULL, L, LEFF, E":  PRINT  VO ,  VM,  VFULL,  L;  " 

■ ;  LEFF,  EE 
8515  IF  VFULL  =  0  OR  LEFF  =  0  THEN  PRINT  "BAD  CONSTS .  " :  GOTO 

8010 
8530  FOR  I  =  1  TO  512 
8540  QX(I)  =  I  *  PER 

8550  CON(I)  =  -LOG(ABS( (QY(I)  -  VO )  /  VFULL))  /  LEFF 
8560  REM  PRINT  CON(I) 
8570  NEXT  I 
8580  FLDATA  =  2 
8600  RETURN 

9800  REM  *****************DISK  DIR*********************** 
9810  CLOSE  #1 

9820  INPUT  "WHICH  DRIVE  ■ ;  DRV$ :  IF  DRV$  =  ""  THEN  DRV$  =  "C" 
9830  DRV$  =  LEFT$(DRV$,  1) 

9840  IF  DRV$  =  "A"  OR  DRV$  =  "B"  GOTO  9870 
9850  IF  DRV$  =  "C"  OR  DRV$  =  "D"  GOTO  9870 
9860  GOTO  9820 
9870  MYCALL$  =  DRV$  +  " : " 
9875  FILES  MYCALL$ 

9880  OPEN  "coml:9600,n,8,l,cs,ds,cd"  FOR  RANDOM  AS  #1 
9890  RETURN 

9892  INPUT  "WHICH  DRIVE";  DRV$ :  IF  DRV$  =  "  "  THEN  DRV$  =  "C" 
9894  DRV$  =  LEFT$(DRV$,  1) 

9896  IF  DRV$  =  "A"  OR  DRV$  =  "B"  GOTO  9905 
9898  IF  DRV$  =  "C"  OR  DRV$  =  "D"  GOTO  9905 
9900  GOTO  9892 
9905  DRV$  =  DRV$  +  " : " 
9910  INPUT  "FILE  NAME";  FILENAME$ 
9915  B$  =  DRV$  +  LEFT$(FILENAME$,  8) 
9920  OPEN  "0",  #1,  B$ 
9925  FOR  I  =  N  TO  J  STEP  5 
9930  PRINT  #1,  QX(I),  ",",  OD(I),  " , " 
9940  NEXT  I 

9945  CLOSE  #1 

9946  GOTO  2106 
9999  END 
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